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SUMMARY 


Tlie  inia^iu'^  properties  of  parametric  upconversion  have  been 
investigated  analytically  and  experimentally  for  the  case  of  monochromatic 
object  waves,  and  photograplis  of  upconverted  10.6  jum  images  obtained. 

Analyses:  Using  geometric  optics  an  expression  was  derived 
relating  the  locations  of  the  object,  the  pump  and  the  upconverted  image. 

It  was  found  that  in  general  the  image  quality  is  degraded  by  a  thickness 
aberration  resulting  from  sum-frequency  contributions  from  each  differen¬ 
tial  slab  of  the  nonlinear  interaction  material.  The  magnitude  of  this  first 
order  thickness  aberration  was  found  to  approach  zero  for  the  case  of  plane 
wave  pump  interacting  with  plane  wave  object  waves.  The  angular  aperture 
of  the  upconverter  was  found  'o  be  determined  by  the  phase  match  conditions. 
The  effect  of  higher  order  factors  such  as  crystal  birefringence,  phase 
mismatch,  and  pump  wave  divergence  were  considered  qualitatively. 

The  analysis  was  then  extended  by  applying  the  techniques  of 
physical  optics.  Using  Fourier  transformer  formalism  an  expression  was 
derived  relating  the  angular  spectrum  of  the  sum -frequency  field  with  that 
of  the  object  and  pump  fields.  Phase  mismatch,  crystal  birefringence,  and 
the  effects  of  a  finite  aperture  were  included  in  this  analysis.  Two  cases 
were  examined  in  detail,  one  using  a  plane  wave  pump  and  the  other  with  a 
Gaussian  distributed  pump  beam.  The  Gaussian  pump  case  was  found  to 
suffer  more  image  degradation  than  the  plane  wave  case  which  approached 
tlie  diffraction  limit. 


A  characteristic  equation  was  developed  to  describe  the  overall 
performance  of  the  image  upconverter.  Object  scene  radiance,  image 
detail,  exposure  time,  quantum  conversion  efficiency,  detecto’  sensitivity 
and  the  interaction  aperture  are  related  to  a  universal  constant. 

Experimental  Results:  A  series  of  image  upconversion  experi¬ 
ments  were  carried  out  with  a  10.6-|im  illuminated  object,  a  1.06-um  pump 
and  an  oriented  proustite  crystal  (Ag^AsS^).  The  line  spread  functions  were 
measured  for  two  upconversion  configurations:  plane  wave  interactions  and 
the  Gaussian  shaped  pump  beam.  Measured  results  on  image  location,  aber¬ 
rations,  the  field  of  view,  and  angular  resolution  agree  reasonably  well  with 
the  theoretical  results.  For  the  plane  wave  pump  case  about  70  x  20  resolu¬ 
tion  elements  were  measured.  Two  dimensional  images  of  a  resolution 
test  chart  have  been  recorded  on  photographic  film. 

Conclusions:  As  a  result  of  this  analytical  and  experimental  in¬ 
vestigation  it  is  concluded  that  the  monochromatic  image  upconversion 
process  can  approach  diffraction  limited  operation  if  the  upconversion  pro¬ 
cess  takes  place  with  plane  wave  object  and  pump  beams  (Fourier  space 
interaction). 

A  parametric  study  has  been  done  for  such  imaging  upconverter 
characteristics  as  conversion  efficiency,  number  of  resolution  elements, 
pump  power,  pump  power  density,  crystal  figure  of  merit,  and  overall  sen¬ 
sitivity:  This  approach  indicates  some  present  limitations  of  the  upconver¬ 
sion  process  and  the  future  possibilities  that  can  be  based  on  the  extrapola¬ 
tion  and  trends  in  the  present  technology. 

Publications :  1.  W.  Chiou,  "Geometric  Optics  Theory  of 
Parametric  Image  Upconversion,  "  Journal  of  Applied  Physics,  Vol  42, 
p  1985,  March  1971. 

2.  W.  Chiou  and  F.  Pace,  "Image  Upconver?  ion  of  Infrared 
10.6-^m  Radiation,"  lEEE/OSA  Conference  on  Laser  Engineering  and 
Applications,  Washington,  D.C,  June  1971. 
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This  is  the  final  rf|K»rl  under  Cunirari  NOOOH-TU-t  -oHil  h.ivnir. 

as  its  objective  the  deinonstratiim  and  o|)erahon  a  10.C>uni  ini.n:«  innvrr- 

ter  system  with  as  high  a  conversion  efficiency  as  |Hisstltle.  :i  ri’MiUiiiim  nl 

- 12 

150  y  150  resolution  elements  and  noise  equivalent  |miW(  r  (NKlM  of  10 
watt/HZ^^^. 

The  (iarametric  imace  uiiconverter  represents  a  new  cla^s  of 
infrared  (IR)  device  technology.  It  combines  ihc  following  characteristics 
hitherto  unattainable: 


•  Permits  detection  of  10. 6-um  radiation  by  conventional 
image  tube  techniques 

•  Obtains  imaging  in  real  time  (no  mechanical  (canning  as 
in  present  FLIR's) 

•  Inherently  low  noise 

•  Operates  at  or  near  room  temperature 

•  Can  be  designed  for  any  wavelength  in  the  1  to  13  ani 
region 

•  Coherent  (upconverter  can  preserve  frequency  and  phase) 

In  the  IR  parametric  upconverter  technique,  a  strong  laser  pump 

at  CO  is  employed  to  upconvert  IR  signal  radiation  u:,o  to  the  sum -frequency 
p  In 

oOs  =  +  ‘*’’p  visible  or  near  visible  region  of  the  spectrum.  The 

frequency  translation  is  achieved  by  propagating  three  waves  (signal,  pump, 
and  sum -frequency)  nearly  colUnearly  and  phase -matched  in  a  crystal  (such 
as  proustite)  having  a  large  value  of  nonlinear  susceptibility.  The  upconver- 
ted  signal  can  then  be  detected  by  a  photoemissive  device  such  as  a  photo¬ 
multiplier,  image  orthicon,  image  intensifier,  or  silicon  \idicon.  The  use 
of  an  Nd:YAG  laser  is  contemplated  as  a  pump  source  due  to  its  CW,  pulsed 
high -power  output  capability,  and  mode  purity.  Moreover,  higher  quantum 
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eflicieiicy  photocatliodes  and  silicon  vidicon  are  presently  under  develop¬ 
ment  at  the  upconvertcr  output  wavelength. 

A.  PROGRAM  APPROACH 


The  approach  chosen  to  achieve  the  program  objectives  includes 
the  following  tasks: 

•  Obtain  analytical  expressions  for  image  upconversion 
analogous  to  the  thick  lens  formula  by  means  of  a 
modified  geometrical  optics  technique 

•  Use  physical  optics  techniques  to  predict  image  up- 
converter  line  and  point  spread  functions 

•  Measure  the  line  spread  function  of  the  upconverter, 
at  various  positions  in  the  fieM  of  view.  Investigate 
effect  of  Gaussian  distribution  of  pump  beam  to  the 
image  resolution. 

•  Investigate  and  compare  the  characteristics  of  two 
generic  upconverter  systems,  the  one  with  the  non¬ 
linear  material  in  the  image  space  and  the  other  with 
the  nonlinear  material  in  the  Fourier  space.  The 
objective  is  to  simultaneously  optimize  conversion 
efficiency  and  image  quality. 

•  Incorporate  all  of  the  previous  findings  into  a  labo¬ 
ratory  device  that  demonstrates  an  image  with  quality 
approaching  the  150  x  150  resolution  element  program 
objective 

•  Measure  and  optimize  image  quality,  conversion 
efficiency,  and  NEP  in  the  desmonstration  setup 

B.  ANALYTICAL  RESULTS 


Qualitative  descriptions  of  parametric  image  upconversion  are 
currently  available  in  literature.  However,  very  little  has  been  published 
on  the  quantitative  considerations  necessary  for  selecticn  and  evaluation  of 
image  upconverter  systems  ix>ssessing  acceptable  image  and  sensitivity 
characteristics. 
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Analyses  and  experiments,  designed  to  fill  this  gap,  were 
carried  out. 

Imaging  properties  of  a  parametric  image  upconversion  process 

are  investigated  both  theoretically  and  experimentally  in  the  present  study. 

The  theoretical  treatment  given  here  is  more  rigorous  than  previous 
(13  14  15) 

work'  ’  ’  The  particular  case  of  tne  nonlinear  interaction  of  a  mono¬ 

chromatic  object  wave  and  a  monochromatic  pump  wave  is  considered.  The 
theory  can  easily  be  extended  to  the  case  of  the  polychromatic  object  waves 
if  one  considers  a  monochromatic  object  wave  as  a  Fourier  component  of 
polychromatic  object  radiation. 

The  geometric  optics  theory  of  the  parametric  image  upconversion 
is  presented  in  Chapter  2.  A  ray  tracing  technique  is  used  in  the  analysis. 

If  the  phase  relation  of  object  waves  is  not  conserved  by  the  parametric  pro¬ 
cess,  the  ray  tracing  technique  is  not  valid.  The  coherent  nature  of  the 
process,  necessary  for  the  effective  interaction  of  the  three  waves  (signal, 
pump,  and  sum  frequency),  transfers  the  phase  of  the  object  waves  to  the 
sum -frequency  waves.  This  consideration  qualitatively  justifies  the  geo¬ 
metric  optics  treatment  on  the  parametric  image  upconversion  process.  A 
more  rigorous  justification  is  given  below.  A  formula  relating  the  locations 
of  a  point  object,  a  point  pump  source,  and  an  upconverted  point  image  is 
obtained  for  a  dispersive  infinitesimal  slab  of  nonlinear  material.  It  is 
assumed  in  the  analysis  that  all  waves  are  phase  matched  and  are  paraxial. 

An  image  produced  by  a  finite  thickness  nonlinear  crystal  is  next  considered. 

A  finite  thickness  crystal  gives  rise  to  thickness  aberrations  because  the 
upconverted  image  position  is  different  for  the  various  differential  slabs  of 
the  crystal.  Thickness  aberrations  are  inherent  in  the  process  and  differ 
from  lens  aberrations  in  ordinary  lenses.  The  amount  of  aberrations  for 
several  different  upconversion  systems  is  analyzed.  It  is  found  that  the  con¬ 
verter  pumped  by  a  plane  wave  pump  and  interacting  with  a  plane  object 
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wave  is  free  of  thickness  aberrations.  The  angular  aperture  of  the  upcon- 
verter  is  shown  to  be  determined  by  the  phase-match  condition.  The  effect 
of  diverging  pump  beam  on  image  quality  is  also  considered.  The  effect  of 
other  less  important  factors,  such  as  crystal  birefringence  and  phase  mis¬ 
match  is  considered  qualitatively. 

The  parametric  image  upconversion  process  is  analyzed  in  Chapter  3 
by  the  physical  optics  technique.  The  purpose  of  this  analysis  is  twofold; 
the  first  is  to  obtain  an  expressio.’  for  the  upconverted  image  field  in  terms 
of  the  infrared  object  field,  the  second  is  to  confirm  the  ray  tracing  tech¬ 
niques  used  in  the  geometric  optics  theory.  The  Fourier  transform  form¬ 
alism  is  adopted  here  because  of  the  complexity  of  the  object  wave  and  the 
image  wave.  A  wave  equation  for  the  angular  spectrum  (Fourier  component) 
of  field  disturbances  with  a  source  term  is  obtained.  The  angular  spectrum 
of  non-linear  polarization  wave  is  considered  to  be  the  source  term  in  the 

wave  equation.  The  wave  equation  is  solved  using  a  technique  developed  by 
(23) 

Kleinman  .  Phase  mismatch  and  crystal  birefringence  are  included  in 
this  treatment.  The  solution  relates  the  angular  spectrum  of  sum-frequency 
field  to  that  of  infrared  object  field.  An  equation  relating  the  positions  of 
upconverted  image,  object,  and  pump  is  obtained  by  using  stationary  phase 
condition.  The  equation  is  identical  with  that  derived  in  the  geometric  optics 
theory  if  crystal  birefringence  is  negligibly  small.  This  result  clarifies  the 
range  of  validity  of  geometric  optic  theory.  The  analysis  also  confirms  that 
aberrations  exist  for  finite  interaction  lengths.  The  result  agrees  with  the 
thickness  aberration  predicted  in  the  geometric  optics  analysis  as  described 
above.  Several  different  image  upconverter  systems  are  analyzed.  It  is 
shown  that  the  amount  of  aberration  for  a  plane  wave  pumped  upconverter  is 
independent  of  the  object  position  if  object  wave  is  not  a  plane  wave.  This 
result  agrees  with  that  of  the  geometric  optics  treatment.  The  special  case 
of  Gaussian  distribution  beam  pumped  upconverter  is  considered.  It  is 
found  that  Gaussian  beam  pump  can,  under  certain  condition,  introduces 
additional  degradation  in  the  upconverted  image. 
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The  characteristic  equation  which  describes  overall  system  per¬ 
formance  of  an  image  upconverter  is  considered  in  Chapter  4.  This  equation 
relates  infrared  object  radiance,  the  amount  of  image  detail,  exposure  time, 
quantum  conversion  efficiency,  detector  sensitivity,  and  the  interaction 
aperture  to  a  characteristic  constant  hi^^,  where  h  is  Planck  constant  and 
is  infrared  frequency  of  monochromatic  object  radiation.  Characterization 
of  an  image  upconverter  enables  one  to  understand  the  limiting  performance 
of  a  particular  system. 

C.  EXPERIMENTAL  RESULTS 

Two  image  upconverter  configurations  have  been  considered  for 
imaging  experiments.  Image  properties  of  two  upconverter  configurations, 
namely: 

•  Plane  wave  object  and  plane  wave  pump 

•  Sperhical  wave  object  and  plane  wave  pump  have 

been  measured 

The  first  configuration  was  chosen  because  of  its  low  thickness 
aberrations.  The  second  configuration  was  chosen  because,  although  it  has 
a  larger  value  of  aberration,  it  provides  a  check  on  this  theory 

Image  properties  were  evaluated  by  measuring  the  line  spread  function 
first  rather  than  recording  an  image  of  a  resolution  test  pattern.  The  line 
spread  function  technique  uses  simple  one -dimensional  measurement  tech¬ 
niques  and  the  data  obtained  can  provide  a  description  of  the  system  in  the 
formalism  of  optical  transfer  functions  employed  in  optical  image  resolution 
analyses. 

Measured  data  on  image  location,  image  size,  and  number  of 
resolution  elements  show  reasonable  agreement  with  theoretically  calculated 
values.  The  field  of  view  of  the  experimental  system  was  12  degrees.  The 
number  of  resolvable  elements  contained  in  this  field  of  view  was  70  (one 
eleme;it/3  milliradians) . 
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Upconverted  mages  at  0.967  ium  of  a  10.6  jum  illuminated  reso¬ 
lution  test  chart  are  then  quantitatively  measured  with  the  scanning  technique 
and  are  recorded  on  photographic  films.  The  images  recorded  contain 
70  V'  20  resolution  elements. 

As  a  result  of  the  theoretical  and  experimental  study,  it  is  found 
that  the  monochromatic  image  upconversion  process  can  be  nearly  diffraction 
limited. 
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CHAPTER  2.  GEOMETRIC  OPTICS  THEORY  OF  PARAMETRIC 
IMAGE  UPCONVERSION 
I.  INTRODUCTION 

Tlie  photomultiplier  detection  of  infrared  (IR)  radiation  by  means 

(1-5) 

of  the  parametric  upconversionprocess'  '  has  stimulated  great  interest  in 

the  next  step  of  upconversionof  IR images^^  the  visible  or  near-visible 

region.  The  theory  of  parame'ric  upconversion  is  well-established^^^ 

however,  image  upconversion  requires  considerations  beyond  those  for 

single -resolution  element  theory.  Firester  analyzed  the  problem  using  both 

Fourier  transform'  'and  geometric  optics  approaches.  Because  of  the 

simplifying  assumptions  in  his  analysis,  some  interesting  image  upcon- 

(15) 

verter  properties  are  not  apparent  in  his  results.  Andrews  has  reported 
resolution  calculations,  also  based  on  a  ray-tracing  m«  hod,  of  image  up- 
converters  utilizing  specific  optical  systems.  His  analysis  applies  for  a 
near -plane  wave  pump  and  therefore  does  not  cover  the  case  of  an  off-axis 
IR  object. 

In  this  paper,  the  discussion  is  limited  to  image  aspects  of  the 
upconverter;  the  problems  common  to  a  single -resolution  element  upconver- 
ter  (such  as  conversion  efficiency,  bandwidth,  and  material  selection)  are 
not  considered. 

A  formula  relating  the  location  of  the  upconvertei  image  to 
location  of  the  pump  and  the  IR  object  is  derived  in  Section  II  for  a  very  thin 
upconverter.  The  result  shows  that  a  point  IR  object  will  be  imaged  into  a 
point.  The  formation  of  the  upconverted  image  by  a  thick  upconverter  is 
considered  in  Section  III.  The  finite  thickness  of  the  nonlinear  material 
causes  aberrations  peculiar  to  the  parametric  image  upconversion. 


7 


II. 


IMAGE  FORMATION  BY  INFINITESIMALLY  THIN  IMAGE 
UPCONVERTER 


This  section  contains  a  discussion  of  the  formation  of  the  sum- 
frequency  image  of  an  IR  object  by  an  infinitesimally  thin  upconverter.  In¬ 
finitesimally  thin  means  that  the  thickness  of  the  nonlinear  material  is  small 
compared  with  other  dimensions  such  as  the  object,  pump,  and  image  dis¬ 
tances  and  the  transverse  dimensions  of  the  crystal,  but  that  it  is  larger 
than  the  wavelengths  involved;  that  is,  large  enough  for  traveling -wave 
parametric  interaction  inside  the  nonlinear  material.  For  this  case,  it  will 
be  shown  that  a  point  IR  object  is  imaged  to  a  point  sum -frequency  image. 

The  image  location  is  obtained  by  tracing  spherical  pump  rays, 
spherical  IR  object  rays,  and  spherical  sum-frequency  rays  that  are  para¬ 
metrically  interacting  in  the  nonlinear  material.  For  simplicity,  the  anal¬ 
ysis  is  carried  out  in  a  two-dimensional  model. 

To  further  simplify  the  analysis,  the  following  assumptions  are 

made: 

1.  The  pump,  the  IR  object,  and  the  sum-frequency  waves  are 
in  near -phase -match  condition  throughout  the  interaction 
volume;  the  slight  phase  mismatch  is  neglected. 

2.  Because  all  rays  are  paraxial,  all  the  angles  are  small  enough 
to  use  the  approximation  tan  0  »  0. 

3.  The  pump  is  a  point  source.  The  pump  beam  generated  by  a 
single -transverse -mode  laser  can  be  considered  as  a  point 
source.  (The  effect  of  an  extended  pump  source,  which  repre¬ 
sents  a  multi -transverse -mode  laser  beam,  will  be  considered 
later. ) 
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4.  The  reflected  components  of  the  pump,  the  IR  object,  and 
especially  the  sum-frequency  and  hit?her-order  scattered 

/ig\ 

waves'  generated  at  the  surface  of  (he  material  are  ignored 
because  these  rays  do  not  contribute  to  the  image  formation. 
The  first  assumption  is  justified  since  the  maximum  allowable 
phase  mismatch  is  small  compared  with  the  other  wave  vectors  involved. 
Assumptions  2  and  3  are  discussed  in  Section  IV. 

To  provide  a  consistent  convention  for  the  use  of  positive  and 
negative  quantities  in  the  analysis,  simple  Cartesian  coordinates  are  used. 
If  a'  and  j/)'  are  defined  as  the  angles  of  the  IR  wave  vector  K.  and  the  sum- 
frequency  wave  vector  K^,  measured  from  the  pump-wave  vector  (Fig¬ 
ure  1)  mside  the  nonlinear  material,  the  phase-match  condition  results  in 
the  folbwing  relationship: 

K. 

sin  i//  =  —  sin  a' 

K 

S 


For  paraxial  approximation, 

s 

a'  =  0  a'  (1) 

s 

Tlie  angular  de magnification  factor  0  =  KVK  is  not  constant  because  I  K,  I 

is  1 

and  |k  I  change  with  direction;  however,  for  a  small-angle  approximation, 
s 

the  change  is  neglected. 

Because  of  the  change  in  the  indices  of  refraction  at  the  air- 
crystal  interface,  the  direction  of  the  wave  vectors  inside  the  nonlinear 
material  will  be  modified  according  to  Snell's  law,  as  shown  in  Figure  2b. 
The  relationships  between  the  angles  are; 
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FIGURE  1.  PHASE-MATCH  CONDITION  FOR  PARAMETRIC 
UPCONVERSION 
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1 


(2) 


9 

_P^ 

n 

P 


n. 

1 


(0  +  O') 

p 


where 


0^(£  =  i,  p,  s)  =  angle  between  U  =  i,  p,  s)  and 
the  reference  axis  normal  to  the 
interface 

n^(£  =  i,  p,  s)  =  index  of  refraction  for  the  wave  designated 
by  the  subscript  I 

Tlie  primed  angles  are  measured  inside  the  nonlinear  material,  and  the  un¬ 
primed  angles  are  measured  in  the  air.  The  angles  ip'  and  a'  are  obtained 
from  equations  1  and  2: 


ip'  =  0 


(3) 


i  *\  i  p 

/ 

TTie  direction  of  the  sum -frequency  wave  exiting  from  the  infinitesimally  thin 
nonlinear  material  is  then  by  Snell's  law: 


e  =  n  0'  =  (0'  +  ip')  n  (4) 

s  s  s  p  s 

Figure  2a  shows  the  geometry  of  upconverter  image  formation.  The  point  IR 
object  is  located  at  (s,  h),  the  point  pump  source  is  located  at  (p,  o),  and  the 
sum -frequency  image  is  formed  at  (s',  h').  IR  ray  1  makes  an  angle  0^  and 
pump  ray  2  makes  an  angle  0  with  the  reference  axis;  they  interact  at  the 
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(b)  ANGULAR  RELATIONS 

FIGURE  2.  FORMATION  OF  IMAGE  WITH  INFINITESIMALLY  THIN 
UPCONVERTER 
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material  located  distance  y  from  the  reference  axis  and  generate  sum-fre¬ 
quency  ray  3.  The  direction  of  ray  3  is  calculated  from  equation  4.  The 
intersection  of  ray  3  with  ray  4,  which  is  generated  at  another  point  in  the 
material,  is  the  sum-frequency  image  and  is  located  at  (s',  h'). 

From  the  geometry  of  Figure  2a, 


h' 


=  e 


8 


y  -  h 
s 


(5) 


and 


Substituting  equations  2,  3,  and  4  into  equation  5, 


1  1 

- 1  +  - 

—  0  1 

L'* 

_s  1 

n  p 
P 


h'  h  "s 

-T - 0 

8  8  n. 


=  0 


This  equation  must  held  for  any  y  for  the  upconverter  to  form  an  image.  Each 
quantity  inside  the  brackets,  therefore,  must  be  equal  to  zero: 


1  /3  1 


s'n  sn  ~  n  p 
s  i  p 


h  «  s 

—  =  M  =  0  - 

h  t  “j® 


(6) 

(7) 


Equation  6  is  the  thin -lens  formula  for  the  iniage  upconverter,  and  equation 
7  gives  the  transverse  magnification.  Equations  6  and  7  and  the  geometry 
shown  in  Figure  2  are  for  a  case  of  a  real  point  pump  source  and  a  real  ob¬ 
ject  located  in  the  right  half -plane.  Similar  equations  can  be  derived  for 
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other  arrangements.  The  phase-match  condition,  however,  limits  the  num¬ 
ber  of  arrangements  that  can  be  used  for  an  image  upconverter;  it  restricts 
the  entrance  angles  0^^  and  0^,  and  limits  the  angular  aperture  of  an  upcon¬ 
verter  (as  discussed  in  Section  IV).  For  instance,  most  upconverters  using 
presently  available  nonlinear  crystals  do  not  allow  an  arrangement  in  wliich 

a  real  pump  source  and  a  real  object  are  located  on  opposite  sides  of  the 
(15) 

crystal.  Equations  6  and  7  apply  to  the  following  cases,  which  permit 
upconversion  to  take  place. 

1.  A  real  point  pumu  source  and  a  real  IR  object  on  the 
same  side  of  the  material  (Figure  2).  The  image  is 
virtual,  erect,  and  on  the  same  side  as  the  object. 

2.  A  real  point  pump  source  and  a  virtual  object  arc  on 
opposite  sides  of  the  material.  The  image  can  be 
either  virtual  or  real. 

3.  A  virtual  point  pump  source  and  a  real  object  are  on 
opposite  sides  of  the  material.  The  image  can  be 
either  virtual  or  real.  This  is  a  dual  of  case  2. 

4.  A  virtual  point  pump  source  and  a  virtual  object  are 
on  the  same  side  of  the  material.  The  in;age  is  on 
the  same  side  of  the  virtual  object  and  real.  This 
case  is  a  dual  of  case  1. 

A  virtual  point  source  represents  a  converging  spherical  wave, 
and  a  real  point  source  represents  a  diverging  spherical  wave. 

It  is  sometimes  convenient  to  describe  an  optical  system  by  a 
paraxial-ray  transfer  matrix.  It  relates  the  position  y,  and  the  slope  (in¬ 
clination)  0  of  the  ray  in  the  input  plane  of  the  system  to  the  position  y^  and 
the  slope  0^  of  the  ray  in  the  exit  plane  of  the  system  as  follows; 
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If  (Vj,  Og)  and  {y^,  0^)  are  defined  as  the  IR  ray  enl ranee  position 
and  slope,  and  the  sum -frequency  ray  exit  position  and  slope,  respectively, 
an  image  upconverter  by  a  ray  transfer  matrix  tan  be  described.  From 
equations  3,  4,  and  5,  for  an  infinitesimally  thin  upconverter; 


llie  determinant  of  the  ray  matrix  is  clearly  not  unity,  as  ex])ected  from  the 

nonreciprocity  of  the  upconversion  process.  The  propagation  of  paraxial 

rays  through  an  image  upconverter  system  containing  optical  components 

such  as  lenses  and  transmission  media  can  now  be  obtained  by  multiplying 

the  ray-transfer  matrix  of  each  optical  element. 

in.  ABERRATIONS  OF  IMAGE  UPCONVERTER  WITH  NONLINEAR 
MATERIAL  OF  FTM TE, THICKNESS 

1.  IMAGE  FORMATION  WITH  THICK  MATERIAL 

Hie  discussion  of  image  formation  by  an  infinitesimally  thin 

image  upconverter  in  Section  n  has  furnished  a  relatively  simple  method 

of  calculating  the  position  and  size  of  the  image.  In  practice,  the  thickness 

of  nonlinear  material  is  neither  infinitesimal  nor  infinite  so  that  the  effect  of 

material  thickness  on  image  formation  must  be  considered.  The  straight 

forward  method  of  attacking  this  problem  is  to  trace  all  rays  interacting  at 

different  transverse  planes  inside  the  material. 
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All  the  assumptions  made  in  Section  11  are  considered  to  be  valid 
throughout  the  discussion  in  this  section. 

Figure  3  illustrates  the  geometry  of  ray  tracing.  The  object, 
image,  and  pump  rays,  which  interact  at  transverse  plane  T,  are  traced 
using  a  procedure  similar  to  that  described  in  Section  11.  Transverse 
plane  T  is  located  at  distance  t  from  the  center  of  material  C,  which  is  taken 
as  the  origin  of  the  coordinate  system.  The  pump,  object,  and  image  are 
located,  respectively,  at  (p,  o),  (s,  h),  and  (s',  h').  The  points  O'  and  P' 
are,  respectively,  the  object  and  pump  locations,  corrected  for  refraction 
at  the  material  surface.  The  points  O',  N,  and  L  are  on  a  straight  line 
representing  an  object  ray  inside  the  material  and  making  angle  9^'  with 
respoc .  lo  the  reference  axis  X.  The  straight  line  P'ML  is  an  equivalent 
pump  ray.  The  straight  line  LNO'  is  an  image  ray,  making  an  angle  9^'  with 
the  X-axis,  inside  the  material.  The  straight  line  KI,  making  an  angle  9^ 
with  the  X-axis,  is  an  image  ray  in  air. 

The  following  relationships  are  obtained  by  geometry: 

From  aO'AL, 


From  A  P'BL, 


From  AOGN, 


From  A  PMF, 


9 

P 


AL 

BL  -  h 

O'A  ■ 

O'A 

9  '  = 
P 

BL 

P'B 

GN 

NF  -  h 

OG  D 

8  -  -r 


FM 

FM 

PF 

D 

P-2 

(10) 

(11) 

(12) 

(13) 
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FIGURE  3.  IMAGE  FORMATION  FOR  PARAMETRIC  UPCONVERTER 
WITH  FINITE  THICKNESS  NONLINEAR  MATERIAL 
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From  AO'GN, 


From  A  P'MF. 


Using  Snell's  law  and  equations  10,  12,  and  14, 


0.' 

1 


BL  -  h 


Similarly,  from  equations  11,  13,  15,  and  Snell's  law. 


By  geometry. 


and 


Q  /  3  _ 

y^  =  K.I=LB.e;(?.t) 


0 

S 


Substituting  the  relationship 

9  '  =  9  '  (1  -  0)  +  06  ' 

S  p  1 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


obtained  from  the  phase -matching  condition  and  after  rearranging  equations 
18  and  19, 


18 


(s '  +  D/2) 


-1 


8) 


1  + 


(D/2  +t)(l-  8) 


(D/2  +  t)  8 


n  (P  -  D/2)  +  (D/2  -  t)  n.(s  -  D/2)  +  (D/2  -  t) 


n  8 
s 


n  (p  -  D/2) 
h' 


+  (D/2  -  t) 


n.(s  -  D/2)  +(D/2  -  t) 


8h 


/  D/2  +  t 


(s  -  D/2)  n.  +  (D/2  -  t) 


n  - 
6 


s  + 


D/2 


/  -• 


=  0(21) 


s '  +  D/2 
where  y  =  BL 

Since  equation  21  must  hold  for  every  y  and  h  for  image  forma¬ 
tion,  each  member  within  the  braces  {  }  must  be  separately  equal  to  zero. 
Equations  expressing  the  image  generated  at  a  transverse  plane  T  are  then 
given  by: 

_ 1 _  8 _ L-1 


4"s  )  '(^•^01 


(22) 


h'  =  h 


(23) 


The  image  location  (s',  h')  is,  as  expected,  a  function  of  both 
D  and  t.  The  images  generated  at  the  different  transverse  planes  of  a  single 
object  are  formed  at  different  locations.  The  dependence  of  the  image  posi¬ 
tion  s'  upon  t  leads  to  a  blurring  of  the  image.  This  effect  is  analogous  to 
lens  spherical  aberration  but  is  peculiar  to  the  image  upconverter.  The 
thickness  aberration  is  due  to  the  finite  thickness  of  the  material,  whereas, 
lens  spherical  aberration  is  due  to  the  obliquity  of  nonparaxial  rays. 
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Spherical  abei  ration  would  be  added  to  the  thickness  aberration  if  exact-ray 
tracing  is  used  rather  than  paraxial-ray  tracing.  Other  monochromatic 
aberrations  (such  as  coma,  astigmatism,  curvature  of  field,  and  distortion) 
would  have  to  be  considered  if  higher-order  theory  were  used  for  the  analysis. 

Equations  22  and  23  have  been  derived  for  a  case  where  both  the 
pump  source  and  the  object  are  located  outside  the  nonlinear  material.  If  a 
pump  source,  either  real  or  virtual,  is  located  inside  the  material,  equa- 


A  similar  modification,  m  =  1,  should  be  made  for  the  case  where  an  IR  ob¬ 
ject  is  inside  the  material.  The  modification  is  necessary  for  this  case  be¬ 
cause  the  refraction  of  the  pump  and  object  rays  at  the  front  surface  of  the 
crystal  does  not  affect  the  imaging  process.  The  refraction  of  the  sum- 
frequency  rays,  however,  must  be  taken  into  account  regardless  of  the  pump 
and  the  object  locations,  because  the  image  will  be  observed  in  air  rather 
than  inside  the  nonlinear  material.  The  existence  of  chromatic  aberration 
is  evident  from  the  dispersion  of  the  material.  Chromatic  aberration  exists 
even  in  an  infinitesimally  thin  image  upconverter.  Since  chromatic  aberra¬ 
tion  is  not  perculiar  to  the  image  upconverter,  even  though  it  may  differ 
somewhat  from  lens  chromatic  aberration,  it  will  not  be  considered  further. 

2.  PLANE -WAVE  IR  BEAM 

Plane -wave  IR  beam  incidence  parallel  to  the  reference  axis  can 
be  presented  by  equating  the  object  distance  s  to  infinity.  The  sum-frequency 
image  position  is  then  given  by: 


"s  ‘‘  ■ 


p 

2 


(24) 
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The  sum-frequency  wave  focal  points,  generated  at  (he  front  and  back  sur¬ 
faces  of  the  material,  are  obtained  by  equating  t  to  r  and  -  respectively: 


(25) 


3.=  "-^) 

b  n  (1  -  0) 


for  t  =  -  I 


(26) 


The  longitudinal -thickness  aberration,  defined  as  the  difference  between 
these  two  focal  points,  is  then: 


=  I V  -  *f' I  =  I  irfc)  I 

s 


(27) 


llie  longitudinal-thickness  aberration  for  plane-wave  IR  beams  depends  only 
on  the  material  properties,  but  not  upon  the  pump  position.  The  lateral¬ 
thickness  aberration  is  given  by: 


Ay' =  As' tan  9  (28) 

s  max 

where  0  is  the  maximum  angle  that  the  sum-frequency  ray  makes  with 
s  max 

the  reference  axis. 

3.  PLANE -WAVE  PUMP  BEAM 

Plane -wave  pump-beam  incidence,  normal  to  the  material  is  con¬ 
sidered  next.  For  this  case,  the  pump  distance  p  becomes  infinite  and  the 
image  is  located  at: 
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'Hu*  l(in(;itu(liii;il -il)i(*ki)C‘S.s  iilK'rratioii  becomes 


A  s'  1  s, '  -  s  '  I  Id  --i'  1  (30) 

1)  f  n 

s 

It  Ls  imtepeiulciit  of  the  III  object  location.  The  lon(;itudinal-tinckness  aber¬ 
rations  of  the  plane -wave  III  l)eam  and  the  plane -wave  pump  l)eam  are  re¬ 
lated  by 


(^s') 


a 


Usl  ,,  ..2 

p  (1  -  (i) 


(31) 


Since  is  usually  less  than  unity  for  imajte  upconvcrsion,  the  Ion(;itudinaI- 
Ihickness  al)erration  is  smaller  for  the  parallel -object -ray  case  than  for  the 
jurallel -pump-beam  case. 

llie  lateral -thickner^^s  aberration  is  given  also  by  equation  28. 


4.  PUMP  SOURCE  AND  IR  OBJECT  IN  SAMLJrJRANSVERSE  PLANE 

The  image  |)osition  is  given  by  equations  22  and  23,  using  the 

relationship  s  =  p.  Since  both  s'  and  h'are  functions  of  t,  the  blurring  of  the 

(14) 

image  could  not  lx?  avoided.  Because  Firester  ignored  refraction,  he 
concluded  that  resolution  for  this  optical  arrangement  was  unlimited.  The 
non- zero  aberration  is  due  to  dispersion  of  the  material.  Indeed,  s'  will 
be  independent  of  t  i:  the  material  is  nondispersive.  However,  h'  is  s'ill 
dependent  on  the  object  iK)sition.  If  p  =  s  and  '^g  "  •'p  ■  equation  22 

and  23.  then 

s'  =  p+d(^-i)  (32) 
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Since  s'  is  independent  of  t,  the  longitudinal -thickness  aberration  does  not 
exist;  also  the  transverse  image  spread  defined  by  ^h'  I  h^/  -  b^^'  1  becomes 

Ah'=0  (34) 

Above  result  indicates  that  the  obliquity  of  the  IR  object  ray  would  not  cause 
image  blurring  if  the  crystal  disix;rsion  can  be  neglected.  Image  blurring 
Ah'  of  a  dispersive  crystal  due  to  the  object  ray  obliquity  is,  however,  not 
zero  in  general.  This  effect  resembles  the  coma  of  an  ordinary  lens.  This 
effect  is  called  thickness  coma  since  it  is  due  to  the  finite  thickness  of  the 
material,  but  it  is  not  caused  by  violation  of  par;L\ial  approximation  as  in 
t.he  case  of  an  ordinary  lens. 

The  longitudinal-thickness  aberration.  As',  of  a  dispersive 
material  for  the  case  s  -  p,  is  generally  not  zero  but  is  minimum. 

5.  PARALLEL-PUMP  AND  PARALLEL-OBJECT  RAYS 
Pump  and  object  rays  can  be  made  parallel  to  the  axis  with  a 

suitable  lens  arrangement.  Since  both  s  and  p  are  equal  to  infinity,  the  image 
will  be  formed  at  infinity.  Thickness  aberration  and  thickness  coma  should 
disappear  for  this  case.  This  arrangement  represents  upconversion  in 
Fourier  space  and  is  an  aberration-free  image  system.  The  resolution  for 
this  arrangement  is  limited  by  higher-order  effects  discussed  later. 

6.  PUMP  AND  OBJECT  AT  CENTER  OF  MATERIAL 

Pump  and  object  can  be  positioned  at  the  center  of  the  material 
with  a  suitable  optical  system.  The  pump  and  object  distances  are  zero  for 


this  case.  The  refractive  indices  and  n^  are  equal  to  unity  because  the 
pump  and  the  object  are  inside  the  material. 

The  image  locations  s'  and  h'  are  given  by 


s'  =  ^ 
®  2 


n 


h' 


(35) 

(36) 


Both  the  iongitudinal  thickness  aberration  and  the  thickness 
coma  are  zero  as  is  evident  from  Equations  (35)  and  (36). 

7.  COMPUTED  VALUES  OF  ABERRATION  AND  COMA  FOR 

SPECinC  CASES 

Consider,  next,  a  specific  upconverter  configuration  using  a 
proustite  (Ag^AsS^)  crystal,  pumped  by  an  Nd  :  YAG  laser  (1.06  pm),  and  a 
signal  from  a  COg  laser  (10.6  pm).  For  this  typical  values  of  thickness 
aberration  and  c^ma  have  been  computed. 

Figure  4  shows  the  longitudinal  aberration  that  occurs  due  to 
the  longitudinal  thickness  of  the  nonlinear  crystal  material.  Longitudinal 
aberration  takes  the  image  of  a  point  and  elongates  it  along  the  optical  axis 
of  the  upconverter  system.  The  values  of  pump  position  (p),  object  position 
(s),  and  aberration  (As')  have  been  normalized  with  respect  to  the  thickness 
(D),  of  the  crystal  material.  Each  curve  corresponds  to  a  different  value 
of  normalized  pump  position  (p/D)  with  respect  to  the  region  of  nonlinear 
interaction.  When  the  apparent  object  distance  (s/D)  is  equal  to  the  apparent 
pump  distance,  that  is,  the  spatial  wavefronts  are  matched,  the  longitudinal¬ 
thickness  aberration  becomes  a  minimum.  Other  requirements  permitting 
upconverter  operation  should  be  set  at  this  value. 
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As  the  pump  distance  goes  to  infinity  (p/D  =  «),  that  is,  the  pump 
wave  becomes  a  plane  wave,  the  normalized  aberration  achieves  a  fixed 
value  of  3.9.  If  both  the  object  distance  and  the  pump  distance  are  infinite, 
the  case  covered  in  Section  III-5,  the  aberration  becomes  zero. 

Figure  5  shows  the  thickness  coma  values  computed  for  a  fixed 
pump  position.  This  curve  is  representative  of  a  family  of  such  curves  and 
indicates  ?  decrease  in  coma  as  the  curvature  of  the  signal  wavefront  is 
decreased.  For  small  values  of  off-axis  displacement  (h),  the  values  of 
coma  will  generally  be  very  small  for  large  values  (  I  s/D  I  >  10)  of  object 
distance. 

IV.  EXTENDED  PUMP  SOURCE.  ANGULAR  APERATURE.  AND 
RESOLUTION 

Degradation  of  image  quality  due  to  an  extended  pump  source,  as 
well  as  the  angular  aperture  of  a  parametric  image  upconverter,  are  con¬ 
sidered  in  this  section.  Formulae  for  resolution  calculation  are  also  pre¬ 
sented. 

1 .  EXTENDED  PUMP  SOURCE 

In  general,  a  pump  source  used  for  an  image  upconverter  is  a 
laser.  A  pump  laser  can  be  operated  either  in  a  single -transverse  mode 
or  in  a  multitrans verse  mode.  The  pump  wavefront  of  a  single -mode  laser 
beam  will  always  be  a  single -spherical  wave  or  ne^r-plane  wave.  The  as¬ 
sumption  of  a  point -source  pump  with  a  spherical  wavefront  used  in  the 
preceding  sections  therefore  applies  to  this  case.  A  multitransverse  mode 
laser  beam  can  be  considered  as  a  composite  of  many  spherical  wavefronts. 
The  pump  source  can  then  be  considered  as  an  extended  source  of  finite 
size.  Hie  pump  radiation  at  each  point  in  the  material  will  be  represented 
by  wave  vectors  with  a  range  of  directions  equal  to  the  divergence  of  the 
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FIGURE  4.  NORMAUZED  THICKNESS  ABERRATION  AS  A  FUNCTION  OF 
AXIAL  OBJECT  DISTANCE,  WITH  PUMP  POSITION  AS  A 


10 


-60  -50  -40  -30  -20  -10  0  10  20  30  40  50  60 


1-3277  S/D 

FIGURE  5.  THICKNESS  COMA  OF  AN  IMAGE  UPCONVERTER  AS  A 
FUNCTION  OF  AXIAL  OBJECT  DISTANCE 
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pump  beam.  The  size  of  an  extended  pump  source  is  related  to  the  diver¬ 
gence  angle  6,  for  a  small-angle  approximation,  by 

a  =  p6  (37) 

where 

a  =  equivalent  pump  size 

5  =  pump  divergence  angle 

The  divergence  of  the  pump  waves  will  cause  a  divergence  in  the 
sum -frequency  waves  for  a  fixed  object-wave  vector.  The  same  sum -fre¬ 
quency  divergence  could  be  caused  by  an  equivalent  object  divergence  if  the 
pump  divergence  6  is  made  zero.  The  equivalent  object  divergence  €  that 
would  give  the  same  sum -frequency  divergence  as  the  pump  divergence  can 
be  obtained  from  the  angular  relationship 


9  =  (1  -  3)  e  +  ^  0 

s  n  p  n.  1 

P  1 


(38) 


The  change  in  0^  due  to  the  change  in  0^,  keeping  0^  constant  is 


30  n 

^  A0.  =  0  P 

1 


30 


(30) 


Similarly 


90 

s 


90. 

1 


(40) 
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The  equivalent  object  divergence  €  is  then  obtained  by  equating 
equations  39  and  40; 


€ 


n. 


1  -  /3 


n 


(41) 


An  image  upconverter  with  a  point  object  and  an  extended -pump 
source  can  now  be  represented  by  an  equivalent  system  consisting  of  a  point 
pump  source  and  a  finite  object.  The  angular  resolution  of  an  upconverter 
with  a  multitransverse-mode  laser  pump  is  therefore  limited  by  the  angle  e. 

2.  ANGULAR  APERTURE 

To  be  useful,  an  image  system  should  have  a  reasonably  wide 
field  of  view.  The  principal  factor  that  restricts  the  field  of  view  of  a  param¬ 
etric  image  upconverter  is  the  phase -matching  condition. 

For  crystal  of  length  D,  phase -matched  upconversion  takes  place 

277 

as  long  as  I  AK  l=IK  -K  ,  where  AK  is  assumed  to  be 

-  — s  -p  — i  D  — 

along  the  direction  of  K  for  convenience,  as  shown  in  Figure  6.  Three - 

— s 

dimensional  phase  matching  instead  of  two-dimensional  matching  should  be 
considered  for  an  image  converter  if  the  anisotrophy  of  the  crystal  is  signi¬ 
ficant.  Only  the  two-dimensional  match,  shown  in  Figure  4,  is  discussed. 

The  phase -match  error  AK  is  given  by 


AK  =  K  -  (K  +  K.  cos  a')  sec  ip' 
s  p  1 


(42) 


where 


K.  sin  a' 

tan  lb*  = - - - 7 

^  K  +  K.  cos  a' 
P  1 


(43) 


The  small-angle  approximation  sin  x~x  reduces  equation  43  to 
the  angular -demagnification  relation.  The  previous  sections  assumed  that 
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FIGURE  6.  ANGULAR  RELATIONS  IN  CRYSTALLOGRAPHIC 
COORDINATE  SYSTEM 
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AK  is  negligibly  small  compared  with  the  smallest  K.  Given  the  values  of 

0',  a',  and  0  (the  angle  between  the  crystallographic  C  axis  and  the  normal 
P 

to  the  crystal  surface),  A  K  can  be  obtained  using  equations  42  and  43. 

Figure  7  shows  the  numerical  results,  using  refractive  index  data  reported 
(17) 

by  Hobden  ,  for  a  1-cm  long  Ag  AsS.  upconverter.  The  IR  object  and  the 

O  O 

pump  wavelengths  are  10.63  pm  and  1.064  pm,  respectively;  they  are  polar¬ 
ized,  respectively,  extraordinary  and  ordinary.  The  angles  and  (p!  are 
related  to  0j^  and  0' ,  respectively,  by: 

6^  =  -  0'  (jt  =  i  or  p)  (44) 

It  has  been  shown  that,  in  general,  noncritical  phase  matching  gives  a  larger 

(18) 

field  of  view  for  the  object  field  than  other  phase -matching  arrangements. 

This  is  true  only  if  nonmonotonical  angular  variations  of  sum-frequency  out¬ 
put  are  not  allowed.  This  situation  is  demonstrated  in  Figure  7  by  the  lines 
marked  a,  b,  and  d.  The  change  in  A K  is  monotonic  when  the  IR  object  inci¬ 
dent  to  angle  0.  is  monotonic  along  lines  a,  b,  and  d.  Consequently,  the 
angular  sum -frequency  output  variation  exhibits  a  single  maximum  where 
I  A  K  I  is  minimum.  It  is  clear  from  Figure  7  that  line  b  has  the  largest 
object  angular  aperture.  The  field  of  view  can  be  increased  slightly  if 
double  maxima  in  thesum -frequency  output  angular  variation  is  allowed. 

Line  c  in  Figure  7  represents  this  case. 

For  a  1-cm  long  Ag  AsS  image  upconverter  with  a  10.6-pm 

U  O 

object  and  a  1.  06- pm  pump,  the  maximum  angular  aperture  lies  within  the 
shaded  area  of  Figure  7.  The  allowable  maximum  I  A  K  1  is  chosen  to  be 
27t  for  convenience;  the  sum-frequency  output  is  zero  at  this  value  of  AK. 
Figure  8  shows  the  field  of  view  of  a  1-cm  long  Ag„AsS  image  upconverter. 

O  U 
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FIGURE  7.  PHASE-MATCH  ANGLES  FOR  Ag,AsS,  10.  6  TO  0. 96  urn 
UPCONVERSION  ^  ^ 
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FIGURE  8.  CRYSTAL  ORIENTATION  AND  ANGULAR  APERTURE  FOR 
AggAsSg  10.6  TO  0.96  m  IMAGE  UPCONVERTER 
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The  crystal  should  be  cut  so  that  tlie  optical  axis  makes  an  angle  of  20.  3  de¬ 
grees  with  the  normal  to  the  surface.  The  angular  aperture  of  other  non¬ 
linear  materials  suitable  for  image  upconversion  has  been  shown  to  be  less 

(15) 

than  a  few  hundred  niilliradians.  Therefore,  the  paraxial  assumption 
used  in  the  preceding  discussion  of  angular  aperture  is  shown  to  be  valid. 

3.  RESOLUTION 

The  resolution  of  an  image  upconverter  can  be  calculated  if  the 
aberrations  and  the  angular  aperture  of  the  system  are  known.  Only  the 
point  pump  source  is  considered  here.  The  modification  to  an  extended  pump 
source  is  straightforward  since  its  affect  on  the  object  field  is  known. 

Figure  9  shows  the  image  of  a  point  IR  object  located  at  (s,  h); 
the  pump  is  located  at  (p,  o).  The  image  is  a  line  connecting  points  and 
which  are,  respectively,  the  images  generated  at  the  front  and  the  back 
surfaces  of  the  material.  Images  generated  on  a  transverse  plane  located 
between  these  extreme  surfaces  lie  along  the  Because  of  thickness 

aberration,  the  lateral  aberration  is  not  zero;  also,  the  thickness  coma 
produces  the  transverse  spread  Ah'.  The  transverse  dimension  of  the  image 
is  the  sum  of  these  two  aberrations: 


H'  =  Ah'  +  As'  0 


s  max 


h;-h;|.|s;-sje 


s  max 


(45) 


Again,  0  is  the  maximum  angle  that  the  sum -frequency  ray  makes  with 
s  tnH.x 

the  reference  axis. 

To  analyze  the  resolution  of  an  image  upconverter,  consider  the 

geometry  of  Figure  10.  Two  object  points,  located  at  (s,  h  )  and  (s,  h  ),  are 

1 

imaged  into  lines  and  I^.  In  order  that  these  two  images  be  resolvable, 

the  point  (s'  ,  h'  )  and  the  point  (s'  h/ _  -  As'  0  )  should  coincide, 

fl  fl  f2  d2  s  max 
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NOTE:  REFRACTION  AT  THE 
MATERIAL  BOUNDARY  IS  NEGLECTED 
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FIGURE  9.  BLURRING  OF  SUM- FREQUENCY  IMAGE  DUE  TO 
ABBERRATIONS  OF  PARAMETRIC  UPCONVERTER 
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f 


148) 

h'  =h'  -As'e  147) 

fl  b2  s  mux 

Equation  46  is  automatically  satisfied.  The  chan(!;e  in  9 

s  max 

due  to  the  difference  in  object  transverse  location  is  neglected.  Usinij;  equa¬ 
tion  23,  equation  47  is  rewritten  as 

h,  M.,  =  h-M,.  -  ^s'e  (48) 

1  tf  2  tb  s  max 

The  minimum  resolvable  transverse  size  of  an  object  located 


at  (s,  h)  becomes 


Ah  =  hl~-1 
tb 


As' 

+  ^  0 
M..  s  max 
tb 


Hie  linear  resolution  of  the  object  field  about  (s,  h)  is  then  given  by 

R  =  “  elements/unit  length  (50) 

The  angular  resolution,  y,  which  is  more  meaningful  than  the  linear  resolu¬ 
tion  for  some  optical  systems,  for  a  small-angle  approximation  is  given  by 


y  =  Rs/Ah  (51) 

The  resolution  of  a  plane-wave  pump  case  is  now  considered  as 

"s 

an  example.  For  this  case,  p  =  <»,  0  =0,  and  0  =  —  /30, 

p  s  max  ni  i  max 

according  to  equation  20. 

The  minimum  resolvable  transverse  size  is  equal  to  the  lateral - 
thickness  aberration,  because  Ah'  is  zero  in  this  case: 

(Ah)  -D^-^0.  (52) 

p='»  n.  1  max 
1 
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The  linear  resolution,  given  by 


R 
P  = 


00 


n. 

1 


D(1  -  /3)G. 


1  max 


(53) 


is  independent  of  the  object  location. 

V.  OTHER  FACTORS  AFFECTING  IMAGE  FORMATION  OF  IMAGE 
UPCONVERTER 

A  nonlinear  crystal  used  for  upconversion  must  be  anisotropic 
for  parametric  interaction  to  take  place.  The  electro-optic  tensor,  which 
is  responsible  for  the  generation  of  nonlinear  polarization  in  the  material, 
is  nonzero  if  the  crystal  does  not  possess  a  center  of  symmetry.  Anisotropy 
and  dispersion  are  also  required  for  phase  match.  The  treatment  presented 
in  the  preceding  sections  does  not  include  the  effect  of  anisotropy.  The  ef¬ 
fect  of  dispersion  was  partially  accounted  for  by  using  different  indices  of 
refraction  for  different  wavelengths. 

The  first  effect  of  anisotropy  is  due  to  the  angular  change  of 
the  index  of  refraction  of  the  extraordinary  wave.  This  change  also  affects 
the  angular  demagnification  factor,  since  it  is  a  ratio  of  the  IR  wave  vector 
to  the  sum -frequency  wave  vector.  The  image  of  a  point  object  generated 
even  by  an  infinitesimally  thin  material  will  not  be  a  point  if  the  angular 
change  of  the  index  of  refraction  is  included. 

The  second  effect  of  anisotropy  is  due  to  double  refraction.  The 
wave  vector  normal  to  the  wavefront  and  the  Poynting  vector  are  not  in  the 
same  direction  in  an  anisotropic  crystal.  The  analysis  given  in  the  preceding 
sections  is  not  rigorous  ray  tracing  because  of  double  refraction.  The  ray 
tracing  inside  the  material  should  be  done  along  the  rays  instead  of  along 
the  wave  vectors.  Snell’s  law  should  be  used  for  rays  instead  of  for  the 
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wavefront,  and  the  angular  demagnification  factoi-  should  bo  modi  lied  lor  ray 
angles.  Also,  the  effective  interaction  thickness  of  the  upconvoi'ter  will  be 
different  from  the  material  physical  thickness  if  double  refraction  is  con¬ 
sidered. 

Another  effect  of  dispersion  is  variation  of  the  angular  demagni¬ 
fication  factor  with  wavelength.  The  natural  consequence  of  dispersion  is 
chromatic  aberration,  which  is  important  if  the  object  spectrum  or  the  pump 
spectrum  is  wide.  Image  degradation  caused  by  chromatic  aberration  will 
not  be  severe  for  most  of  the  presently  available  optical  upconverters  because 
of  narrow  infrared  spectral  bandwidths,  ty]5ically  in  order  of  10  ^  cm.  It, 

however,  may  become  predominant  for  a  large  bandwidth  upconverter,  such 

(19) 

as  the  one  reported  by  Midwinter. 

Another  factor  affecting  the  image  formation  is  the  three-dimen¬ 
sional  effect.  The  ray-tracing  analysis  presented  in  previous  sections  is 
carried  out  in  two  dimensions.  The  reference  axis  normal  to  the  crystal 
surface  from  the  point  pump  is  not  a  symmetry  axis.  This  means  that  a 
parametric  image  upconverter  is  not  circular  symmetric.  The  angular 
demagnification  factor  and  the  refractive  indices  depend  on  the  transverse 
coordinates  (x,  y)  (not  the  same  as  the  crystallographic  coordinates),  where 
the  interaction  takes  place.  This  is  a  three-dimensional  generalization  of 
the  double-refraction  effect  previously  discussed.  If  the  nonlinear  crystal 
was  oriented  so  that  the  crystallographic  C  axis  (optic  axis)  is  normal  to  the 
entrance  and  exit  surface,  and  the  angular  aperture  was  spread  almost 
evenly  about  the  C  axis,  the  image  upconverter  would  be  circular  symmetric. 
However,  the  requirement  for  phase  matching  and  maximum  light  transmis¬ 
sion  necessitate  crystal  cuts  that  do  not  permit  circular  symmetry.  A 
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three-dimensional  field  of  view  is  therefore  not  a  circular  cone.  The  Lack 
of  symmetry  introduces  distortion  of  the  upconverted  image.  A  practical 
image  upconverter,  therefore,  must  not  use  highly  anisotropic  material  for 
high-quality  imaging. 

The  discussion  of  angular  aperture  presented  in  Section  IV 
applies  to  a  spherical  pump  beam;  that  is,  the  pump  intensity  remains  con¬ 
stant  over  the  wavefront.  In  practice,  however,  the  pump  beam  is  near 
Gaussian;  its  intensity  decreases  with  the  distance  from  the  beam  axis. 

The  brightness  of  the  image  of  a  finite  object  decreases  accordingly  with  the 
transverse  distance  from  the  reference  axis.  This  means  that  the  image 
does  not  have  clear  definition  at  its  edge.  The  fading  of  the  image  also 
implies  the  reduction  in  the  field  of  view  of  an  imaging  system.  The  fading 
of  the  image  occurs  even  along  the  reference  axis,  where  it  is  well  within 
the  field  of  view  of  the  object  field.  The  resolution  calculation  becomes 
ambiguous  with  a  Gaussian  pump  beam. 

VI.  CONCLUSION 

Formation  of  a  sum-frequency  image  by  means  of  parametric 
upconversion  has  been  analyzed  for  general  pump  and  object  configurations 
by  applying  the  ray-tracing  metiiod.  A  formula  relating  the  image  location 
to  the  object  and  the  pump  locations  has  been  derived.  A  geometric  optics 
theory  of  aberrations  for  parametric  image  upconversion  has  been  described. 
The  aberrations  are  due  to  the  finite  thickness  of  the  nonlinear  material. 

This  phenomenon  is  peculiar  to  the  parametric  image -upconversion  process. 
Aberrations  due  to  the  finite  thickness  of  the  material  vanish  when  both  the 
pump  and  the  object  waves  are  parallel -plane  waves.  The  longitudinal 
aberration  is  minimum  when  the  pump  and  the  object  are  located  on  the 


same  transverse  plane.  The  limiting  resolution  of  an  image  upconverter 
can  be  calculated  when  the  aberrations  for  the  particular  image  upconverter 
are  known.  A  method  for  calculating  the  resolution  has  been  described. 

The  primary  factor  that  limits  the  field  of  view  of  a  parametric 
image  upconverter  is  the  phase-match  condition.  Other  factors,  such  as  dis¬ 
persion,  anisotropy,  and  nonparaxial  rays  degrade  the  image,  but  are 
second -order  effects. 
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CHAPTER  3  PHYSICAL  OPTICS  THEORY  OF  PARAMETRIC  IMAGE 
UPCONVERSION 

I.  INTRODUCTION 

Imaging  theory  of  parametric  image  upconversion  was  treated 
by  the  geometric  optics  approach  in  preceding  chapter.  An  equation  relating 
locations  of  object,  pump,  and  image  was  derived  and  the  existence  of  image 
aberrations  due  to  the  finiteness  of  the  nonlinear  material  thickness  was  made 
clear.  In  geometric  optics  treatment  of  parametric  image  upconversion, 
some  higher  order  factors  that  would  affect  imaging  properties  were  ne¬ 
glected.  It  was  considered  also  in  the  geometric  optics  theory  that  upcon- 
verted  image  is  formed  by  constructive  interference  between  the  waves 
radiated  fi’om  each  differential  slab  of  polarization  in  the  material.  This 
concept  lead  to  thickness  aberration  of  the  image  upconversion.  To  verify 
the  concept  of  constructive  interference  and  to  understand  the  higher  order 
effects;  such  as,  anisotropy,  phase -mismatch,  and  finite  interaction  aper¬ 
ture,  the  problem  must  be  treated  in  the  formalism  of  physical  optics.  It 
is  the  purpose  of  this  chapter  to  treat  the  problem  in  the  framework  of  wave 
theory.  The  problem  considered  here  is  to  determine  the  nature  of  the  image 
field  arising  from  a  dielectric  polarization  generated  at  the  sum-frequency 
by  interaction  of  the  object  field  and  the  pump  field  in  the  nonlinear  material. 

The  imaging  properties  of  the  parametric  image  upconversion  has  been 

(13)  (22) 

treated  by  Firester  in  the  formalism  of  Fourier  optics  that  resembles 

closely  to  the  treatment  of  linearly- invariant  filter  system.  The  advantage 

of  such  an  approach  is  that  the  process  of  upconversion  of  the  complex  non 

planar  objects  field  to  the  complex  non  planar  Image  field  can  be  studied  by 

investigating  the  interaction  of  planar  angular  spectrum  components  in  the 

nonlinear  material. 
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Tlie  results  obtained  by  Fourier  optics  treatment  describe  more  quantita¬ 
tively  the  imaging  properties  of  upconversion  process  and  can  be  tied  di¬ 
rectly  to  the  formalism  of  optical  transfer  function  description  of  an  imaging 
system. 

General  formulation  of  image  upcor.  version  process  taking  into 

account  of  phase -mismatch  and  the  crystal  anisotropy  is  presented  in  Sec- 

(23) 

tion  II.  The  treatment  presented  follows  Kleiman  closely  and  is  more 

(13) 

rigorous  than  Firester  .  The  basic  equations  that  describe  the  angular 
spectrum  of  image  field  in  terms  of  object  and  pump  angular  spectrum  are 
obtained.  The  concept  of  constructive  interference  of  radiation  from  dif¬ 
ferential  slabs  is  also  verified  in  Section  II.  Imaging  properties,  such  as 
angular  spectrum  of  image  field,  image  position,  and  the  point  spread  func¬ 
tion  of  several  image  upconversion  systems,  namely  i)  planar  wave  pump 
system,  ii)  spherical  wave  pump  system,  and  iii)  Gaussian  wave  pump  sys¬ 
tem  are  considered  in  Sec.  III.  The  results  obtained  agree  with  the  geometric 
optics  results.  The  effect  of  finite  interaction  aperture  on  upconverted 
images  is  also  discussed  in  Sec.  III.  The  effective  interaction  aperture  at 
the  sum-frequency  is  smaller  than  the  aperture  at  pump  and/or  object  fre¬ 
quencies. 

II.  GENERAL  FORMULATION 

1.  Basic  Equations 

Mathematically,  the  formulation  is  based  upon  Maxwell's  equa¬ 
tion  for  a  lossless  anisotropic  medium  free  of  charges  and  currents  but  con¬ 
taining  a  prescribed  source  in  the  dielectric  polarization.  All  linear  induced 

polarization  is  assumed  to  be  included  in  the  dielectric  constant  tensor.  The 

(23  24) 

analysis  follows  very  closely  to  Kleiman' s  treatment  ’  of  second 
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harmonic  generation  except  that  the  analysis  will  be  carried  out  in  Fourie; 
spatial  frequency  space  rather  thiui  in  spatial  coordinate  space.  It  is  shown 
in  Appendix  A  that  the  Fourier  decomposed  angular  spectrum  of  the  field 
disturbance  can  be  considered  as  plane  waves  propagating  in  different  direc¬ 
tions.  Parametric  interaction  of  complex  (non  planar)  waves  in  the  non¬ 
linear  material  can  be  understood  by  studying  the  interaction  of  their  angular 
spectrum  components.  The  Fourier  optics  formalism  puts  the  problem  in 
the  framework  of  linear  transverse -spatial  invariant  system. 

The  theory  of  parametric  upconversion  involving  three  plane 
waves  has  been  treated  in  considerable  detail^^’  Under  matching  condi¬ 

tions  the  plane  wave  signal  field  E.(z)  and  the  output  plane  wave  sum-fre- 

^  (1) 

quency  field  E  (z)  satisfy  the  relations 
s 

and 

E,(z)  =  E.(0)  cos 


where  the  characteristic  length  £  for  power  conversion 


2a), o)  d 
1  s 


C  (K.K  ) 
1  s 


172 


(56) 


with  K  the  wave  number,  o)  the  angular  temporal  frequency,  d  the  appropriate 
component  of  the  second-order  polarization  tensor,  and  E  (0)  the  pump  elec  - 
trie  field  at  input  surface  of  the  nonlinear  material.  In  the  limit  z  «  £,  the 
sum-frequency  field  is  given  by^^\ 
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12 


(57) 


E  (z)  ~  27rd 
s 


U) 

E  (0)  E  (0)  z 
C  1  p 


(23' 

As  shown  by  Kleinman'  ,  the  linear  growth  of  E  (z)is  ascribed  to  constructive 

interference  between  a  f3rced  wave  generated  by  the  polarization  and  a  free 

wave  generated  at  the  surface  by  the  forced  wave.  It  is  equally  valid,  as  we 

shall  see  later,  to  ascribe  the  growth  to  constructive  interference  between 

the  waves  radiated  from  each  differential  slab  of  polarization  in  the  material. 

We  developed  geometric  optics  theory  according  to  latter  view.  Eq.  (57) 

is  the  solution  for  the  matching  case  satisfying  E  (0)  =  0.  Actually  Eq.  (57) 

s 

does  not  give  the  correct  sum-frequency  field.  The  value  of  E  (0)  is  related 

s 

to  the  amplitude  of  nongrowing  waves  in  the  medium  and  a  reflected  wave 
at  the  surface^^^’  These  waves  are,  however,  many  order  of  magni¬ 

tude  weaker  than  the  growing  wave. 

The  rigorous  analysis  of  upconversion  process  is  based  upon  the 

(10,23,24) 

inhomogeneous  vector  wave  equation 

/  u)  \  2  /  (j)  \  2 

VXVXE^(r)  4;  E^(r)  =  47r  P  (l)  (58) 

where  ^  (x^^^  is  the  linear  susceptibility  tensor)  and  E^  are 

respectively  dielectric  constant  tensor  and  electric  field  at  the  frequency 

0/  =  CD,  +  o)  .  The  plane  polarization  wave  P  (r,  t)  is  related  to  the  signal 

S  1  {3 

electric  field  E.(r)  and  pump  field  E  (r)  by  the  relation 
-1  -  ^  p  - 

P  (r,t)  =  P  (r)  e“^‘^s^=d  ;  E.(r)  E  (r)  ^  (59) 

--  --  =-i - p- 
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For  non  planar  signal  and/or  pump  fields,  which  is  the  case  for 
imaging,  it  would  be  advantageous  to  solve  the  inhomogeneous  wave  equation 
for  the  angular  spectrum  component  waves  rather  than  solving  Eq.  (58) 
directly.  The  wave  equation  for  Fourier  components  wave  is  (Appendix  B) 

/a'  \  2 

V  X  V  X  6'  {k  ,  z)  -  €  :  (So  (x  ,  z)  =  4771  -^  )  (P  (k  .  z)  (60) 

-s  -s  =e  -s  \C/  “  -s 


It  is  apparent  that  operation  V  on  angular  spectrum  only  has  --  z^  operation. 
In  this  analysis  Cartesian  coordinates  x,  y,  z  are  defined  such  that  z  is  in 
the  propagation  direction  normal  to  the  surface  of  the  crystal,  and  the  optic 
axis  (consider  only  uniaxial  crystals)  -  is  in  the  xz  plane  (Fig.  11).  If  both 
signal  and  pump  fields  are  ordinary  wave  polarized  with  the  electric  vectors 
perpendicular  to  the  xz  plane  (along  the  y  direction)  the  effective  component 
of  the  polarization  wave  and  the  sum-frequency  electric  field  are  polarized 
along  the  x  axis.  For  this  simple  case  Eq.  (60)  becomes 
,2  \  /  u)  \  2 


(k  ,z)  =  4n(-^)  (p  (k  ,z) 
\  ^z^  e/  sx  -s  \  C  /  X  -s 


(61) 


where 


€  = 
e 


11 


=  d  6’ 


(k  ,z) 
ly  -s 


(See  Appendix  C) 


S  (k  ,  z) 

py  -s 


The  general  solution  of  Eq.  (60)  consists  of  an  inhomogeneous 
wave  (or  forced  wave)  which  is  some  convenient  particular  solution  of  in¬ 
homogeneous  equation,  and  suitable  free  waves  which  are  solutions  of  the 
homogeneous  equation  obtained  by  setting  ^  -  0.  The  free  waves  are  re- 
quired  so  as  to  satisfy  the  boundary  conditions  at  the  crystal  surface.  The 
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FIGURE  11. 


RELATION  BETWEEN  THE  COORDINATES  (x,  y,  z) 
AND  THE  CRYSTALLOGRAPHIC  AXIS  (X,  Y,  Z) 
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free  waves  have  the  form 


1/2 


z 


(62) 


where  the  refractive  index  n  is  a  function 

s 

direction  s  (unit  vector)  of  the  wave  vector 

<j) 


of  the  free  wave  propagation 


(63) 


The  relation  between  the  wave  vector  Kof  free  wave  and  the  wave  vector 

K  of  the  polarization  wave  P  (r)  will  be  considered  later.  The  polarization 
-s  -  s  — 

direction  01  satisfies  the  relation 


where 


:‘u=  0 


O) 


“s  = 


s 


(nI-€)-Kx 
2  9=  ='  -s-s 


(64) 

(65) 


where  I  is  the  unit  dyadic  and  €  is  the  dielectric  constant  tensor  at  sum -fre¬ 
quency  (j)  .  Equations  (64)  and  (65)  determine  both  n  and  the  polarization 
s  s 

direction  <11  .  We  define  a  quantity  n',  called  effective  index,  and  a  unit 
vector  a  by  the  relation 

U) 

Kg  =  -^n'a  (66) 

In  general  n'  ^  n  (a),  but  n'  =  n  (a)  along  the  phase -match  direction, 
s  —  s  - 

When  n'  ^  n  (a),  the  polarization  wave  is  said  to  be  mismatched, 

®  (23) 

the  mismatch  is  described  quantitatively  by  (f  in  the  relation 

nV  -  n  s  =  <^  N  (67) 

-  s—  — 
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or 


K  '  K  /  v*  N 
»  C- 


m) 


Mfirri'  ilta*  (ultii  %rt  i<ir)  mtfMMl  It*  iH«»  inrMirfif  fitirlurfi  oi  ltt«» 

tIrIrfMIlHr®  IfcMh  v*  Jtfltl  P  IM  llT  aiilKiriitliM  ITllh 

a  |»»ro»  »M^|  K  Hr»r  it  p  p)irr||)i  4  Ifipr  i*4%r  r%)iirfi#««t  Uy  ICq.  (62) 

a  p 

%lt|rh  fimai  rMai  Itt  I  hr  itir^rttrr  *4  lltr  fitrirti  »r4%  if  )h  urtkr  l«  MllUly 
Uiuntt4r>  niitiiiii<iit  4i  ilir  (rnri*Ml  iiurf4rt»  ttn?  iftHnituni  c4  pr#ii#nr^ 

«4  liiriH  iitfml  4M  irrr  '»4vr#  in  ihr  prailurfian  «<  KruKirinif  ufavr  lltmuith  in- 


irrlrirfurr  *4  4  trr**  nnvr  andttirrrHtrjivr 


(23) 


W)tm  V*  -  0  i)»r  wavtn  jirr  n^it) 


Itluftr  nuirltrvi  ami  ihr  uftmim*  in  itl%‘rn  by  (57).  ^nnnrnlly 

(23) 

umlrr  nrarl^  ntaicltiitt;  ritmliiiitnn  lollmirini;  KIrinnian  II  ranb(*ii)»Qvnlhtl 
ittr  iirtminy.  navr  in  glvnn  by 


»5  (*  .  /)  If  G(2  I  c)  >  :  (t*  (*  .  *) 
-n  n  -  -n 


wlicrr 


r«(;i) 


I  •  c 


4*  «*’**‘^ 


and 


(60) 


(70) 


2-'  ^  ^  n  Kj,  -  K  I 


Tlie  dyadlr  >  |s  diflmHl  bv 

«# 


(23) 


M.' 


)  -2-  i  |(N -Tl  )(>.‘U)  -  N  •  a  r^Oi  'U 

"  C(n  IV)  -  -  -  - 

s  i 


(71) 


(72) 
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Tit  i$  a|)|)rti\im4i|tHt  ii  tan  ttriiit  n 

V 

>  “  I  ,  ^ 

rtn  n  r  * 

«  I 

|tr«>V|t||;tl  Itairl  nral  Iv  IkutllMl  •«•  lltr  ogt'arr  !  »•»  Ihr  •  aor  ml  ^n 

uftlinJiry  |tllin|t  an  nftlinaf  )  olt^nal  «aVr,  aiMl  an  r%lta<>rt||iut v  »unt 

lm|nrnr)f  watr,  il>r  i;f«ti»int?  in  |4|  tCOl  Itmunr# 

»*  .#)  2f\  .  <f  i;»2  u  ./» 

t  <n  n  J 
n  I 

Ivquaiimiii  (69)  and  (74)  ran  l*r  niriiirn  in  m«rr  (ainillar  Inrm*  l»>‘  n*»iins.*  ilui 


C.2K)-r-'‘  "-S'' 

it 

11-2  KVALUATiQN  QK  MISMATCH  FACTOR 

llclorc  Kq.  (69)  can  be*  atifilird  it  ih  ncrcfinnry  in  ilu* 

nimmnlch  (iiclor  *.  In  termn  «(  wave  venom.  Fq.  (68)  |tivi*ii  lhal 


a' 

/  (K  ‘  K)  •  N  (70) 

C  -*«  - 

As  sho'  I  In  Fluurc  12(hewavr  vector  Kol  tree  waves  l.s  constrained  to  lie 

“’s 

on  the  wave  vector  surface  —  n  (s)  for  the  sum-frequency  w  ;  the  wave 

vector  K  of  polarization  wave  is  eiven  bv  K  Ki  •  K  and  i.s  arbitrarv. 

-s  8  *1) 

Wc  shall  denote  the  wave  vector  of  polarization  wave  in  the  matchiiv^  direc¬ 
tion  byK  •  Wc  shall  assume,  for  .simplicity,  that  the  surface  normal  N  is 
— M 

a  matching  diiection  and  also  that  K.. is  in  the  direction  N.  Figure  12 
shows  geometrically  the  mismatch  relation  expressed  by  Eq.  (76).  We  have 
taken  the  z  direction  to  be  along  N.  the  optic  axis  in  the  xz  plane, 
and  the  y  direction  normal  to  the  plane  of  figure.  The  tangent  plane  to  the 
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FIGURE  12.  GEOMETRIC  REPRESENTATION  OF  'MISMATCH 
RELATION 
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vt^rltir  ill  K..  li>  hIhiwiu  ami  ilit'  mtritiai  i<>  ihm  iilam*  tu.iKf  ilir 

{  w)ih  iIh^  4  atiiik.  ti  oluHihi  lir  mMml  ihai  ilir  tint'll^  {  rriirt  iiriiiii 
iuiim*(r*>|iv  III  ili*»  rr^iiial  ai  ih*’  Irrtiurwry  tli*'  4i»iiM»irii|iv  aH^lr  i  in 

Kivtrtt  Uy  (Af>|irmiu  t)). 

3  /  I  *  \ 

lim  {  -  n  I**  I  nm  *►  mn  **  I  J  ■  i  j  iTJ) 

m  m  m\«  tt  *  f 

u  r 


irhtrf#  #  -  Afliilr  ti»1iir«>«>n  Ihr  ttfilir  aailn  amt  Ihr  jihaM'  tiulrh  ilitrrlim 

fW 

filii  I  -  itide^  til  r«ilrart|i4t  alumi  ai  * 

m  -  M  # 

n  -  0itfrm»r«lin»fy*»i»¥t>  imlr%  ol  frlr4nii4i  41  « 

p  • 

n  -  af(lin4ry*«r4V»  iml#\  »!  rflrartum  41  •* 

o  p 

Thin  aotln  alfo cftvrrnii  thr douttlr  rHrarliofirflrri.  in  tkhirh ittai h  ihr  ortlinar ^ 
and  rxtraordinary  aravmi  arr  ai  iKr  namn If M|upnry.  In thi*  ana)y«.«  11  iraa  an* 
aumad  lhal  Ihr  nifinal  llrld  and  Ihn  iminpllrld  arr  Imlh  fmlar Urd  aa  ordinary  aravra 
andlhpaum-lroqunnt*ylip|d  la  polarixrdaa  cxlranrdlnary  aavra.  The  mlamairh 
lacior  V  arrordliiK  lo  Kq  176)  la  Ih©  dialanr©  from  lo  Ih©  wav©  vreiar 
iurlac*  meaiured  ftarallel  lo  N.  Thia  la  afiparvnl  from  the  boundary  coitdi* 
lion  lhal  requires  conlinuily  of  iransvers  wave  verinra  of  both  forred  and 
free  waves  across  ihe  boundary  surface. 


K  •  K  '  *  ♦  a  (70) 

-s  -p  -I 

This  condilion  assures  also  the  conlinuily  of  an  angular  »i)eclra  with  spalial 

an{(ular  frequency  k  across  the  crystal  boundiry.  It  is  somewhat  icilious 

""  s 

but  can  be  shown  lhal  (Appendix  D)  ^  _ 


U.’ 

S 


C 


K 

sz 


K  tan  C 
sx 


K 

2 


M  sx 


K  C 
sv 


(79) 


yie- 
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hi*  1 1- 


4I«>| 


*  *  4  >  «.  ^ 

a  11  ny  II 

*»  ^  n 

K  '  tK  *  •  *  *) 

»#  «  • 


jiiMrr  ihr  il  I*  igniMlh  «•r4rlv  anti  ihr  quanUly  in  Ihn  ftnuarn 

hrarkri  I*  arruntl xtrtirr  in  (K^  •  K)  an»l  nttMlIrr  lh4n  ll»n  linraf  t|U^nllfy.  II 
I#  |«r«t|trr  lit  itur  •implt^i  (<trtn 


ft 

-j  • 

C 


♦  K 


« 


(101 


\lr  lu%r  n»r<|  fan  «  -  (  in  lin  (10)  mnrr  (  I. 

U-3  Hi;i*l*LTAN*T  f  MrQt  r*  CY  HHU) 

Thn  iium*(rrqtimcy  ani^Ur  n|»erirom  fitm|w»rni  JS>  t)  uivrn 
by  Kq.  (69)  in  itinmirtl  by  a  iiMnr  nfjive  imlaruallon  £  (»^.  ().  which  in  a 
Kinirinr  rrtmniirK*ni  ol  inial  |>tilarU4iicm  W4%‘«*  P(i).  Thr  renullnnl  num-lre* 
qucnry  (irlH  privlurril  by  P(r)  lit  Fourirr  Inimtiorm  «l  6  («  .  it) 

Ml  M  -  II  <M  II 


K  (r)  I  ^  U  ,  !:)♦**•*♦  ’  (81) 

H  -  It  -It  “1* 

SubHilluliiti:  Kq.  (69)  into  Kq.  (81)  Ihr  rcnullnnl  (irltl  can  lie  written 

(r)  ^  Jn  J  >  •  (Ik  (82) 

where  unc  h.i.s  liorii  made  of  intoftral  reiircscntation  In  Kq.  (70). 

2 

For  paraxial  as.sumption  the  qundraic  term  «  /2K,,  can  lx* 

s  M 

iicclcctod.  It  then  follow.s  from  Fqs.  (82)  ;ind  inverse  Fourier  integral 
representation  of  iI’(k.  z)  that 
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*■  (a,i*  fo  •*'>/•'•.'  *  '  «  ' 


..  -**  *»» 


/  il  |»*  >  IMf  *J  «>  # 


TV  initcmilon  mrt>r  •.  ihi>  («riur 

6  U  •**•<»<#)&  <y  •  y*) 

TV  miftr«iicitt  o¥#r  (»*  £.!¥## 

iMr)  '  >  /q  i  ^  P  t%  •  »i»f#.y.  #)  (§4) 

Then  the  held  ni  the  eiui  tturfatf  ol  the  nonlinetr  *  rynial  («  -  tP  ran  he 
wrillen 

tjlp.  0)  -  D  /q  d  I)  e’***^'M^  ^  Pin .  »kD.  y.  1)1  tW) 

U  we  let  (I  -  f))  D  *  «  and  moke  a  at  a  dummy  variable  of  infeyralion  (IS) 
becomes 

E  (p.  D)  •  di  e**^M‘  >  :  P(  a  •  (O  .  (  a.  y.  D)  (M) 

-g  ^  'Q  «• 

This  expression  shows  that  the  resultant  sum ’frequency  field  Is 

produced  by  constructive  Interference  of  the  w&^es  radiated  from  each  dif- 

fcrential  slab  of  polarisation  In  the  material.  It  Is  this  concept  that  was  used 

in  the  geometric  optics  theory.  Double  refraction  (or  anisotropy  effect)  of  the 

sum  frequency  electric  field  Is  included  In  Eq.  (86)  through  the  parameter 

(x  -  £  D  4  C  z).  Now  refer  to  Figure  13  the  line  at  angle  L,  which  is  normal 

to  the  tangential  plane  to  the  wave  vector  surface,  represents  the  direction 

(27) 

of  energy  flow  of  the  sum-frequtncy  field  .  It  is  intuitive  to  expect  that 
the  path  of  integration  to  obtain  resultant  field  should  be  along  the  energy 
flow  direction  rather  than  the  z  direction.  To  understand  the  implication  of 
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\  lAAlCf  l«tiA4,  (hlAIWl 


OiRECTlOH  OF  FOYNTihC 
VECTOR  OF  EXTRAOROtNARV 
SUM  frequency  WAVE 


FIGURE  13.  ILLUSTRATION  OF  DOUBLE  REFRACTION  OF  SUM- 
FREQUENCY  WAVE 
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lilt*  %  •  {  M  •  t  /  |ttl  M»I  ((riiotr  |ki«|ltlo  lit  lljf  {•••I.M  UulliMi  »  .iVf 

|Ik$  fi«turt'p  Itmintii,  liy  K^,  V^,  ami  |*ii|ttla  lit  lltr  ■•  um  Irr^ltinu  %  m.i^r 

Ihtf  t4»fi<;rVtTr  |*>  V,.  fri;^#il  I  q  ‘IMjI  40  la  ia!r|;t4l 

iHi-rr  amiirr  |»i|nlii  {flli  iaf;  Ihc  iiMiii  <if rqMriir  V  lit  l»l  al  4t»  •■Iturr  %  «'»  |>*<»lMl  Vtv 
pImII  tflcafHitr  lltr  Miyrrr  jxttalw  Ifttm  Ihr  ri«it*|)iK4lr  oruta  \%t*  •Irllt^r 
far  ihr  tiJ|»*rfirrr  imiltil*  la  •lirli  a  aa|r  Hial  ^  ®*  '***'*  ****  ****'  '^1***' 

»  r 

m^Mfll^  ait^lr  (  Hrllh  ||*r  t  4\|ft  tWn  ami  ll»r  |*ulli  tj  |aie^r4l|«M 

•  « 

trtfrf  ihr  iMturre  mtinia  l»»  i»|»rc|l|r*l  liy  ^  f ,  >.  v  »|ih 

I  ?  I  I  ?  ? 

hnlH  r«in#i4tti  in  ihr  i  tir^r’aliait  ntnirr  |i  irirr*  lltr  (>|*i>r rt  4i|>>«t  |»<iini  tHra 
irhrn  nrr  hipI  (  -  I),  trt  r4il4in  Ihr  iniri^f  aiKtn  |Mia*i»rir»  % 

Thin  ran  hr  inirrutrriml  ihai  ihr  nuin  >lrr<|urnrv  lirhl  4i  iltr  irtinu  1%.  v.DI 
in  utii4inr<l  hy  inireriiiltin  ahme  ihr  mrrisy  Ihm  linr. 


Thr  Aninil^r  n|>rrifam  ol  ihr  nuni*l»r<|umrv  lirhl  41  ihr  r\ii 


nurfarr  ol  ihr  crynial  in  nlmply 


*  .  0)  f)  r  2  4  4*  U  .  D) 

n  n  -  n 


"  (1") 


whore 


uK  K  -  K  . ( ^ 

n/  n  nx  2K. 


Rq.  (87)  indicalen  lh;il  Ihc  ani:ular  H|H'rlrum  al  llu*  exit  surf  are  Is  nol  a  sum¬ 
mation  over  (he  contributions  o(  the  Fourier  com|N)nenls  of  the  iKilari/alion 

(13) 

in  fhe  differential  slab.  FIrcstor  ralruinti'd  the  ani:ular  s|H*clrum  at  the 
exit  surface  by  summini:  over  the  anpilar  spectrum  of  same  spaliai  frequency 
radiate  i  in  each  differential  slab.  By  doing  so  he  olilained  the  faiMor 


/JL 

'  4  V  K.  ■  K  /  /  4  VK  ,■  K  /  ' 

L  1  S  -J/  L  i  S  -J 


md  considered  it  as  a 
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f*  iti  iiiiMii.tii  li.  •'tMili'ary  lu  Inn  iitilial  iin* 

<•!  Ml  |tlur>(;  ihmIi  M  I'tuMluiiui.  IW'iMuno  u(  iitin  li#  Mrrivetl 

mI  iitt'tx  rrt  I  lot  III  lltr  ii|i)rrl  Mini  littt 

III.  AMii  I  All  M‘i  t  rill  M  or  iMAtii,  nt  i  J)  mil  m  vKHAi.iMAur 

l‘l•<'u^\  I  ll^Kr^•  m.h 

till  fr»4tla  ifUjItirvI  III  |trrt'r«ilti^  arriluti  Mill  |jip  «Brt»  |u  rr«|i|«lr 
lli(r  altj^wlMf  efirrltYlltl  ol  Ihe  *11111 « | r<r«)t|rytry  ItlM^r  ||r|ll  luf  «|*rtt»|lit;rtrr 

•yairtii*  |i  Bhiiiilil  l»r  lanininl  u»i*  hrfr  1H4I  I**).  i60l  Inlmi  l#|lw  arruutii  t*| 
tt«|i>4ttMirh  4tMl  MtilM4ru|*|i  |i  1*111  l*r  BNonn  Ikal  Ihr  limit?  liHrrnrlimi  !i|?rr- 
lufr  tiiMy  U  im  lwlc*|  iti  Ihr  ralniUlIxiti  l<y  KiirtMlurliiig  n|»rrlurr  |utirl|a(i> 

|l  «  |||  |*r  nlw^tl  !l|  IhlB  Brel  linn  IHmI  Ihr  imA^r  ^«bII|u4|  f  rljlIMl  tlcrlml  lit  Ihr 

igruittrfrir  uftiirB  ihn*f'y  I*  uhRMiiird  Iimhi  M4iH*«ury  |ihm»r  rumnd^nilt^  To 
rtiiti|»Mrr  Ihr  friiiilin  »!  ihifi  nrrinio  oilh  IhnI  u|  igrottiririr  oftHm  Ihrvtry.  H 
nap  Mnnumrtl  1I14I  ihr  mmlinrar  mmrrinl  fllU  ihe  n^tjirr  Umttdrd  hy  {tnnillrl 
Itlmiri*  /  •M2  and  /  |t  2  iKli;  I3t.  ll  *sjt  alwt  anntitttrti  Ihm  all  llrW 

ttuaitiiUrM  tuvr  drlimic  inilari/a'ion.  nlnr?*  |Mramnrir  uftrfmvfrnion  in  a 
nttniinrar  rr>‘nial  Kfnrr.ill>  rof|ulrr»  drfinlir  itularMnllon  lor  phanr  malrhlnc. 
For  mniplimv  wr  Mitaln  roimidrr  ihai  ihr  infrared  ohim  and  ihr  itump  wavra 
arc  ptilari/cti  an  ordinary  wa^'rn  and  ihc  itunt*frcc|urnry  lirlda  arc  |ailarixrd 
an  c  xlraordln.ir>'  wavr.  For  ihin  ranr  Ihr  Held  and  Ihc  lailarirailon  wave  ran 
Ih'  dcnrrllH*<l  lt>*  nrahir  t|uanllllcn. 

II I  - 1  PI.ANK  WAVF  PUMP  SYSTK  M 

For  a  phitu*  puiup  wave  proftaiial tin:  In  Ihc  dircrilon  ol  Ihc  axin 

IK  7 

K  (r)  K  c  P  (80) 

-|)  -p 

its  angular  sped  rum  is  simply  a  Iwu*dimcnsional  della  fund  ion  multiplied 
by  a  phase  fac  lor 

£  (k  ,z).  E„cV«(‘  (90) 

-P -P  P  -P 
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IniiUlt;  Ihr  |iuim|)  itiHtulitr  nimMrufti 

t  f*  c  3  *  K  •  ’i) I 

*  («.  ,  /I  -  T  5  61*  I#  ^  2  u  2 


(on 


*h»trr 


I)  .  0 
2  1 


K  ~  Krjivr  ntttnWr 

^  P 

^ —  '  «rjiw»  fH»mt«rr  in  IH*»  cftMiil 

P  * 

P 


•n*} 

T  ’  irAAiiititffikm  rortllrirfii  kir  iHp  pump 

p  r-  *• 

«iir«>  •!  ih»  mairr lAt  Mirfare 

IIk^  Q(t)fcf  (lH4  durtrlbyliofi  trrtiMi  Ihr  |4anr  «  -  «  » 
br  dnried  bjr  K^(p,  -  #).  ai^Ur  >pe<inim  inaHk  ihr  nonllnr^r 

rryiMal  la  eimi  ti)r 


*’  *  V|t|- 


I  » 


irhrr* 


<  »  < 


I) 

2 


(02) 


SU. 

1^1 


a) 


iifTi 


/k,Ip. 


•aVp 


K 


I 


2» 

* 


Irrr  iifMirr  «ravr  mimbrr 
*  watt*  mimbrr  inauir  ihr  rryMal 


(aJ) 


% 


iransmisHion  coefflclont  for  the  angular 

spectrum  at  the  material  surface. 

Transmission  coefficient  T,  depends  on  the  vector  k.  ,  since  k, 

I  -i  -i 

determines  the  incidence  angle  of  a  spectral  component  and  transmission 
coefficient  depends  on  the  incidence  angle.  For  the  paraxial  assumption, 
which  is  valid  for  the  most  of  upconverters,  we  can  neglect  the  k  .  dependence 
of  triuismission  coefficient  and  approximate  it  by  a  constant. 

The  angular  spectrum  of  polarization  wave,  convolution  of 
S  (k  ,  z)  and  5.(k  z),  is  given  by 

p-p  i  -i 


for  -  D  •  z  -  D 
2  2 

-1  CK  ^  -  K  (z  +  5.)) 

Ofe  .  z)  -  dT  TS  S.(i  ,  -  s)e  Pa2  p  2 

~S  p  i  p  I  s 

2  2  1/2,  D,  _2  2.1/2  j, 

ll(K|^  )  (s-^)+(Kj  -K.  )  (z+-)) 

e 


(94) 


(P  (k  ,  z)  =  0 
-s 


for  -  —  >  z  >  — 
2  2 


From  phase  ma*ch  condition  K  =  K.  +  K  and  a  =  k  .  +  k  and  K  =  K  z 

sz  iz  pz  s  ~i  ~p  — p  p  0 

it  can  be  shown 

K  =K 

sz  p  I  -s  s  s  s 


where 


K  =  K  i  +  K 
— s  sz  o  -s 


Substituting  this  relation,  Eq.  (94)  can  be  rewritten  as 

D  2  2.  1/2  .  Dx  2  xl/2D. 

iK  -5-  Jl.(K:„  -K  )  (z+=-)-(K.„  -K  )  -^J 


(P  (k  ,  z)  =dT  T.(5  S.(k  ,  -  s) 
-s  pip  i-s 


e  pa  2  e  ia  -s 


la  -pa 


2  2  1/2 
i(Kia  -  Ks  )  s 


D  „  D 

fni*  •  “““  ^  Z  ^  — 

tor  2  2 


(95) 
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The  angular  spectrum  of  the  sum-frequency  field  is  obtained  from  Eqs.  (74), 
(75)  and  (95), 


2nco  dT  T.S 

&  (k  ,z)  =  i  C(n  n  )1/2  i-s  2 

s  ~s  si 


Sin 

Ak(z  +  y) 

2 

Ak(z  +  y) 

L  2  J 

xe 


AK(z+?-  .r,  2  2,^/2  D  ^  D  2  2.1/2  2  2  1/2 

^  2  -if  K.  -K  )  -z- +K  —  -  {K.  -K  )  s-(K  -K^)  (z  + 

—  '■  ia  s  2  pa  2  la  -s  s  s 

e 


^ - T 


D  „  D 
. —  ■<  z  ^  — 

2  2 


(96) 


where  AK  is  given  by  Eq.  (88) 

The  angular  spectrum  of  upconverted  filed  beyond  the  exit  surface  of  crystal 
is  given  by 


n 

S  {k  ,  z)  =  T  S  {k  ,  —  )e  ^ 
s -6  s  s -s  2 


1/2 


-  K 


sa 


(z-y) 


Z  >5. 
2 


(97) 


where 


Tg  =  transmission  coefficient  at  the  exist  surface 

K  -  free  space  wave  number  at  w 

sa  X  s 

s 

It  follows  then  from  Eqs.  (96)  and  (97) 
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m|0 


SinC"^-^)  27r.a»  T  T  T  d5  D 

n/  .nr  .  2  2  ISSipP 

^  “t'-s  ■  ■  C(n  n.)l/2 

o  Si 


2  2  1/2  n  D  2  2  1/2  2  2  1/2  D 

-iUK.  )  (^  -  s)  .  K  y  -  (K  -  «  )  D  -  (K  -  O  (z-^  ] 

la  s  2  pa  2  s  s  sa  s  2  " 


•e 


where 


(98) 


D 

/  .  IK  •  p  iK.  z, 

For  the  case  of  a  plane  object  wave  E^r )  =  E^e  -  -  e  iz 
its  angular  spectrum  are 

S,{h  ,  -  s)  =  S .6  (k  -  k) 
l-s  I  S  - 

The  resultant  sum -frequency  field,  which  can  be  calculated  by 
taking  Fourier  transform  of  Eq  (98),  becomes 


E  (p,  z)  =  i 
s  - 


27rw  dT.T  T  E  E.D  Sin(^^)  .AKD 
s  1  p  s  p  1  2  -1 

C(n  ■■■■■■  ^  ^ 

SI  2 


.  K  D  .  (k/  -  .  (kJ  .  -  1)1 

•  0  A 

for  z  > -^ 


(99) 


The  sum -frequency  field  is  also  a  plane  wave  and  its  amplitude 
is  reduced  by  the  mismatch  function  Sin  ( )  /(— —^  ).  For  plane  object 

U  u 

wave  propagating  along  the  z  axis  k  is  equal  to  zero  and  the  mismatch  function 
becomes  unity  because  of  AK  =  0,  and  Eq.  (97)  reduces  to  Eq.  (57)  if  we  ignore 
the  phase  change  expressed  as  exponential  in  Eq.  (99). 
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We  shall  next  consider  the  relation  of  inia^!,e  location  to  the  object 
location  by  investigating'  the  phase  factor  of  Eq,  (98).  Using  Eq.  (88)  foi 
AK  and  the  approximation 


2  2.  1/ 

-  K  )  ~  K 


2 

K 

2K 


which  is  valid  for  paraxial  assumption,  the  phase  factor  of  S  {k  ,  z)  is  given 

s  s 

by 

2  2  2 

(^sa  - 1  )  -  ^^ia  -2^-^^ -  ^  ^ 

sa  la  s 


1 

-  [K  -  K  ^  +  tK  +  -^(-— •  )|~ 

s  M  -sx  2  K  K  '  2 

M  s 

Since  the  constant  phase  factor  that  is  independent  of  x  is  common  to  all  the 
angular  spectrum  and  represent  a  constant  phase  shift,  it  can  be  ignored. 
The  phase  factor  that  depends  on  k  is  then 


D  Xs  ,  z 


D 


D 


D 

+  — 


^(x  )  =  -  fgx'  — +  ^^ (• 

-S  '’-SX2  2  K  2K  2K.  2K 

sa  sa  la  M 


k.  2K 
la  s 


(100) 


Optimum  focusing  occurs  when  4>(x  )is  identically  equal  to  zero  for  all  x  . 

“s  ~s 

If  this  condition  is  satisfied,  the  phase  angle  of  all  angular  spectrum  becomes 

identically  equal  and  the  resultant  field  at  the  sum  frequency  will  form  a 

clearly  defined  image  across  a  plane  z  =  -  s'  which  makes  $  ( x  )  =  0. 

— s 

Unfortunately,  the  condition  4>(x  )  -  0  for  all  x  can  not  be  satis- 

s  s 

fied,  since  we  have  one  equation  and  two  independent  variable  x  and  x  . 

“SX  “sy 

There  is  no  unique  solution  z  -  s'  that  makes  tl’(x  )  =  0  for  all  x  . 

s  -s 

There  is,  however,  a  solution  z  =  -  s  'for  a  particular  x  .  This  can  be  in- 

s 

terpreted  that  different  components  of  angular  spectrum  focus  on  different 
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transverse  plane.  The  sum  frequency  field  produces,  therefore,  a  blured 
ima^e.  The  origins  of  this  image  degradation  are  crystal  anistropy  and 
phase  mismatch.  It  is  natural  to  expect  that  the  upconverter  is  not  circular 
symmetric  because  of  crystal  anisotropy. 


When  the  linear  phase  factor  ^  k  D  is  ignored,  which  is  valid 

-sx 

assumption  for  the  most  of  upconverters  because  of  paraxial  assumption  and 

small  anistropy,  there  is  a  unique  solution  z  -  -  s'  that  makes  )  =  0  for 

s 

all  «:  =0  for  all  k  provided  K  IK.  +  K  I  is  a  constant, 

-s  -s  s  -1  -p 


K 


s  =  s 


sa 


DK 


K. 


sa  .  1  1 

i-rrr  + 


la 


s  M 


--- ) 
K.  K  ' 
la  sa 


(101) 


For  infinitesimally  thin  crystal,  D  ~  0  and  Eq.  (101)  reduces  to 


s'  = 

1  s 


(102) 


This  equation  is  identical  to  the  result  obtained  by  the  geometric  optics  theory. 
(Eq.  (6)  of  Chapter  2  with  p  =  «). 

The  assumption  K  =  constant  is  not  valid  since  K  is  related  to 
s  s 

-i  ^  ~  -s^  by  the  equation 

K^  =  x^+(K  +K)^ 

G  s  iz  p 

When  both  pump  and  object  waves  are  ordinary  waves,  it  is  true 
2  2  2  2  2 

that  K.  +  K  {=  K  .  +  K.  +  K  )  is  a  constant,  but  K  is  obviously  not  a 
i  p  -I  iz  p  ’  s  ^ 

constant.  The  unique  solution  z  =  -  s'  that  makes  ‘J>(x  )  =  0  for  all  k  does, 

s  s 

therefore,  not  exist  even  for  an  upconverter  of  isotropic  crystal.  It  is  clear 
iiow  that  thickness  aberration  introduced  in  the  geometric  optics  theory  is  due 
to  non-uniqueness  in  the  position  of  stationary  phase  plane  and  amplitude  re¬ 
duction  of  angular  spectrum  of  the  image  field  caused  by  the  mismatch. 
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The  angular  aperture  of  the  upconverter  can  be  calculated 
from  the  condition  Sin  =  0.  This  relation  will  give  the  extreme  values 

for  X  and  consequently  the  maximum  object  wave  incidence  angle.  This 
s 

procedure  is  carried  out  in  the  geometric  optics  theory  and  the  *  suit  is 
presented  in  Sec.  IV- 2,  Chapter  2. 

The  amount  of  longitudinal  spread  of  upconverted  image  may  be 
evaluated  from  Eq.  (101)  by  substituting  extreme  values  of  obtained  from 
the  condition  Sin(^^)  =  0,  The  electric  field  distribution  of  sum-frequency 
image  of  a  point  object  located  at  z  =  -  s  will  be  considered.  We  assume 
that  crystal  anisotropy  is  small  so  that  £  =  0.  The  image  field  distribution 
will  be  evaluated  over  the  plane  z  =  -  s'  given  by  Eq.  (101).  The  angular 
spectrum  S.(k,  -  s)  of  a  point  object  field  E.6(p  -  p.)  is 

1  _  j 

ix  •£ . 

S.(k  ,  -  s)  =  <?.e  ®  '  (103) 

1  -s  1 


We  set  p .  =  0  for  simplicity. 


The  angular  spectrum  over  the  plane  z  =  -  s'  is,  neglecting  constant  phase 
factor,  then 


S  (x  , 
s  -s 


-  s')  =  <5 


Sln(^) 

1m 


(104) 


where 


S  -  i 
s 

AK  ^ 


27ru}  T  T.T  d<§  SD 
s  s  ip  p  i 

c(n  n.)l72 

2  ®  ^ 

X 

■— —  (See  Appendix  F) 


1 

The  image  field  is  Fourier  transform  of  Eq.  (104) 


(105) 
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E 


(P.  -  SM  =  -  5 


ifc 

-s  -s  . 
e  (Jk 


Using  polar  coordinate  expression  k  =  k  and  p  =  pe^^,  we  can  convert 

—  s  s  — 

above  integration  to  Hankel  transform 

2. 

Sinf 

E  (p,  -s')=-27r5  /  - S - (xpkdx  (106) 


(jlR\ 

’'”\4K.  / 

m 


where  J  is  zeroth  order  Bessel  function.  E  (p,  -  s)  can  be  expressed  in  the 
o  s 

form  (Appendix  E) 


477K  5  K  p^ 

E^(p,  -s') - ^[£.Si(-^)l  _  (107) 

where  sine  integral  function  is  given  by 

.  /-x  sin  t 

Si(x)  =  j  — —  dt 

^  o  t 

KiP^ 

The  function  [  r  -  Si(  ~ — )|  has  the  same  character  as  the  fimction  describing 

u  L/ 

difiraction  from  a  round  hole.  The  image  field  distribution  E  (p,  -  s'), 

s 

consequently  the  extent  of  aberration,  is  independent  of  the  object  location. 
This  result  agrees  with  that  of  the  geometric  optics  theory. 
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It  should  be  noted  that  the  spread  of  sum-frequeney  iiDa^o  in  Kq.  (107) 

is  caused  by  thickness  aberration.  The  spread  is  not  due  to  a])erture  dif¬ 
fraction,  since  wo  assumed  an  infinitely  wide  crystal.  It  is  easily  seen  from 
Eq.  (107)  that  the  field  spread  becomes  wider  for  the  tliickor  nonlinear 
material.  It  also  can  be  shown  that  Ey(p,  -  s')  approaches  to  the  delta  func¬ 
tion  as  the  material  thickness  D  approaches  zero. 

The  sum -frequency  electric  field  expressed  in  Eq.  (107)  is  a 
point  spread  function  of  the  system  under  consideration.  The  minimum  re¬ 
solvable  image  spot  size  according  to  Rayleigh  criteria  may  be  calculal^d  by 

TT 

noting  that  the  function  g  “  Si(x)  possesses  first  zero  at  x  2.  The  radius 
of  resolvable  spot  size  is  then 


1/2 

p  =1.4(D/K.) 
res.  1 


(108) 


ni-2  SPHERICAL  WAVE  PUMP  SYSTEM 

We  shall  consider  the  upconverter  pumped  by  spherical  waves. 
A  point  pump  source  that  radiates  spherical  wave  is  located  at  the  point 
(x  ,  y^,  -  p).  The  pump  field  distribution  across  the  plane  z  =  -  p  is  ex¬ 
pressed  by 

E  ip,  -  p)  =  E  a(x-x  )6(y-y  ) 

p  ~  p  p  p 


The  angular  spectrum  of  this  point  pump  source  is  then 


S  {k  , 

p  -p 


-IK  •  p 

-p)-S  e 

P 
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where 


’  i^p/(2'rrr  P  ^  X  X  +  y  y 

p  -p  p  o  p  o 


We  assume  a^ain  that  an  object  is  located  on  the  plane  z  =  -  s  and  its  angular 
spectrum  is  given  by  "  s).  The  angular  spectral  of  pump  field  and 

object  field  inside  the  crystal  arc  respectively 

o  p  1/2  -pv  9  9  1/2  n 

-i[K  -0  -(k:“  -X  )  (P-|)-(K  (z  +  ^)] 

^  .c  n  -P  -P  P^  P  2  p  p  2  ■* 


S  (k  ,  z)  =  T  (S  e 
P“P  P  P 


u  .(x  . ,  z)  ^  T  o  \  -  s)  e 

1-1  i  i  -1 


,  D  D 

for  -  TT^  z  ^  — 

2  2 

(109) 

i[(K./-x.  )  (s  -  .  (k/  -  x.^)  (z.y)] 

la  i  z  i  i  z 

- —  ^  z  ^  — 

2  2 


The  angular  spec^um  of  polarization  wave,  calculated  again  by  taking 

convolution  of  5  (x  ,  z)  and  (S.(x  . ,  z)  is  given  by 
P“P  1-1 

9  9  n  D  2  9  1/2  n 

i[(K  -X  )(z  +  — )  +  K.  (s -- ) +(K  -  X  )  (p-^)-x  -P  ] 

<H(.  ,  z)  =  dT  T  5  e  "  2  P**  ®  2  -P 

s  1  p  p 

,2  ®  ■  y  p  ■  J  P  ■  f 

■  f  Ak^^SAk',  -s)e  pa  pa 


2  ^  2 


(P{k  ,  z)  -  0 
s 


D  D 

—  >  z  >  — 

2  2 


(110) 
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2.1/2 


The  phase  match  condition  (K  -  k  )  +  (K.  -  k  .  ) 

p  p  1  1 


2,1/2 


-  (K  “  - 


2,1/2 


.2  2,1/2  _  K‘ 


and  the  approximation  (K  -  k  )  ~  K  -  —  are  used  to  arrive  at  Eq.  (110). 

2iv 

The  angular  spectrum  of  the  sum -frequency  field  beyond  the  exit  surface  of 
crystal  is 

2tw  T.T  T  d5  D  sin 

g(s,z)  =  i  - 

s  -s  c(n  n.)  1/^  ,  AKD , 

s  1  {  — ) 

2  X  ^  X  ^ 

i[(K^-^)D  .  K.^(S-|)  .  (K  -  ^  )  (p  -  p  )(z-5-)l 

•  e  2K  pa  sa 

s  ^ 


D  D 


D 


•  f  dK  'S .(k  ' ,  -s)e 
-  1  - 


2  ®“2  P'T  P”? 

pa 


2K.  2K 
la  pa 


D 

"  ^2 


(111) 


Above  expression  may  be  evaluated  analytically  for  a  few  simple 
object  angular  spectrum  S  (x  -  s).  The  most  interesting  case  is  one  in 

S  1 

which 


-IX  .  •  p. 

S  (k  . ,  -s)  =  S .  e  '  ' 

1  -1  i 

Above  is  angular  spectrum  of  a  point  object  located  at  the  point  (p  -  s). 
The  integral 
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can  be  calculated  by  using  the  relation 


2  . 


ibx"  -iax  .  ,  Tr_  (a^/4b) 
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We  let  x^  =  y^  =  0  for  simplicity  in  following  calculation.  The  angular 
spectrum  of  image  field  produced  by  a  point  object  and  a  point  pump  becomes 


S  {k  ,z)  =  S  expi 
s  -s  s 
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where 

4770)  T.T  T  dDS  S 

S  =  -  —-g._L.P_s _ i_2 
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la  pa 


s.n(  — ) 

(m 


(114) 
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The  phase  factor  that  dependH  on  n  ^  4s  found  to  In* 


If  we  neglect  anistropy  factor  f:  and  let  *  Vj  =  0,  we  obtain  an  equation 
relating  the  image  location  (z  =  -  s')  to  the  positions  of  point  pump  and  point 
object. 


K 

sa 


2 


=  0 


where  z  =  -  s'  makes  4>(k  )  =  0 

—  s 

For  an  infinitesimally  thin  material  D  =  0  and  we  obtain 


(116) 


71 


(117) 


K  K.  K 
sa  la  _  pa 

s'  ‘  s  ~  p 


Since  K.  =  n.K.  (where  j  =  s,  i,  p) 

]  ]  ]a 

We  may  rewrite  Eq.  (117)  in  the  form 


K  K.  K 
s  1  p 


s'n  sn,  pn 

s  1  P 

and  recognizing  /3  =  Kj/k  and  Kn/K  =  I-/3  we  obtain 

s  s 


1  ^  _  1-/3 


s'n  sn.  pn 
s  1  p 

(118) 


This  relation  is  identical  to  Eq.  (6)  of  Chapter  2. 

Longitudinal  spread  of  the  upconverted  image  by  anistropic  nonlinear  crystal 

with  finite  thickness  is  again  apparent  from  lack  of  unique  solution  z  =  -  s' 

that  makes  Eq.  (115)  identically  zero  for  all  k  .  The  effect  of  off-axis  waves 

s 

is  accounted  by  terms  which  are  linear  in  the  transverse  object  vector  p .  in 

Eq.  (115).  This  effect  leads  for  thickness  coma  of  upconverted  image. 

If  £  =  0  and  K  =  K,.  =  constant,  a  good  approximation  for  near  isotropic 
s  M 

material,  Eq.  (116)  can  be  reduced  to 


1 


1 


/  /  D  , 

(s  +  -r)  n 
2  s 


D 


(119) 
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Above  equation  is  equal  to  Eq.  (22)  of  Chapter  2  if  we  set  t  =  — 

in  latter  equation.  It  is  tempting  to  conclude  from  the  result  that  a  point 

object  at  z  =  -  s  produces  a  point  image  at  the  sum -frequency  at  the  point 

z  =  -  s'  if  crystal  anisotropy  and  mismatch  are  negligibly  small.  This 

conclusion  is  in  error  because  the  condition  K  ==  K,,  =  constant  does  not 

imply  the  matching  condition  AK=0.  Under  the  matched  isotropic  operating 

condition,  uncertainty  in  the  location  of  stationary  phase  plane  is  introduced 

through  K  dependence  of  K  .  In  the  geometric  optics  theory,  we  attributed 
s  s 

thickness  aberration  as  a  consequency  of  interference  of  sum -frequency  fields 
radiated  from  each  differential  slab  of  polarization  wave.  It  is,  however, 
shown  in  this  chapter  that  the  real  causes  of  thickness  aberration  are  mis¬ 
match  and  crystal  anistropy. 

We  shall  next  evaluate  the  image  field  of  a  on -axis  point  object 
under  the  assumption  S  =0.  We  also  assume  that  Eq.  (119)  is  the  stationary 
plane.  The  angular  spectrum  of  im^e  field  then  follows  from  Eq.  (114) 
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(See  appendix  F) 

The  electric  field  distribution  of  image  field  is  given  by  (Appendix  E) 
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It  is  apparent  that  the  point  spread  function  Eq.  (122)  depends  on  location  of 
the  point  pump  and  the  point  object.  This  result  agrees  with  the  result  of  the 
geometric  optics  theory. 
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Ill-  3  GAUSSIAN  BEAM  PUMPED  SYSTEM 


We  consider  the  system  in  which  the  nonlinear  crystal  is 
placed  in  the  near  field  of  the  fundamental  TEMqq  mode  of  a  laser  output 
beam.  The  electric  field  distribution  of  this  beam  is  Gaussian  and  the  dif¬ 
fraction  is  minimum.  According  to  theory  the  radial  distribution  of  electric 

f281 

field  of  the  fundamental  mode  is  ' 


E(p,z)  =  Eoe'^  ^*/(l  + 


(123) 


where  Eq  =  constant.  The  spot  size  at  distance  z  is  given  by 
w  =  a(l  +  ^2^1/2 


where  ^  is  the  normalized  longitudinal  distance  from  the  beam  focus.  At 
the  beam  focus  the  spot  size  is  minimum  and  is  denoted  by  a.  When  the  non¬ 
linear  crystal  is  placed  in  the  near  field  ^  »  1  and  the  beam  is  nearly  paral¬ 
lel.  The  pump  electric  field  distribution  just  inside  the  entrance  surface  of 
the  material,  for  the  case  considered,  is  given  by 


Ep(p,  -  f  )  =  E„ 

The  angular  spectrum  of  pump  field  is 

R  D  ^  /4) 

D 


(124) 


(125) 


The  angular  spectrum  z  -  -  -  of  object  field  emanating 

Ct 


from  z  =  -  s  is  again 
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Then  the  angular  spectrum  of  polarization  wave  inside  the  crystal  becomes 


(P  (Kg.  7.)  =  a^‘^pdTi  eil(K 

a2  2 

•/  ‘5.(Ks- Kp,-s)e  -4  '^p  e 


-a%p/4 

e 


(126) 


The  angular  spectrum  of  the  upconverted  image  field  is  obtained  following 
Eq.  (69)  and  (97) 
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For  a  point  object  located  at  the  point  (^j,  -  s),  the  angular  spectrum  of  the 
image  field  is 
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where  integral  I  is  given  by 
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When  pj  =  0  and  a2  » _ ^  which  corresponds  to  the  Fresnel  zone  of 


la 


Gaussian  aperture,  the  expression  for  the  angular  spectrum  of  upconverted 
image  field  becomes. 
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where 

z  D  D  D  D  s 
♦  ("s)"-!®  Itsa  2Ksa’"  2K^  ~  2Kia  *Kia 

,-ax^D 
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Above  equation  is  identical  to  Eq,  (100)  obtained  for  plane  wave  pump  case. 
The  discussion  on  aberration  for  plane  wave  pump  system,  therefore, 
directly  apply  to  the  Gaussian  wave  pump  system.  The  effect  of  Gaussian 
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0^22  2 

field  distribution  appears  as  an  attenuation  factor  exp[  -  (s-^  )  k  /a  ] 


It  can  be  shown  that  Eq.  (120)  becomes  equal  to  Eq.  (98), 


if  the  Gaussian  attenuation  factor  exp 


D 

®‘2  2 


la 


is  neglected  in  Eq.  (120) 


and  ^^(^s ,  -  s)  in  Eq.  (08)  Is  set  equal  to  ,  which  is  Fourier  component  of 
point  object  located  at  =  0  and  z  =  -  s.  This  indicates  that  the  Gaussian 
beam  pumped  upconverter  can  be  considered  to  be  equivalent  to  the  plane 


D  X  D  X 

wave  pumped  system  if  (s- «)-^«a.  The  condition  (s-—)-^  «a  implies 

2  Kia  2  Kia 

that  the  beam  size  of  Gaussian  pump  beam  is  considerably  greater  than  the 


interaction  aperture  determined  by  the  object  field  of  view. 


As  shown  for  the  plane  pump  wave  case  the  resolving  power 
of  the  Gaussian  wave  pumped  system  remains  constant  for  the  range 


D  K 

a  >(s--5-)-^  .  The  resolving  power  of  the  Gaussian  wave  pumped  system 
^  '^ia 

D  K 

would  decrease  with  the  beam  size  for  the  range  a<  (s-t-)-^^  due  to  the 

^  *Sla 

Gaussian  field  distribution. 

ni-4  FINITE  INTERACTION  APERTURE 

Finite  interaction  aperature  resulting  from  the  finite  trans¬ 
verse  dimension  of  the  material  can  be  accounted  for  by  introducing  aperture 
function 
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pi' 

-’2  [0  p  outside  A 


(130) 


If  we  express  pump  field  at  z  =  by  E  (p,  then 

^  \r  ^ 


effective  pump  field  distribution  is 


_eff  „  ,  D  \  n/  Dv 


(131) 


The  angular  spectrum  of  effective  pump  field  is, 


(132) 


D 


where  (B(k  ,  -  «  )  is  angular  spectrum  of  aperture  function  Eq.  (130). 

r  2 

Similarly  the  angular  spectrum  of  effective  object  field  is  given  by 


^  ‘-i’  ■  f '®  ®  '-1  ■? '  ' 

The  introduction  of  a  finite  interaction  aperture  due  to  the 

material  dimension  or  other  aperture  stops  has  the  effect  of  broadening 

the  angular  spectrum  of  both  pump  and  object  field  disturbance.  The  smaller 

the  aperture,  the  broader  the  angular  spectrum.  Since  angular  spectrum  of 

the  sum-frequency  image  field  is  proportional  to  convolution  of  Eq.  (132)  and 

(133),  image  angular  spectrum  will  be  a  function  of  convolution  of  aperture 

spectrum  ;(B  0  (B.  This  indicates  that  the  effective  aperture  for  the  sum- 

frequency  field  is  smaller  than  the  material  aperture  A.  For  the  case  of 

uniform  illuminated  plane-wave  interaction,  the  effective  interaction  area 

for  the  sum-frequency  field  is  one  half  of  the  material  aperture.  This  effect 

(29) 

has  also  been  pointed  out  by  Boyd  and  Ashkin. 
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«f^K,-f)  =  «i(£l,-f  )  (133) 

The  introduction  of  a  finite  interaction  aperture  due  to  the 

material  dimension  or  other  aperture  stops  has  the  effect  of  broadening 

the  angular  spectrum  of  both  pump  and  object  field  disturbance.  The  smaller 

the  aperture,  the  broader  the  angular  spectrum.  Since  angular  spectrum  of 

the  sum-frequency  image  field  is  proportional  to  convolution  of  Eq.  (132)  and 

(133),  image  angular  spectrum  will  be  a  function  of  convolution  of  aperture 

spectrum  ;<B  @  <B .  This  indicates  that  the  effective  aperture  for  the  sum- 

frequency  field  is  smaller  than  the  material  aperture  A.  For  the  case  of 

uniform  illuminated  plane-wave  interaction,  the  effective  interaction  area 

for  the  sum-frequency  field  is  one  half  of  the  material  aperture.  This  effect 

(29) 

has  also  been  pointed  out  by  Boyd  and  Ashkin. 
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CHAPTER  4  SENSITIVITY  OF  IMAGE  UPCONVERSION  SYSTEM 
I.  INT^DUCTION 

An  important  consideration  in  defining  the  performance  of  any 

image  system  is  the  sensitivity.  One  way  of  describing  the  sensitivity  is 

noise-equivalent-power  (NEP)  per  single  picture  element.  This  is  a  useful 

criterion  of  some  devices  and  it  is  usually  a  true  measure  of  sensitivity  that 

does  not  depend  on  the  optics  of  system,  exposure  time,  and  amount  of  de- 

(31) 

tail  in  the  reproduced  image  .  Another  way  of  comparing  sensitivities  is 

(31) 

in  terms  of  the  threshold  scene  radiance  .  Threshold  scene  radiance  de¬ 
pends  explicitly  upon  the  exposure  time,  system  optics,  and  the  image  de¬ 
tail. 

The  sensitivity,  as  expressed  by  NEP,  of  photoelectronic  image 
tubes  and  photographic  films  considered  by  themselves  is  only  a  function  of 
quantum  efficiency  of  photosensitive  surface  and/or  noise  due  to  dark  leakage 
and  the  amplifier  which  follows  the  photo -sensitive  device.  NEP  of  the  image 
system  utilizing  parametric  upconversion  process  is,  however,  implicitly 
dependent  upon  amount  of  image  detail  as  is  shown  by  the  following  considera¬ 
tion.  For  a  given  pump  power,  set  either  by  available  laser  power  output  or 
by  the  power  handling  capacity  of  the  nonlinear  material,  attainable  conver¬ 
sion  efficiency  of  upconversion  process  depends  on  the  interaction  aperture. 
Number  of  resolution  elements  of  the  diffraction  limited  image  upconverter 
increases  with  the  interaction  area.  Conversion  efficiency  of  upconversion 
is  inversely  proportional  to  the  interaction  area.  It  then  follows  that  the  con¬ 
version  efficiency  decreases  as  more  detail  is  required  in  the  reproduced 
image.  The  reduction  in  conversion  efficiency  directly  leads  to  increase  in 
NEP.  Sensitivity  of  the  image  upconverter  expressed  in  terms  of  NEP  is. 
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therefore,  no  more  true  than  the  sensitivity  description  in  terms  of  threshold 
scene  radiance.  The  NEP  criterion  is  for  a  single  picture  element  and  is 
defined  for  unity  signal-to-noise  ratio,  therefore  it  can  not  be  used  to  dis¬ 
cus  5  the  system  that  must  see  half  tone  pictures.  The  NEP  criterion  is 
normalized  for  unit  bandwidth  or  unit  observation  time  and  hence  it  is  in¬ 
dependent  of  exposure  time.  The  above  considerations  suggests  the  desir¬ 
ability  of  expressing  the  system  sensitivity  in  a  more  intrinsic  fashion  so 
that  the  threshold  scene  radiance,  image  detail  including  resolution  and  con¬ 
trast,  conversion  efficiency,  and  the  interaction  aperture  would  be  related 
to  some  universal  constant.  It  is  tie  purpose  of  this  chapter  to  obtain  a 
characterization  of  the  image  upconversion  system  which  will  involve  both 
sensitivity  and  picture  detail.  Specific  attention  is  paid  to  the  system  using 
a  combination  of  an  upconverter  and  a  photoelectro nic  image  recording  device. 

An  equation  characterizing  of  an  ideal  image  upconversion  sys¬ 
tem  consisting  of  a  quantum  noise  limited  upconverter  and  a  perfect  picture 
pick-up  device  is  derived  in  Section  11.  The  threshold  quantum  number  of  the 
photoelectron  camera  tubes  is  considered  in  Section  III.  Threshold  quantum 
number  of  picture  tubes  is  defined  as  minimum  required  photo  quanta  per 
picture  element  per  integration  time  to  produce  signal-to-noise  ratio  of 
imity  at  the  output.  The  characteristic  equation  of  a  practical  upconversion 
system  using  an  imperfect  pick-up  device  is  also  considered  in  Section  in. 
Finally  the  NEP  expression  is  presented  in  Section  IV. 

II.  CHARACTERISTIC  EQUATION  OF  AN  IDEAL  IMAGE  UPCONVERTER 
SYSTEM  (QUANTUM  FLUCTUATION  LIMITED  PERFORMANCE) 

Consider  a  square  element  of  area  of  side  length  "a"  on  the  up- 
converted  image.  It  is  assumed  that  a  picture  pickup  device,  with  an  ideal 
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photosensitive  surface,  isusedto  record  the  upconverted  image.  By  "ideal” 

is  meant  that  the  photosensitive  surface  will  record  single  incident  quantum. 

2 

Let  this  element  of  area  a  absorb  on  the  average  Q  quanta  of  wavelength 

s 

X  in  the  exposure  time  t  alloted.  Because  at  a  room  temperature  the  absorp- 
s  (32  33) 

tion  of  near  infrared  and  visible  photo  quanta  obeys  Poission  statistics  ’  , 

the  average  absorption  number  Q  will  have  associated  with  it  deviations 

s  1 

from  the  average  whose  root  mean  square  value  is  Qg^-  This  deviation  set 

a  limit  to  the  accuracy  with  which  the  average  number  Q  may  be  determined. 

s 

The  smallest  change  in  Q  which  can  be  detected  is,  therefore,  of  the  order 
1  s 

ofQ  ^ 
s 

Next  consider  two  adjacent  square  elements  A  and  B.  The  dif¬ 
ference  in  average  number  of  quanta  arriving  at  tliese  two  elements  can  be 
distinguished  provided: 


1/2 

(Qo .  -  *  I  .  .  =  [  (QgA  * 


This  relation,  indicates  that  the  smallest  intensity  difference  between  adjacent 

two  elements  (denoted  by  AQ  )  which  can  be  detected  is  also  of  order  of  the 

1/2® 

statistical  fluctuations  Q  ,  where  Q  is  an  average  quantum  flux. 

®  ®  2 

The  irradiance  H  ,  in  watt/cm  ,  is  related  to  Q  in  the  following 
s  s 


way 


hi/  Q 

TT  ®  ® 
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(134) 


We  can  now  introduce  the  threshold  contrast  C  as 

th 
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For  most  of  any  application  H 


C 


th 


SA 


H  H  ,  therefore 
SB  S 


-1/2 


(135) 


From  (134)  and  (135) 


or 

Of^  =  const.  (136) 

S  th 

? 

where  a  is  the  angle  subtended  by  the  element  a  at  the  nonlinear  crystal. 
(Fig.  14) 

Radiance  H^,  of  the  infrared  object  which  is  upconverted  is  pro¬ 
portional  to  the  image  irradiance  H  ,  since  the  input  electric  field  (power 

O 

density)  and  the  output  electric  field  (power  density)  are  linearly  related  and 

(37) 

the  upconversion  process  is  noiseless  .  Eq.  (136),  therefore,  can  be  re¬ 
written  as: 

H.C..  =  const  (137) 

1  th 


Eq.  (137)  is  the  characteristic  equation  of  an  ideal  image  upconverter  system. 
The  constant  term  on  the  right  depends  on  the  particular  system  under  con¬ 
sideration.  The  meaning  of  Eq.  (137)  is  this:  let  any  two  of  the  quantities 

(H,,  C.. ,  u()  be  specified,  then  Eq.  (137)  sets  the  threshold  value  for  the  third, 
i  th 

To  evaluate  the  constant  term  on  the  right  of  Eq.  (137),  we  con¬ 
sider  the  upconverter  system  shown  in  Figure  14. 
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1) 


To  simplify  the  analysis  the  following  assumptions  were  made: 

The  thickness  of  the  nonlinear  material  is  very  small  compared 
to  other  parameters,  such  as  object  and  image  locations,  so 
that  aberrations  due  to  finite  thickness  are  negligible. 

2)  Phase  mismatch  factor  for  all  interacting  rays  is  negligibly 
small  so  that  the  conversion  efficiency  is  constant  over  whole 
angular  aperture.  This  assumption  implies  that  the  angle  sub- 
tanded  by  the  infrared  object  at  the  converter  is  very  small 
compared  to  the  field  of  view  of  the  system.  It  also  assures 
that  all  rays  considered  are  para.xinl  and  small  angle  approxima¬ 
tion  can  be  used. 

If  the  object  radiance  is  such  that  quanta  at  the  wavelength 

are  radiated  per  unit  area  per  unit  time,  the  total  number  of  quanta  radiated 

by  the  infrared  element  with  the  side  length  "b"  for  exposure  time  t  will  be 
2 

N^b  t.  Of  this  number  only  a  fraction  will  be  intercepted  by  the  nonlinear 

crystal  of  active  aperture  diameter  A  situated  distance  &  from  the  infrared 

object.  The  infrared  quanta  intercepted  by  the  nonlinear  crystal  in  time  t 

2  2 

is,  lor  a  Lambart  distribution  equal  to  N^b  t  Sin  0.  where  20  is  the  angle 

substended  by  the  nonlinear  crystal  at  the  object.  These  amount  of  quanta 

will  interact  with  the  pump  quanta  and  together  generate  quanta  at  wavelength 

A  .  The  number  of  quanta  Q  generated  at  wavelength  A  is  given  by 
s  s  s 

Q  =  nN,b^t  sin^  0  (138) 

s  i 

Where  r;  is  photon  (quantum)  conversion  efficiency  of  the  upconversion  process 
and  is  assumed  to  be  a  constant  for  a  given  system. 
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All  generated  quanta  Q  will  be  directed  and  focused  on  the  cle- 
2  ** 

ment  of  area  a  at  distance  S'  from  the  nonlinear  crystal.  From  the  u'com- 
etry  of  Figure  14  and  small  angle  approximation, 


Sin  ~  tan  0  =  — 

Transverse  magnification  factor  a/b  is  obtained  in  Chapter  2  as 

.  Pn  S' 
a  ^  s 

b  ■  n.S 
1 


(139) 


(140) 


where  n  and  n.  are  indices  of  refraction  of  nonlinear  material  at  A  and  a. 
si  si 

respectively. 

/3  is  the  angular  demagnification  factor  defined  in  Chapter  2. 
Substituting  Eq.  (139)  and  Eq.  (140)  in  Eq.  (138),  become 


/n.A  \  2 

\  s  / 


(141) 


a  2 

where  o  =  -,  is  the  angle  substended  by  the  element  of  area  a  at  the  non- 

d 

linear  crystal. 

Using  the  relation  we  obtain  from  Eq.  (141) 

hp  / 2pn 

ww  i  I  s  1  ^ 


(142) 


Combining  Eq.  (134)  and  Eq.  (142),  we  obtain  the  characteristic 


equation; 


4hv.  /  ^  V 

*  th  .  .2  \  n,  / 

tr/A  \  1  / 


(143) 
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Threshold  contrast  defined  in  Eq.  (134)  is  over  the  upcon- 
verted  ima^^c  field,  but  it  also  represents  threshold  contrast  of  infrared 
object  IxHause  of  the  linearity  of  the  upconversion  process.  The  angular 
size  O'  in  the  left  hand  of  Eq.  (143)  is  for  the  upconverted  image.  From 
Eq.  (140)  the  angular  size  )  of  the  infrared  object  element  is  related  to  a 
by 

a  n 

s 


cv  -  r 


n. 

1 


The  final  form  of  the  characteristic  equation  which  describes  the 
ideal  performance  of  image  upconversion  system  is  then 

hr'. 

(144) 


2  2 

H  J.,  V 
1  th 


i 

A\2 


or 


The  si.\  parameters  needed  to  specify  an  ideal  system  are  in  left 
hand  of  Eq.  (144).  When  any  five  of  these  are  arbitrarily  specified,  the 
threshold  value  for  the  sixth  is  given  by  Eq.  (144).  It  should  be  made  clear 
at  this  point  that  the  parameters  > ,  tj,  and  A  can  not  be  arbitrarily  specified 
independently  from  one  another  as  will  be  shown  later. 

Ill .  PERFORMANCE  OF  A  NW -IDE AL  IM AG E  U PCON VE RTE R 

SYSTEM 

The  characteristic  equation  (114)  describing  the  performance  of 
an  ideal  image  upconverter  system  was  obtained  by  assuming: 

1)  Upconversion  process  is  noise less^^*^^ 
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2) 


All  the  radiation  intercepted  by  the  photosensitive  surface  of 
the  pickup  device  arises  from  the  radiation  coming’  from  the 
object.  -  i.e.  no  background  noise 
3)  A  noise -free  image  detector 

The  second  assumjjtion  of  background  noise -free  ((jndition  is 
usually  met  in  practice  since  the  background  thermal  radiation  at  the  sum- 
frequency  and  signal  frequency  wavelength  is  small  enough  to  be  neglected 
at  room  temperature.  The  pump  radiation,  whose  wavelength  is  very  close 
to  the  sum -frequency  wavelength,  can  practically  be  eliminated  by  inter¬ 
ference  filters.  The  use  of  optical  filters  in  the  system  attenuates  the  sum- 
frequency  radiation.  In  the  following  discussion  this  loss  will  iDe  treated  by 
considering  the  filters  as  an  integral  part  of  the  image  detector.  A  noise - 
free  image  detector  is  not  achievable  at  present,  i.e.  the  photosensitive  sur¬ 
face  requires  more  than  one  quanta  to  record  incoming  information.  An 
equation  describing  the  performance  of  a  practical  image  upconverter  sys¬ 
tem  must  therefore  include  the  properties  ol  the  image  detector. 

Any  image  detecting  device,  such  as  a  photoelectronic  camera 
tubes,  a  photographic  film  and  even  the  human  eye  can  be  used  in  a 

practical  image  upconverter  system.  In  the  following,  we  shall  consider 

specifically  the  sensitivity  of  photoelectronic  picture  reproducing  devices. 

If  the  target  of  a  photoelectronic  pick-up  system  is  read  off  by 

means  of  a  scarming  electron  beam  with  frame  time  t  ,  and  if  the  number 

o 

of  picture  elements  is  M,  the  sampling  time  t  of  each  picture  element  is 

We  will  assume  that  the  highest  frequency  in  the  image  signal  will 
be  roughly  F  =  1/2  t,  where  t  is  the  picture  element  samping  (dwell) 
time.  If  the  average  number  of  photo-quanta  absorbed  in  one  picture  element 
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of  area  a  between  storage  time  intervals  t  is  Q  and  the  quantum  effi- 

s 

ciency  of  the  target  is  a,  an  electron  charge  oQ  e  is  stored  in  a  picture 

s 

element  between  two  scans*.  Since  this  charge  will  be  read  off  in  time  t, 

the  average  signal  current  i  for  one  picture  element 

s 


aQ  e 
s 


(145) 


(33,  34) 

The  shot  noise  associared  with  this  signal  current  is 


(^i  )  ^  =  2  ei  F  (146) 

s  s 

If  it  is  necessary  to  amplify  the  signal,  the  noise  current  of 
the  first  amplifier,  (  and  thermal  noise  of  the  amplifier  input  resis¬ 

tance  (Ai^)^  have  to  be  added  to  the  shot  noise.  In  general,  the  dark  cur- 

- p 

rent  noise,  (ni^)^,  generated  either  by  thermal  origin,  charge  leakage,  or 
scanning  beam  current  exists  in  any  pick-up  devices,  and  must  be  added  to 
the  above  noise.  The  total  mean  squared  noise  current  of  the  system  is, 
therefore,  given  by 


(ui  )2  =  b(ai  )^  +  (Ai  )2  +  (Ai  )^  +  (Ai  )^  (147) 

n  s  d  R  A 

The  constant  b  is  introduced  in  above  expression  to  take  into  account  of 

(34) 

noise  figure  of  the  pickup  tubes  .  The  signal -to -noise  ratio  is  then  given 
by 


^  _ 

N  ■  V 


(148) 


*  The  storage  time  interval  t  is  not  in  general  equal  to  the  frame  time 

interval  t  . 
o 
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The  threshold  quantum  number  q  (or  noise  equivalent  quantum  input)  per 

s 

picture  element  per  exposure  time  is  defined  as  that  signal  which  gives  rise 

to  a  unit  signal -to -noise  ratio. 

We  shall  next  evaluate  the  threshold  quantum  number  for 

several  different  classes  of  photoelectronic  camera  tubes. 

1.  Storage  type  with  electron  multiplication 

The  iconoscope,  the  image  orthicon,  the  vidicon-multiplier  belong 
(35) 

to  this  class  .  The  electron  multiplication  is  practically  noise¬ 
less  so  that  the  amplifier  noise  and  the  noise  of  signal  resistance 

(34,  36) 

are  negligible  compared  to  other  noise  .  The  signal-to- 

noise  ratio  for  this  case  becomes 


If  the  device  is  quantum  noise  limited,  that  is  to  say  (Ai  )2  » 

- 2  ® 

(Ai^)  ,  the  signal -to -noise  ratio  is,  using  (145)  and  (146) 

I  =  .  (149) 

N  S 


The  threshold  quantum  number  q  is  given  by 

s 

q  =  1/(7 
s 

2  2 

If  the  device  is  dark  iioise  limited,  (Ai ,)  »  (Ai  )  , 

Q  S 

quantum  number  is  given  by: 

^s  eu  '^(Ai  )2 
d 


(150) 

the  threshold 

(151) 
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2. 


Storage  type  with  an  amplitier 

The  vidicon  and  the  orthicon  belong  to  this  class.  In  these  devices 

(35) 

there  is  no  direct  photoelectronic  multiplication  .  The  amplifier 
noise  and  the  signal  resistance  noise  are  now  dominant  noise 
sources  and  the  threshold  quantum  number  is  given  by: 


s  ea 


(152) 


3.  Nonstorage  Type  with  Multiplier 

(35) 

The  image  disector  tube  belongs  to  this  category  .  This  type 
is  less  sensitive  than  the  storage  type  with  electron  multiplication, 
because  the  photoelectric  charge  is  generated  during  only  a  pic¬ 
ture  element  sampling  time  t  instead  of  a  storage  time  t. 

Let  Q  '  denotes  the  average  number  photo -quanta  absorbed  per 
s 

picture  element  in  the  sampling  (read  out)  time  t.  The  average  signal  cur¬ 


rent  per  picture -element  is 


oeQ  ' 
s 


Fluctuation  noise  current  associated  with  this  current  is  again 


(Ai  )^  =  Aei  F-/ 

S  V  s  ^  T 


therefore  the  signal  to  noise  ratio  for  quantum- limited  case,  is 


'2ei  F 
s 


If  incoming  photon  flux  remains  the  same  as  the  previous  devices  discussed, 

Q  '  Is  related  to  Q  by  the  expression 
s  s 
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Q  '  =  -  Q 
t  ^  s 


so  that  the  signal-to-noise  ratio  becomes 

S/N 


v/  IV 


M 


(153) 


(154) 


where 


From  the  above  we  see  that  the  threshold  quantum  number,  q  ,  is  -  .  There - 

s  o 

fore  the  threshold  quantum  number  is  greater  than  that  encountered  in  storage 
devices  by  a  factor  of  M. 

As  shown  above,  threshold  quantum  number  per  picture  element 

per  storage  time  of  a  photoelectronic  camera  tube  follows  two  different  type 

of  relations  depending  upon  whether  it  is  operating  in  quantum  limited  region 

or  amplifier  (or  dark)  noise  limited  region. 

The  characteristic  equation  of  a  practical  image  upconvertor 

system  can  now  be  derived.  The  number  of  quanta  at  sum-frequency  per 

picture  element  required  to  produce  unity  signal-to-noise  ratio  is  q  and 

s 

the  corresponding  minimum  number  of  quanta  at  the  signal  frequency  is  then 

q  /rj.  The  characteristic  equation  (144),  for  a  practical  image  upconvertor 
s 

system  can  now  be  written  as 


H.c,, 

1  th 


2  2 

r  t 


(1)^ 


V 


hu. 


(155) 
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IV. 


NOISE  EQUIVALENT  POWER  OF  AN  IMAGE  UPCONVERTER 


The  noise  equivalent  power  of  a  practical  image  upconverter 
system  is  derived  below  for  comparison  with  the  characteristic  equation  de¬ 
rived  in  preceding  section.  The  noise  equivalent  power  of  a  detection  sys¬ 
tem  is  defined  as  the  signal  input  power  that  gives  rise  to  unit  signal-to- 
noise  ratio  in  unit  bandwidth. 

Threshold  power  at  sum-frequency  per  picture  element  of  an 

image  detector  with  storage  time  t  is  equal 

q  hi; 
s  s 

where  q^  =  threshold  quantum  number  of  the  image  detector  as  defined 
before . 

h  =  Planck’s  constant. 

V  =  Sum-frequency 
s 

threshold  power  at  infrared  signal  frequency  is  then 

q^hi;./T? 

s  1 

where 


T]  -  quantum  conversion  efficiency 
i;^  =  infrared  signal  frequency 


The  above  is  the  threshold  power  for  a  unit  signal -to-noise  when 
the  IF  amplifier  bandwidth  is  F  =  Since  the  noise  in  this  type  of  detector 
is  proportion  to  the  square  root  of  the  IF  bandwidth,  the  threshold  power  for 
unit  bandwidth,  i.  e.  the  NEP,  is 


where 


NEP 

F 


q  hi;. 

NEP 

77VF 

noise  equivalent  power  in  watt/Hz 
1/2t 


1/2 


(156) 
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T  =  sampling  time  per  picture  element  and  is  related  to  the  frame 

scan  time  t  by 
o  ^ 

T  =  to/M 

where 

M  =  is  number  of  picture  element  in  a  picture  frame.  Quantum 

conversion  efficiency  rj  is  inversely  proportional  to  the  interaction 

.  (1) 

aperture  area. 

r?-KVA^  (157) 

where 

A  =  diameter  of  interaction  aperture 

K'  =  constant  for  fixed  available  pump  power 


For  a  diffraction  limited  image  upconversion  system  the  maximum 


number  of  picture  elements  obtainable  is,  using  Rayleigh  criteria,  given  by 

2 


M 


2  0  A 
'  max 

1.22X 


(158) 


where 

2  =  is  angular  aperture  of  the  upconverter. 

Quantum  conversion  efficiency  may  be  expressed  in  terms  of  number  of 
picture  elements. 


(159) 
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Above  expression  is  a  characteristic  equation  of  diffraction  limited  image 
upconversion  process. 

The  constant  K  is  characteristic  of  a  given  upconverter  depending 
only  on  material  constant  and  the  pump  power  and  is  given 


K 


44.8  P  Id  V(b 
_ p'  eff  max 

3 

n.n  n  A  X, 

1  p  s  s  1 


(160) 


The  final  expression  of  NEP  is 

1/2 

hy.q  /Mt  \ 

NEP=:-^(-;^) 

Kt  ^  2  ' 


(161) 


The  storage  time  t  and  the  frame  scanning  time  t^  is  not  neces¬ 
sarily  the  same.  They  are  equal  for  the  storage  type  photo-electronic  image 
device,  but  they  are  not  equal  for  nonstorage  type  pickup  tubes. 

The  NEP  of  the  image  upconversion  system  using  a  storage  type  pickup 
device  is  therefore 


NEP  = 


hi/.q 

IS 

K 


M 

w 


1/2 


(162) 


The  image  upconversion  system  with  a  nonstorage  type  pick-up  device  be¬ 
comes 


NEP  = 


hi/.q 

IS 

K 


M 


3/2 


(2t  ) 
o 


■1/2 


(163) 
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CHAPTER  5  IMAGE  UPCONVERSION  EXPERIMENTS 
I.  GENERAL  CONSIDERATION 

An  upconverter  using  single  crystal  proustite  (Ag  AsS„),  a 
Nd:YAG  laser,  and  an  infrared  10.  6  Mm  illuminated  object  as  its  primary 
components  was  used  for  the  image  upconversion  experiment.  The  Nd:YAG 
laser  radiating  at  1.06  Mm  was  chosen  because  of  its  moderate  CW  power 
and  stable  single  mode  output.  The  objects  were  illuminated  with  the  10.  P 
Mm  output  of  a  CW  COg  gas  laser.  The  1. 06  Mm  pump  radiation  and  infra¬ 
red  10.  6  Mm  object  radiation  were  polarized  as  ordinary  and  extraordinary 
waves,  respectively.  The  resulting  upconverted  image  wave  at  0.  967  Mm 
was  polarized  as  an  extraordinary  wave.  The  proustite  crystal  was  prepared 
in  such  a  way  that  its  two  polished  faces  were  nearly  perpendicular  to  the 
phase- match  (near  collinear)  direction.  The  angular  relations  of  the  crystal 
optic  axis,  the  crystal  surfaces,  and  the  directions  of  beam  propagation 
that  were  used  for  optimum  upconversion  have  been  discussed  in  detail  in 
Section  IV- 2,  Chapter  2. 

The  quality  of  an  imaging  system  is  usually  assessed  either  by 
its  limiting  resolution  or  by  its  transfer  function.  The  classical  limiting 
resolution  criterion  is  convenient  and  easy  to  apply,  but  it  does  not  provide 
the  complete  quantitative  data  required  to  characterize  the  system.  The 
basic  measurement  required  for  the  transfer  function  description  of  a  system 
is  the  spread  function  measurement.  Since  the  spread  function  of  a  system 
is  defined  as  the  system  response  to  a  delta  function  input,  the  Fourier 
transform  of  the  spread  function  is  the  system  transfer  function.  Measure¬ 
ments  of  spread  functions  can  be  performed  either  in  two-dimensions  or 
one-dimension;  that  is,  point  or  line  spread  function.  Tlie  line  spread 
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function  technique  was  chosen  for  the  initial  experiments  because  of  ap¬ 
proximate  symmetry  and  the  ease  of  implementing  the  Fourier  transform 
computation  of  a  line  spread  function  to  obtain  an  optical  transfer  function. 

The  line  spread  functions  were  measured  for  two  types  of  image 
upconverter  configurations.  The  first  configuration  chosen  was  the  Fourier 
space  configuration,  in  which  both  the  pump  and  object  waves  are  nearly 
plane  waves.  The  other  system  utilizes  a  Gaussian  shaped  pump  beam  and 
a  cylindrical  object  wave.  These  measurements  were  made  to  obtain  data 
on  image  location  and  the  amount  of  optical  aberrations.  On-  and  off-axis 
spread  functions  of  the  Fourier  space  configuration  were  measured  to  de¬ 
termine  the  limiting  resolution.  The  primary  consideration  in  choosing  a 
nearly  plane  wave  pump  in  the  experiment  is  the  small  pump  acceptance 
angle  of  the  upconversion  process  using  a  proustite  crystal.  As  shown  in 
Chapter  2,  the  pump  angular  aperture  of  the  1.06  jum  pumped  proustite  is 
less  than  0.  2  degrees.  Effective  usage  of  available  pump  power,  therefore, 
requires  minimum  pump  beam  divergence. 

The  intensity  distribution  was  measured  for  the  upconverted 
image  of  a  resolution  test  bar  chart,  which  was  illuminated  by  10.  6  jum 
radiation.  The  data  was  in  rough  agreement  with  the  predicted  line  spread 
function. 

Measurements  of  the  line  spread  function  and  the  intensity 
distribution  of  the  test  chart  image  were  taken  only  along  the  crystallo¬ 
graphic  a-axis.  This  is  justified  because  of  small  birefringence  of  proustite 
crystal.  Measurements  along  two  orthogonal  transverse  directions,  however, 
becomes  necessary  if  the  crystal  anisotropy  is  appreciable  enough  to  intro¬ 
duce  large  asymmetry. 
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The  upconverted  image  of  a  10.  6  Mm  illuminated  resolution 
test  chart  was  recorded  on  photographic  film.  This  result  provides  a  visual 
display  of  an  upconverted  image  and  agrees  with  the  measurements  men¬ 
tioned  above. 

n.  EXPERIMENTAL  ARRANGEMENT 

The  apparatus  shown  schematically  in  Figure  15  was  arranged 
so  that  the  nature  of  pump  and  IR  object  waves  may  be  readily  changed  with 
the  minimum  physical  alteration  of  components.  A  2.  5  cm  in  diameter  and 
0.  6  cm  thick  proustite  crystal  was  used.  The  crystallographic  a-axis  lies 
in  the  plane  of  paper  and  is  parallel  to  the  front  surface.  Its  optic  axis 
lies  in  a  plane  perpendicular  to  the  paper  and  makes  an  angle  of  20  degrees 
with  the  normal  to  the  polished  front  and  back  surfaces.  A  single  mode 
CW  Nd'.YAG  laser  radiating  an  electric  field  polarized  parallel  to  the  plane 
of  paper  was  used  as  a  pump  source.  Four  glass  lenses,  labeled  1,  2,  3, 
and  4,  were  placed  along  the  pump  beam  path  so  that  the  apparent  location 
of  the  source  and  the  beam  size  can  be  altered  by  adjusting  the  relative 
positions  of  those  lenses.  A  COg  gas  laser  radiating  with  an  electric  field 
polarized  perpendicular  to  the  plane  of  paper  was  used  to  illuminate  the 
transparent  objects.  Six  BaFg  lenses,  No.  5  through  No.  10,  were  used  to 
obtain  the  desired  illumination  and  to  change  the  apparent  object  size  and 
location. 

The  pump  wave  and  the  object  wave  are  combined  by  a  dichroic 
mirror  which  has  a  high  reflectance  coating  for  the  1.  00  Mm  pump  radiation. 
The  insertion  loss  of  the  dichroic  mirror  to  the  IR  10.  6  Mm  radiation  is 
about  3  db.  The  10.  6  Mm  object  field  incident  on  to  the  proustite  crystal  is 
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lens  mirror 
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parametrically  upconverted  to  0.967  pm  image  waves.  The  upconverted 
image  waves  are  then  refocused  by  the  glass  lens  No.  11  to  form  a  real 
image.  The  upconverted  real  image  is  recorded  by  a  image  pick-up  device. 

The  image  pick-up  unit  consisting  of  a  cooled  (60°C)  photo¬ 
multiplier,  a  0.96  Mm  band-pass  filter  and  a  narrow  sampling  slit,  forming 
an  integral  detection  unit,  is  mounted  on  a  translation  platform.  The  up- 
converted  image  is  scanned  by  the  sampling  slit  and  the  outijut  from  the 
photomultiplier  was  amplified  and  fed  to  the  vertical  drive  of  an  X- Y 
recorder.  An  electrical  signal  indicating  the  transverse  position  of  the 
sampling  slit  was  applied  to  the  horizontal  drive  of  the  X- Y  recorder, 
ni.  MEASUREMENT  OF  LINE  SPREAD  FUNCTIONS 

Spherical  lenses  No.  1  and  No.  2  are  used  as  a  pump  beam  ex¬ 
pander.  They  are  positioned  so  that  the  pump  beam  focuses  at  infinity. 
Magnification  of  beam  size  equals  to  the  ratio  of  focal  length  of  input  lens 
No.  1  to  that  of  output  lens  No.  2.  Beam  divergence  is  reduced  also  by  that 
factor.  The  lenses  No.  3  and  No.  4  were  omitted  in  the  line  spread  function 
measurement.  The  spherical  lenses  No.  5  and  No.  6  formed  a  beam  ex¬ 
pander  for  the  CO-  laser  ou^ut  beam.  The  lenses  No.  9  and  No.  10  were 
also  omitted  in  this  experiment.  The  10.  6  Mm  line  source  necessary  for 
the  measurement  was  generated  by  the  lenses  No.  7  and  No.  8  and  a  narrow 
slit  placed  between  these  lenses.  Lenses  No.  7  and  No.  8  are  cylindrical 
and  were  oriented  such  that  their  surface  curvature  lies  in  the  plane  parallel 
to  the  plane  of  paper.  The  narrow  slit  was  placed  at  the  focus  of  lense  No.  7 
and  its  long  edge  was  perpendicular  to  the  paper.  Width  of  the  slit  was 
adjusted  narrower  than  the  beam  waist  of  focused  10.  6  Mm  beam.  The 
position  of  cylindrical  lens  No.  8  was  adjusted  so  that  the  image  of  the 
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10.  6  Mni  illuminated  slit  could  be  positioned  at  desired  location. 
in-1.  Results  on  Fourier  Confijruration  Upconverter 

An  imaj^e  upconverter  system  with  a  plane  pump  wave  and  plane 
object  wave  has  been  shown  to  be  free  of  thickness  aberration.  Its  limiting 
resolution  is  governed  by  higher  order  effects  such  as  birefringence,  dis¬ 
persion.  mismatch,  lens  aberrations,  and  the  aperture  diffraction.  It  is, 
therefore,  essential  to  investigate  the  image  properties  of  this  system  to 
assess  the  optimum  image  quality  attainable  by  the  parametric  image  up¬ 
converter. 

Experimentally  the  plane  wave  pump  condition  is  achieved  by 
placing  the  proustite  crystal  in  the  near  field  of  a  fundamental  mode  laser 
beam  of  large  beam  size.  The  beam  size  of  the  pump  beam  must  be  con¬ 
siderably  larger  than  the  interaction  aperture  determined  by  the  infrared 
field  of  view.  The  10.  6  jum  object  wave  radiating  from  the  object  slit  is 
converted  into  a  plane  wave  by  lens  No.  8.  The  distance  between  the  object 
slit  and  lens  No.  8  is  equal  to  the  focal  lengb-  of  lens  No.  8.  The  10.  6  pm 
radiation  emanating  from  lens  No.  8  is  a  truncated  plane  wave  provided  the 
object  slit  diffracts  the  incoming  10.  6  pm  beam.  The  beam  width  of  the 
collimated  10.  6  pm  radiation  impinging  on  the  proustite  crystal  was  about 
4  mm  in  the  direction  perpendicular  to  the  object  slit.  The  1.  06  pm  pump 
beam  was  approximately  8  mm  in  diameter. 

The  sum-frequency  radiation  emanating  from  the  proustite  crys¬ 
tal  is  also  collimated.  This  was  checked  by  transversely  scanning  the 
0.967  pm  radiation  at  different  distances  from  the  crystal.  A  typical  in¬ 
tensity  profile  of  0.967  pm  radiation  is  shown  in  Figure  16.  Neglecting 
slow  decaying  skirt,  the  beam  width  of  0.  967  pm  radiation  was  approxi¬ 
mately  2.  8  mm.  This  value  agrees  reasonably  well  with  the  expected  beam 
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FIGURE  16.  INTENSITY  PROFILE  OF  SUM-FREQUENCY  BEAM  AT  THE 
IMPUT  SURFACE  OF  IMAGE  LENS  NO.  11. 


size  of  2.  66  mm  calculated  by  the  formula 


L  L  ^  L 

A.  ^  A  A 

1  p  s 

where  A.,  A^,  and  A^  are  respectively  beam  size  of  infrared  object,  pump, 
and  sum-frequency  beams. 

The  line  spread  function,  which  is  a  real  image  of  object  slit, 
is  obtained  by  refocusing  the  0. 967  jum  beam  by  lens  No.  11.  Figure  17 
and  Figure  18  are  the  line  spread  functions  measured  by  scanning  slit  detec¬ 
tion  technique.  A  76.  2  mm  focal  length  lens  was  used  as  lens  No.  11  for 

the  line  spread  function  shown  in  Figure  17.  The  focal  length  of  lens  No.  11 

is  152. 4  mm  for  the  line  spread  function  shown  in  Figure  18.  The  slow 
deca3dng  of  the  right  skirt  of  the  measured  spread  function  is  due  to  the 
long  time  constant  of  the  signal  amplifier  employed.  The  line  width  of  the 
measured  line  spread  functions  is  compared  with  the  diffraction  limited 
line  width.  The  diffraction  limited  line  width  of  uniformly  illuminated 
aperture  is: 


b  =  2  ~  (164) 

A 

where 

b  =  diffraction  limited  line  width 
f  =  focal  length  of  the  lens 
A=  linear  dimension  of  aperture 

For  the  case  considered  here,  A  is  equal  to  the  sum-frequency 

beam  width  A  . 

s 

Table  I  shows  the  measured  and  calculated  line  width  for  the 
two  lenses  used. 
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FIGURE  17.  LINE  SPREAD  FUNCTION  OF  A  FOURIER  CONFIGURATION 
UPCONVERTER  SYOTEM 


FIGURE  18.  UNE  SPREAD  FUNCTION  OF  A  FOURIER  CONFIGURATION 
UPCONVERTER 


Table  I  Theoretical  and  Measured  Line  Width  of 
Line  Spread  Function 


f  mm 

Calculated  by  (164)  Mm 

Measured  Mm 

Percentage 

Error 

76.2 

52.  62 

63.1 

20 

152.4 

105. 24 

129.4 

22 

In  measuring  the  line  width  of  spread  functions  shown  in  Figure 
17  and  Figure  18,  the  slow  decay  was  neglected. 

The  discrepancy  between  measured  and  calculated  line  width 
is  small  enough  to  conclude  that  the  system  tested  is  practically  diffraction 
limited.  The  small  discrepancy  is  believed  to  be  caused  by  the  following 
reasons: 

1)  Diffraction  limited  line  width  of  nonimiformly  illuminated 
aperture  is  wider  than  that  of  uniformly  illuminated 
aperture. 

2)  Finite  width  of  the  sampling  slit  tends  to  make  measured 
line  width  appear  wider  than  the  actual  width. 

III-2.  Results  on  Upconverter  with  Gaussian  Beam  Pump  and  Cylindrical 
Object  Wave 

Experimental  arrangement  and  the  position  of  pertinent  com¬ 
ponents  are  shown  in  Figure  19.  The  1. 06  Mm  pump  beam  was  expanded 
to  8  mm  and  was  collimated  by  lenses  No.  1  and  No.  2.  The  10.  6  Mm 
infrared  beam  was  also  expanded  and  collimated  by  lenses  No.  5  and  No.  6. 
An  125  Mm  wide  object  slit  was  placed  at  the  focus  of  lens  No.  7.  The  lens 
No.  8,  with  92  mm  focal  length,  was  positioned  130  mm  from  the  object 
slit.  The  distance  between  lens  No.  8  and  the  proustite  crystal  was  360  mm. 
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FIGURE  19.  EXPERIMENTAL  ARRANGEMENT  TO  REALIZE  AN  UPCONVERTER 
WITH  GAUSSIAN  BEAM  PUMP  AND  CYLINDRICAL  OBJECT  WAVE 


Lens  No.  11  was  located  219  mm  from  the  crystal. 

The  infrared  object  wave  is  cylindrical,  since  lenses  No.  7 
and  No.  8  are  cylindrical.  The  10.  6  jum  image  of  the  object  slit  was  found 
to  be  at  319  mm  from  lens  No.  8.  The  distance  between  this  infrared  line 
image  and  the  proustite  crystal  was  therefore  41  mm.  The  infrared  image 
of  the  object  slit  can  be  considered  to  be  the  apparent  infrared  object  of 
the  image  upconverter.  Limitation  imposed  by  physical  size  of  components 
such  as  carriages,  holders,  and  mounts  necessitates  usage  of  an  apparent 
object  rather  than  positioning  a  real  object  at  desired  position.  The  infra¬ 
red  object  aperture  at  the  crystal  may  be  calculated  by  the  following  re¬ 
lation: 


A  =  s  d, 

P  1  max. 

Using  measured  value  s  =  41  mm  and  computed  value  9  = 

1  max 

0. 18  rad. ,  it  was  found  that  the  object  aperture  is  about  74  mm.  As 
discussed  in  Chapter  3,  Gaussian  distribution  of  pump  field  cannot  be 
neglected,  since  Gaussian  pump  aperture  and  object  aperture  are  compara¬ 
ble. 


On-axis  line  spread  functions  at  different  distances  from  lens 
No.  11  were  measured  and  are  shown  in  Figure  20.  The  numerical  figures 
on  top  of  the  spread  function  indicate  the  distance  between  lens  No.  11  and 
the  sampling  slit.  The  spread  functions  are  seen  to  remain  reasonably  un¬ 
changed  for  a  large  axial  range.  The  spread  functions  taken  at  17,  16.  73, 
16.4  and  16. 14  cm  appear  similar  except  for  the  difference  in  the  peak 
intensity.  Variation  of  peak  intensity  of  those  spread  functions  is  believed 
to  be  caused  by  long  term  system  instability.  The  line  width  of  those 
spread  functions  is  approximately  400  Mm. 
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FIGURE  20.  LINE  SPREAD  FUNCTIONS  MEASURED  AT  DIFFERENT  DISTANCE 
FROM  LENS  NO.  11  (GAUSSIAN  BEAM  PUMP  AND  CYLINDRICAL 
OBJECT  WAVE  CONFIGURATION) 


To  compare  the  experimental  results  with  the  theory,  the 
axial  position  and  the  line  width  of  the  sum-frequency  image  are  calculated 
by  Eq.  (6)  of  Chapter  2  and  the  lens  formula.  Axial  position  s'  of  the  up- 
converted  virtual  image  is  related  to  the  axial  position  s  of  infrared  object, 
for  p  =  ,  by 


n. 


s’ 


p  =  oo  ^ 


(165) 


It  was  shown  in  Chapter  3  that  the  above  relation  also  applies 

to  the  Gaussian  beam  pumped  upconverter.  Using  parameter  values 

/3  =  0.  0866,  n  =  2.  80  and  n.  =  2.  677,  we  found  s’  I  -11.0s  and  s’  -  -431 
s  1  ' 

mm  for  s  =  -41  mm.  The  axial  position  of  the  real  sum-frequency  image 
formed  by  lens  No.  11  was  calculated  by  the  lens  formula.  The  distance 
between  the  virtual  0.967  pm  image  and  lens  No.  11  is  650  mm,  which  is 
the  sum  of  s’  and  the  distance  between  lens  No.  11  and  the  crystal.  The 
computed  distance  from  lens  No.  11  to  the  real  sum-frequency  image  is 
then  170  mm.  The  measured  mean  distance  was  166  mm.  Agreement  be¬ 
tween  the  experimental  result  and  theoretically  computed  value  is  reason¬ 
ably  well. 

The  width  of  sum-frequency  line  image  was  computed  by  the 

relation: 


W  =  W.  M, ,  M 
s  1  8  11  u 


where  W  and  W.  are  respectively  the  width  of  sum-frequency  image  and 

S  1 

object  slit  width,  M  and  M,,  are  respectively  the  magnification  factor 

of  lenses  No,  8  and  No,  11,  is  transverse  magnification  factor  of  up- 

conversion  process.  Using  known  parameter  values  W.  125  um,  M- 

i  8 

31.9/13,  M,,  =  17/65,  and  M  -  1,  we  found  that  W  is  87  um,  which  is 
IX  \x  s 
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more  than  four  times  less  than  the  measured  value  of  400  um.  It  must  be 
noted  that  the  above  computation  was  made  without  taking  into  account  of 
thickness  aberrations.  The  amount  of  thickness  aberration  in  the  plane  wave 
pumped  upconverter  will  be  considered  next. 

It  was  shown  that  thickness  aberration  is  independent  of  object 
location  if  the  pump  wave  is  plane  wave  and  the  object  is  at  finite  distance 
(Chapter  2).  Longitudinal  aberration  I  As*  I  for  this  case  is  given  by  Eq. 
(30),  Chapter  2.  For  a  1  cm  long  proustite  crystal  used  in  the  experiment, 

I  A  s'  I  was  found  to  be  376  mm.  Longitudinal  spread  I  A  s"l  of  the  real 
sum-frequency  image  formed  by  lens  No.  11  is  related  to  I  A  s*  I  by  the 
relation 

1  A  s”  I  -  I  As’  I 


It  was  found  that  I  As"  I  is  equal  to  2.  57  mm.  The  minimum  resolvable  line 
width  may  be  calculated  by  the  relation 


(W  )  .  =  -  I  As'l  (0  ) 
s  min  2  s  max 


Using  I  As*  I  =  3.  76  cm,  {0  )  =  {9,)  A  /x.  =  0.016  rad.,  and  M,,  = 

s  max  1  max  s  i  11 

17/65,  we  obtain  that  approximately  85  Mm.  Neglecting  the 

effect  of  thickness  coma,  the  line  width  of  plane  wave  pumped  upconverter 

should  be  172  um,  which  equals  the  sum  of  aberration  free  image  width  and 

the  minimum  resolvable  width  (W.)  .  .  The  results  of  the  above  computation 

i  mini 

and  measured  results  are  tabulated  in  Table  II. 
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Table  11  Comparison  of  Measured  Results  and 
Theoretically  Computed  Results 


Computed 


Aberration 

Free 

Plane  Wave 

Pumped  Sys. 

Measured 

Measured  Value/ 
Plane  wave 
Pumped  Value 

Line  width,  pm 

87 

172 

400 

2.  33 

Longitudinal 
spread  of 

0.  967  um 
image,  mm 

0 

2.  57 

i 

8.  6 

2.  21 

Additional  width  broadening  and  longitudinal  spread  over  the 
plane  wave  pumped  case  is  due  to  Gaussian  distribution  of  pump  field  as 
discussed  in  Chapter  3. 

The  above  computation  of  aberrations  of  the  plane  wave  pumped 

system  was  made  with  the  results  obtained  in  the  geometric  optics  theory. 

The  minimum  resolvable  spot  size  of  sum-frequency  virtual  image  can  also 

be  calculated  by  Eq.  (107),  Chapter  3.  Using  parameter  value  D  =  1  cm, 

4  -1 

K.  =  2  7r  n.  /A,  =  1.  58  x  10  cm  ,  we  obtain  the  resolvable  spot  size  2  p 
111  ^  '^res 

equals  to  223  um.  The  geometric  optics  theory  gives  the  resolvable  spot 

size  of  virtual  sum-frequency  image  as  1/2  I  A  s’  I  0  ,  and  it  is  found 

s  max 

to  be  293  um.  Agreement  between  these  values  is  reasonable. 

IV .  MEASUREMENTS  OF  RESOLUTION  AND  INFRARED  OBJECT  FIELD 
OF  VIEW 

The  number  of  linear  resolution  elements  was  measured  with 
ihe  Fourier  configuration  upconverter.  The  Fourier  configuration  was 
chosen  because  it  is  nearly  diffraction  limited.  Measurement  described 
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below  allows  us  to  obtain  information  on  resolution  and  field  of  view.  The 

experimental  arrangement  is  shown  schematically  in  Figure  21.  Lenses 

No.  1  and  No.  2  were  again  used  as  a  pump  beam  expander.  Lenses  No.  5 

and  No.  6  formed  a  beam  expander  for  10.  6  um  beam.  Lenses  No.  7  and 

No.  8  and  a  narrow  slit  were  used  to  generate  a  collimated  10.  6  pm  ribbon 

beam.  Mirrors  No.  1  and  No.  2  were  adjusted  so  that  the  1. 06  um  beam 

travels  along  the  direction  which  makes  desired  angle  with  the  z-axis.  The 

infrared  10.  6  um  beam  was  guided  by  mirror  No.  3  to  propagate  in  the 

direction  of  the  z-axis.  Mirror  No.  2  was  arranged  so  that  infrared  10.  6 

um  beam  misses  it.  Mirror  No.  3  was  mounted  on  a  sliding  carriage  so 

that  it  can  be  moved  along  the  direction  parallel  to  the  y-axis.  The  vertical 

angle  a  shown  in  Figure  21  was  adjusted  for  the  maximum  object  angular 

aperture  as  discussed  in  Chapter  2. 

Pump  beam  width  A  was  made  considerably  larger  than  object 

P 

beam  width  A.  so  that  the  effect  of  Gaussian  distribution  is  minimized. 

1 

Angle  0^  that  the  wave  front  normal  of  infrared  object  wave 
makes  with  the  z-axis  is  related  to  the  transverse  position  of  mirror  No.  3 
by  the  relation 

0.  =  tan  ^  (y./f) 

where  C.  is  the  distance  from  the  proustite  crystal  to  mirror  No.  3. 

The  wave  normal  of  sum-frequency  beam  emanating  from  the 

crystal  makes  an  angle  0  with  the  z-axis. 

n 

9  /3  0.  (166) 

s  n.  1 
1 

The  sum-frequency  image  formed  by  lens  No.  11  will  be  dis¬ 
placed  a  distance  y  from  the  z-axis 

s 
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FIGURE  21.  EXPERIMENTAL  SETUP  FOR  MEASUREMENTS  OF  RESOLUTION 
AND  FIELD  OF  VIEW 


(5-4) 


where  is  the  focal  length  of  lens  No.  11 

The  measured  on- axis  and  off-axis  line  spread  functions  are 
shown  in  Figure  22.  The  numerical  value  on  top  of  each  spread  function 
indicates  transverse  position  of  mirror  No.  3.  Measured  and  theoreti¬ 
cally  calculated  transverse  location  y  of  sum-frequency  image  are  tabu- 

s 

lated  in  Table  III. 


Table  III  ^asure^a^  Computed  I'ransverse  Location  of 
0.967  pm  image  (jf  =160  cm  and  -  76.  2  mm) 


y.  cm 

Computed 

Ys 

Measured 

ys  Mm 

6^  Computed 

Degree 

Radian 

20.  32 

893 

862 

+  7.3 

0.129 

15.24 

660 

616 

+  5.47 

0.0953 

4.45 

193 

192 

+  1.59 

0.0278 

0 

0 

0 

0 

0 

-3.05 

-137 

-134 

-1.  09 

0.0191 

-13.97 

-606 

-585 

-5.02 

0.0874 

The  width  of  line  spread  functions  shown  in  Figure  22  appears 
wider  than  actual  width  because  of  the  fast  scan  speed  and  the  long  time 
constant  of  the  amplifier.  Fast  scan  speed  was  required  in  order  to  make 
overall  measurement  time  interval  short  to  avoid  the  effect  of  system  instability 
and  electrical  interference.  Separate  measurements  of  an  off-axis  line 
spread  function  shown  in  Figure  23  indicates  that  the  line  width  is  approxi¬ 
mately  70  pm,  which  is  slightly  larger  than  on-axis  line  width  shown  in 
Figure  17. 
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FIGURE  22.  ON- AXIS  AND  OFF-AXIS  SPREAD  FUNCTION  OF 
FOURIER  CONFIGURATION  UPCONVERTER 


TRED  OFF-AXIS  LINE  SPREAD  FUNCTION 
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The  nuinbfi’  of  linear  resolution  elements  in  (he  10.  6  /Jin  an[j;ular 
aperture  now  enn  be  calculated  with  measured  line  width  and  the'  amount  of 
transverse  dimension  of  sum-frequency  image  field  sliown  in  Figure  22. 
Using  the  Rayleigh  criterion,  we  found  that  the  number  of  linear  resolution 
elements  is  41  lines  over  tne  ol)ject  angular  aperture  of  12.  3  degrees. 

The  number  of  linear  resolution  elements  in  the  angular  aperture  of  dif¬ 
fraction  limited  upconverter  can  be  theoretically  estimated 

0.  7?  A  A 

1  max  s  s 

^  "  1.22  X  n. 

s  1 

where  N  is  the  number  of  linear  resolution  elements.  The  above  relation 
is  obtained  by  combining  (164),  (166),  and  (167).  Assuming  A,  =  2.8  mm 
and  using  known  parameter  values,  we  found  that  N  is  39  lines,  which  is  a 
good  agreement  with  the  measured  result. 

V.  SUM- FREQUENCY  IMAGES  OF  RESOLUTION  TEST  CHART  WITH 
FOURIER  CONFIGURATION  UPCONVERTER 

Exijerimental  investigation  of  0.967  pm  image  of  an  10.6  pm 
illuminated  resolution  test  chart  is  presented  in  this  section.  The  experi¬ 
mental  arrangement  shown  schematically  in  Figure  15  was  used.  Lenses  No.  1 
and  No.  2  constitute  the  1.06  urn  pump  beam  expander.  Lenses  No.  3  and 
No.  4  are  cylindrical  and  were  used  as  a  one  dimensional  beam  compressor. 
The  pump  beam  emanating  from  lens  No.  4  is  collimated  and  has  nearly 
rectangular  cross-section.  Tlie  long  side  of  the  beam  cross-section  is 
about  8  mm  and  is  in  the  direction  parallel  to  the  plane  of  paper.  The  short 
side  of  the  beam  cross-section  is  about  3  mm.  Lenses  No.  5  and  No.  6 
constitute  the  10.  6  um  beam  expander.  Lenses  No.  7  and  No.  8  are  cylin¬ 
drical  and  were  used  as  one  dimensional  beam  compressor.  The  10.  6  um 
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beam  onuiiiatini;  fi'oni  Ions  No.  8  wa.s  collimated  and  also  was  nearly  rect- 
ani;:ulai'  sliapc.  Tlie  lon^  side  of  the  10.  6  um  beam  cross-section  was  in 
the  direction  parailei  to  the  piane  of  paper  aiKl  was  about  8  mm.  The  short 
side  dinien.sion  was  about  2  mm.  Hie  reason  for  shapin(;  both  the  1.06  ^im 
pump  and  th«'  10.  G  pm  beam  into  rectangular  shape  is  to  utilize  available 
1. 06  pm  pump  power  and  10.  6  pm  illumination  power  efficiently.  Hie 
resolution  test  bar  cliart  shown  in  Figure  24  was  placed  between  lens  No.  8 
and  lens  No.  9.  Tlie  test  chart  was  oriented  such  that  horizontal  bars  are 
in  the  direction  perpendicular  to  the  paper. 

The  resolution  chart  used  is  the  standard  Air  Force  1951  target. 
The  change  in  pattern  size  is  in  a  geometric  progression  based  on  the  six 
root  of  2.  The  number  of  line  per  millimeter  doubles  with  every  sixth 
target  element.  An  element  of  target  contains  two  patterns  of  three  lines 
each  at  right  angles  to  each  other.  The  widtli  of  line  to  width  of  space  ratio 
is  a  nominal  1:1.  The  line  length  to  line  width  ratio  is  5:1.  Table  IV  lists 
the  number  of  line-pair  per  millimeter  of  each  element. 


Table  IV  Number  of  Line -Pair  per  Millimeter  of  Each  Element 


of  1951  Air  Force  Resolution  Target 

Group  Number 

1 

Element  Number 

2  3 

4 

5 

6 

0 

1 

1. 123 

1.  26 

1.415 

1.  589 

1.782 

1 

2 

2.  246 

2.  52 

2.  83 

3.178 

3.  563 

2 

4 

4.492 

5.04 

5.  66 

6.356 

7. 125 

3 

8 

8.984 

10.  08 

11.  32 

12.  712 

14.  25 

4 

16 

17.968 

20. 16 

22.  64 

25.  424 

28.  5 

5 

32 

35.  936 

40.  32 

45.  28 

50.85 

57 

120 
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Tabic  rv  (Continued) 


Element  Number 

Group  NumlM'i-  1  2  3  4  5  6 

6  64  71.872  80.64  00.56  101.7  114 

7  128  143.744  161.2  181.2  203.4  228 

Lens  No.  9,  with  152.  4  mm  focal  length,  was  used  as  an  object  mangi- 
fler.  The  magnified  10. 6  pni  image  of  the  target  patterns  was  formed  in 
the  front  focal  plane  of  lens  No.  10.  The  prousUte  crystal  was  positioned  so 
that  the  back  100.11  plane  of  lens  No.  10  Is  located  at  the  center  of  the  cryt- 
tal.  Lens  No.  11  was  positioned  so  that  the  distances  from  Lens  No.  11  to 
the  proustite  crystal  and  to  the  sampling  slit  were  exactly  equal  to  Its  focal 
length.  The  object  field  entering  the  nonlinear  crystal  is  Fourier  trans¬ 
form  (angular  spectrum)  of  the  real  object  field  formed  across  the  front 
focal  plane  of  lens  No.  10.  Upconverted  image  field  emanating  from  the 
crystal  is  inverse  Fourier  transformed  by  lens  No.  11  and  the  upconverted 
real  image  is  formed  across  the  back  focal  plane  of  lens  No.  11.  The  lens 
arrangement  described  above  allows  upconverslon  to  take  place  In  Fourier 
space  and  introduces  minimum  of  aberrations.  Focal  length  of  lenses  No.  10 
and  No.  11  is  120  mm  and  120.  7  mm,  respectively. 

The  resolution  test  chart  and  lens  No.  9  were  positioned  so  that  the 
10. 6  Mm  image  of  the  bar  pattern  under  study  is  approximately  25  mm. 

This  arrangement  makes  full  usage  of  10.  6  um  object  angular  aperture, 
since  the  maximum  wave  front  tilt  of  angular  spectrum  of  object  field  at  the 
crystal  is  approximately  sk  degrees  which  is  comparable  with  measured 
half  angular  .iperture  of  6. 15  degrees. 
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V - 1 .  Intensity  Profile  of  0.967  um  Images 

Measured  intensity  profile  of  Upconverted  images  arc  shown  in 
Figures  25,  26,  27,  and  28.  In  this  measurement  an  arrangement  was 
made  so  that  only  those  test  bars  oriented  perpendicular  to  the  plane  of 
the  paper  were  illuminated  by  the  10. 6  um  beam.  Figure  25  and  Figure  26 
are  the  image  of  resolution  target  elements  1/1  through  1/6.  Figure  26 
was  obtained  with  more  10.  6  um  illumination  on  the  element  1/6.  These 
figures  clearly  show  that  the  upconverter  is  capable  of  resolving  the  target 
element  1/6.  Figure  27  is  the  image  intensity  profile  of  the  elements  1/4, 
1/5,  1/6,  iUid  0/6.  Figure  28  is  the  image  of  the  elements  0/4,  0/5,  and 
0/6.  Salient  parameters  of  10.  6  um  object  and  0. 967  um  image  of  these 
figures  arc  tabulated  in  Table  V. 


Table  V  Parameters  for  Figures  No.  25,  26,  27,  and  28 


Fig. 

No. 

Total 

Width  of 
Targets 
Illumi¬ 
nated 
in  mm 

Magnifi¬ 
cation 
of  Lens 
No.  9 

Total  Width 
of  10. 6  um 

Image  Across 
the  Front 

Focal  Plane 
of  Lens 

No.  10 
in  mm 

Measured 

Total 

Width  of 

0. 967  um 
Image 
in  mm 

Theoreti¬ 
cally 
Computed 
Width  of 
0.967  um 
Image 
in  mm 

25 

8 

3.18 

25.4 

2.3 

2.32 

26 

6.8 

2.  76 

18.  75 

1.71 

1.72 

27 

7.  5 

2.7 

20.3 

2.13 

1.86 

28 

6 

3.2 

19.2 

1.73 

1.76 

Numerical  values  listed  in  last  column  of  Table  V  were  computed  by 
the  relation: 
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w 

b 


w. 

1 


n.  f., 

s  H 

"i  'll' 


wIk'I'c*  Mg  is  the  »n;i(inifitali(tn  factor  of  leiis  No.  9  and  and  are  the 
focal  l(‘n^th  of  l('ns  No.  10  and  lens  N»).  11.  respectively.  A^jreement  be¬ 
tween  theoretically  computed  total  image  width  and  measured  results  is 
reasonable.  Ol)servation  of  these  figures  shows,  however,  that  tlie  image 
width  of  the  individual  bar  at  the  outer  edge  tends  to  be  narrower  than  the 
e.xpected  width.  'FIjc  intensity  reduction  at  the  edge  is  due  to  tapered  10.  6 
pm  illumination  and  |)hase-mismatch. 

The  number  of  linear  resolution  elements  contained  in  the  0.067  um 
image  is  70  if  one  t:ikes  the  element  1/^6  as  the  smallest  resolvable  element. 
The  angular  resolution  corrcspondii.g  to  above  linear  resolution  is  3.0 
miliiradiaus  (Fig.  26). 

The  intensity  profile  of  defocussed  0.  967  um  image  is  shown  in 
Figure  29.  Relative  positions  of  the  proustite  crystal,  lens  No.  10,  lens 
No.  11,  and  the  sampling  slit  were  arranged  such  that  they  are  at  the  lens 
conjugate  planes.  Lens  No.  9  was  removed  from  the  setup.  TTie  resolution 
targets  1  ^1  through  1  were  illuminated  with  10.  6  um  laser  beam  and  were 
positioned  92  cm  from  lens  No.  10. 

The  effects  of  multi-mode  pump  were  also  investigated.  Relative 
locations  of  test  targets,  the  lenses,  the  crystal,  and  the  sampling  slit  were 
the  same  as  that  ised  to  obtain  Fig-ure  25.  The  1.06  um  Nd:YAG  laser  output 
beam  was  purposely  made  multi-mode.  The  pump  beam  was  collimab  and 
shaped  as  described  before.  The  intensity  profile  of  0.  967  um  image  of  the 
target  element  1/1  through  1/6  is  shown  in  Figure  30.  The  resolving 
power  of  upconverted  image  was  poorer  than  TEM  mode  pump  case  as  was 
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Ill  III  III  III  Ml 


FIGURE  25.  INTENSITY  PROFILE  OF  0.967  /urn  IMAGE  OF 
10.6  Mm  ILLUMINATED  RESOLUTION  TARGET 
ELEMENTS  1/1  THROUGH  1/6  . 
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FIGURE 


86  f/m 
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FIGURE  27.  INTENSITY  PROFILE  OF  0.967  IMAGE  OF  10  6 
ILLUMINATED  RESOLUTION  TARGET  ELEMENTS  1 
1  5,  16  AND  0  1 


INFAREO  OBJECT 

(0/6  TARGET  WIDTH  =  0,28  mm> 

RESOLVING  POWER  TEST  TARGE 
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FIGURE  29.  INTENSITY  PROFILE  OF  0.967  DE FOCUSSED 

IMAGE  OF  10.6  ILLUMINATED  RESOLUTION 
TARGET  ELEMENTS  1/1  THROUGH  1/6 
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FIGURE  30. 


expected.  Quantitative  evaluation  of  resolution  degradation  due  to  pump 
beam  divergence  was  not  attempted. 

V-2.  Recording  of  0.967  um  Image  of  10.  6  um  Illuminated 
Resolution  Test  Chart  on  Photographic  Films 

The  experimental  arrangement  described  in  Section  V-1  was  employed 

except  lenses  No.  3,  No.  4,  No.  7  and  No.  8  were  deleted.  Cioss-scction 

of  both  pump  beam  and  10.  6  um  illuminating  beam  is  circular.  A  lens  with 

500  mm  focal  length  was  used  as  lens  No.  11.  This  is  required  in  order  to 

obtain  reasonable  size  0.  967  um  image.  Kodak  infrared  spectroscopic 

films  Type  I-M  were  used  to  record  the  0.967  um  images.  A  combination 

of  0.967  um  bandpass  and  1.06  um  rejection  filters  is  placed  in  front  of  the 

film.  Figure  31  is  an  enlarged  picture  of  the  recorded  0.967  um  image  u! 

resolution  target  elements  1/1  through  1/6.  The  density  at  the  edge  of  the 

picture  is  low,  but  it  clearly  shows  that  the  target  element  1/6  is  resolved. 

The  actual  size  of  the  image  recorded  is  6.  2  mm  long  and  3.  2  mm  wide. 

The  measured  magnification  (image  size /target  size)  of  1.14  agrees  with 

the  expected  magnification  1. 2  calculated  by  the  relation: 

n  f  . 

M  =  M-  fl  —  —  =  3. 18  X  0. 091  X  7^  =1.2 
9  n,  120 

The  reduction  of  density  at  the  picture  edge  is  again  due  to 
illumination  taper  and  phase -mismatch.  A  portion  of  elements  1/1  and  1/6 
were  missing  from  the  picture.  It  was  checked  later  th  U  the  10.  6  pm  beam 
was  not  illuminating  those  portions  of  the  pattern.  The  contrast  of  the 
picture  is  poor  because  of  long  exposure  (180  sec. )  required  and  some  back¬ 
ground  light  leakage. 
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FIGURE  31.  0.967  am  IMAGE  OF  10.6  urn  ILLUMINATED 
RESOLUTION  TEST  TARGET 
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VI.  IMAGING  UPCONVERTER  SENSITIVITY 

As  a  result  of  analyses  on  upconversion  carried  out  in  Chaptc'r  4, 
parametric  relationships  have  been  developed  which  are  used  in  this  stction 
for  evaluating  the  iX)tentiai  sensitivity  of  present  and  future  upcoiiverters. 
The  equations  used  in  the  calculations  arc  first  giviui.  The  results  are  then 
presented  and  discussed. 


A.  Upconversion  Efficiency 


A  key  iiaraineter  in  determining  sensitivity  is  the  iMiwer  upcon¬ 
version  efficiency  which  for  small  crystal  thickness  is  given  by  the  classical 
expression 


13.04  I  Pp 


n  .  n  n  X 
1  s  p  s 


(1) 


where  =  effective  nonlinear  coefficient  for  upconverter  material 

n  i,  n  n  =  refractive  indices  for  input  signal,  output  signal  and  pump 
respectively 

f  =  thickness  of  upconverter  material 
X.,  X  -  wavelengths  of  input  signal  and  output  signal  resix;ctivcly 

1  o 

Pp  =  pump  power 

A  =  interaction  area  on  upconverter  material 


For  a  Type  II  matched  proustUe  crystal,  a  Nd;YAG  pump  laser,  and  a  lO.Gum 
input  signal  the  following  parameter  values  hold. 

0  =12  degrees  -  0.  2  radians 
X  =  10.6  um 

Id  ,,1  =  66  X  10"^  esu  based  on  d^n  -  75  x  10"9  esu  and  matching 
angle  of  20° 
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'  0.  6  cm 

n.  2.677 

w  2.804 

s 

n  2.817 

0.  967  uin 


Fi)r  these  parameters  we  compute 

-12  2  -1 

K  2.95  10  m  W 


B.  Number  of  Resolution  Elements 


It  can  be  shown  that  (or  diffraction  limited  imaging  upcon version 
in  both  a  Fourier  space  device  and  an  image  space  device  an  imaging  upcon - 
version  figure  of  merit  may  be  defined  for  a  circular  image  using  equation  l 
as 


T)  n  -  O^P  K. 
'u  pi 


(2) 


where 


-  conversion  efficient  of  upconverter 


n 


number  of  Rayleigh  resolution  elements  in  upconverted  image 
input  signal  acceptance  angle 
pump  jxiwer 


4  K  '2  _2 

K  -  ±  -  -  8.07  V  10  ^  u; 

^  "  (1.22  x.r 

X.  wavelength  of  input  signal 


C.  Sensitivity 

The  noise  equivalent  |X)wer  of  the  upconverter  system  (NEP) 

s 

may  be  computed  a.s  o.lhtws. 
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A  10.6  Jim  detection  system  is  employed  in  which  the  upconverter 
is  followed  by  the  Varian  Associates  image  intensifier  being  developed 
under  ONR  sponsorship,  followed  in  turn  by  a  detection  device  with  a  res¬ 
ponse  peaked  to  the  yellow  green  output  of  the  P-20  phosphor  at  the  image 
intensifier  output.  Then 


(NEP)g 


2h  1/ 

s 


T)  T)  - 

'u  'll 


+ 


(NEP)^ 


'u 


where 

h  =  Planck's  constant 

1/  -  frequency  of  output  signal 
s 

-  quantum  efficiency  of  image  intensifier 
=  power  gain  of  image  intensifier 
(NEP)j  =  noise  equivalent  power  of  visible  light  detection  device 

The  first  term  is  due  to  quantum  noise  effects  and  the  second 
term  is  the  noise  contributed  by  the  detector  stage.  The  following  set  of 
values  is  employed: 


100  lvalue  obtained  from  Varian) 

nji^o.io 

2hi/g  -  4  X  10*  watt 


(NEP)^  = 

10 

watts  for  a  1  Hz  bandwidth 

We  obtain 

(NEP)g  = 

1 

r4  X  10*'®  10*^®] 

^u 

1  0.1  ’  10^  J 

1 

IV  5  >10'^** 

(4  +  1 )  X  10  =  - 

^u 

L  J 

(3) 
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Note  that  under  the  above  conditions  quantum  noise  limited  opera¬ 
tion  has  been  approached  and  a  furthe*'  increase  in  image  intensifier  gain 
will  do  little  to  improve  overall  system  performance  (although  an  improve¬ 
ment  in  quantum  efficiency  would  still  be  desirable). 

D.  Results 

We  obtain  for  the  throe  principal  equations: 

rj  -  1.03  •  10'^^  P  /A 
'u  p 

n.r»  3.23  •  10"^  P 
u  p 

(NEP)  -  i— -  watt/Hz 

Note  that  the  first  equation  yields  an  upper  bound  on  upconversion 
efficiency  based  on  crystal  power  handling  capability  (pump  power  density). 
The  second  equation  gives  the  relation  between  number  of  resolution  elements 
and  total  pump  power.  The  third  equation  gives  the  sensitivity. 

Using  the  aliove  values  several  interesting  cases  of  imaging  up- 
conversion  have  been  tabulated  in  Table  I  for  two  materials: 

1.  Proustite  (present) 

2.  Chalcupyrite  (future) 

The  first  one  is  a  fair  approximation  of  the  present  imaging  up- 

converter  nut  including  the  NEP  estimate.  This  NEP  estimate  depends  on 

the  Varian  Associates  image  intensifier  that  is  expected  GFE  on  the  followon 

-8  1/2 

progr  »m.  The  present  system  NEP  is  about  10  w/Hz  since  no  image 

intensifier  is  used,  although  on  a  related  single -resolution  element  10.6  4m 

-9  1/2  -10 

upconverter  a  value  of  1. 1  x  10  watt  sec  measured  (3.7  x  10  watt 

1  ^2 

sec  corrected)  was  obtained,  which  is  the  best  value  reported  to  date. 
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TABLE  VI.  CHARACTERISTICS  OF  IMAGING  IIPCONVERTERS 
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2.25x10^  100  cw  140  CdGeP,  8  2x10"^  6.1x10'^^  Shows  future 

(150x  150)  potential 

sensitivity 


Experience  both  at  AIL  and  elsewhere  has  indicated  that  a 
2 

power  density  near  100  w/cm  the  proustite  crystal  is  damaged.  Therefore 

the  first  two  cases  in  Table  IV  are  safe  and  the  two  shown  for  300  and  140 
2 

w/  cm  are  not  realistic.  From  the  last  equation  above  a  power  density  of 

6  2  -5 

10  w/m  gives  a  maximum  conversion  efficiency  of  1.03  x  10  .  This  rep¬ 

resents  a  material  limit  which  can  be  approached  but  not  safely  surpassed 
on  a  cw  basis.  A  pulsed  upconverter  with  significantly  high  peak  power  but 
lower  avera'  e  power  is  capable  of  higher  values  of  conversion  efficiency. 

For  example  a  1  megawatt  pulse  on  a  150  x  150  upconverter  may  be  capable 
of  achieving  up  to  1.  4  x  10  ^  conversion  efficiency. 

The  Q-switched  case  in  Table  VI  shows  a  very  great  improve¬ 
ment  over  the  cw  cases.  However,  this  improvement  is  real  only  during  the 
pulse  duration.  In  order  to  obtain  effectively  an  average  improvement  in 
system  performance  the  system  noise  during  the  interpulse  period  must  be 
suppressed  by  a  gating  technique  that  isolates  the  output  detector  from  the 
system  nclse  source.  Such  noise  sources  include  dark  current  from  the 
image  intensifier  and  laser  lamp  leakage.  A  gating  pulse  on  the  power  supply 
for  the  imsge  intensifier  or  an  electro-optic  shutter  would  serve  this  purpose. 
Also,  the  output  detector  must  provide  integration  of  the  upconverted  signal 
during  the  pulse.  This  sample  and  hold  approach  should  result  in  significant 
overall  system  performance  although  the  maximum  usable  information  band¬ 
width  has  been  reduced  to  the  periodic  sampling  frequency.  In  many  applica¬ 
tions  a  100  frames/second  imaging  device  would  nevertheless  be  extremely 
useful. 

An  important  aspect  of  imaging  upconverter  development  is 
its  future  potential  for  growth.  At  several  U.  S.  A.  laboratoriei^  and  in  the 
U.  S.  S.  R.  work  is  proceeding  to  develop  better  upconverter  materials. 
Significantly,  in  the  U.S.  A.  Standford  University  and  Bell  Telephone 
Laboratories  are  investigating  the  chalcopyrite  crystal  group  for  nonlinear 
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optics  application.  CdGeAs  has  shown  good  results  for  second  harmonic 
generation  and  CdGePg  has  a  calculated  figure  of  merit  57  times  greater 
than  proustite  for  a  Nd:YAG  pump  and  a  10.6  /urn  signal.  The  power -handling 
capacity  of  this  material  is  comparable  to  GaAs  and  therefore  very  much 
better  than  proustite.  Attempts  to  grow  optical  quality  usable  size  crystals 
of  CdGePg  are  being  made  at  Stanford  University.  They  predict  a  usable 
sample  within  one  year. 

The  computations  of  Table  VI  are  shown  for  proustite  and 
are  then  extrapolated  to  include  this  new  material  with  which  the  long  term 
potential  of  imaging  upconverters  becomes  additionally  attractive. 
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APPENDIX  A 


FOURIER  REPRESENTATION  OF  COMPLEX  FIELD  DISTRIBUTION 
If  the  complex  field  distribution  across  any  plane  is  Fourier 
analyzed,  the  various  spatial  Fourier  components  can  be  identified  as  plane 
waves  traveling  in  different  directions.  The  field  at  any  other  point  in  space 
can  be  calculated  by  summing  the  contribution  cf  these  plane  waves,  taking 
due  account  of  the  amplitude  change  and  the  phase  shifts  they  have  undergone 
in  propagating  to  the  point  in  question. 

Let  the  field  across  a  plane  at  z  =  Zy  be  represented  by  1^^ 

(x,  y,  Zj).  This  field  can  be  represented  as  a  Fourier  integral 

Uj(x,  y,  Zj)=  / z^)  e^”  ^  dx  (Al) 

The  component  field  amplitude,  called  the  angular  spectrum  of  y,  z^), 

IL  expressed  as  inverse  transform 

1 

=(2^)2  /  -dx  (A2) 

Now  recall  that  the  equation  for  a  unit -amplitude  plane  wave  propagating  with 

2  XT 

direction  cosine  is  simply  exp  [i - (ax  +  fty  +  yz)]  where 

2  2  2  ^ 
a  +  /3  +  y  =1.  Thus  across  the  plane  z  =  z^,  a  complex  exiDonential  func¬ 
tion  exp  +  x^y  +  rz)  \  may  be  regarded  as  a  plane  wave  propagating 

with  direction  cosines 

“  '  ''x-  ^  -y  -  02 

Consider  now  the  angular  spectrum  of  U  across  a  plane  parallel  to  the  xy 

plane  at  a  distance  z  from  the  transverse  plane  z  .  Let  the‘li  _(x,  z  ) 

^  1  ^  ^  2 

be  the  arigular  spectrum  of 
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-5*5.'  's'  =  (2„)2  ^  "s'E'  's' 


(A3) 


where 


K  =  K  X  +  K  y 
X  o  y  o 


p  =  X  X  +  y  y 
-00 


The  effect  of  wave  propagation  on  the  angular  spectrum  can  be  found  easily  by 
substituting  Fourier  integral  representation  of  UgCp,  2^  into  the  wave  equa¬ 
tion.  For  the  case  of  source  free  isotropic  space  it  can  easily  be  shown 
that'UglK,  Zg)  and‘Uj(K,  z^)  are  related  by 


..^2  2.1/2.  . 


(A4) 


2  2 

This  result  shows  that  when  k  satisfies  K  -  k  >0  the  effect  of  propagation 
over  a  distance  (Zg  -  z^)  is  simply  a  change  in  the  relative  phases  of  the 
various  angular  spectrums.  Since  each  plane-wave  component  propagates 
at  a  different  angle,  each  travels  a  different  distance  to  reach  a  given  trans¬ 
verse  plane  and  relative  phase  delays  are  thus  introduced. 

We  note  that  the  field  U2(a  2.^  can  be  written  in  terms  of  initial 
angular  spectrum  by 

.  ,„2  2.1/2  (z-  -  z J  . 

U  2^p,  Zg)  ^  /  ‘Uj(k,  Zj)  e^  -  e^-  ^dx  (A5) 
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APPENDIX  B 


WAVE  EQUATION  FOR  ANGULAR  SPECTRUM 
In  this  section  the  wave  equation  that  describes  interaction  of 
angular  spectrum  in  a  lossless  non-linear  material  is  derived.  In  aniso¬ 
tropic  medium  the  electromagnetic  fields  are  related  by  Maxwell's  equation: 


V  X  E  (p,z,  t)  = 


1  SB  (p,  z,  t) 
c  9t 


(Bl) 


apje,  Z,  0 

^  X  H  (p,  z,  1)  =  -  »t  (B2) 

where  p  is  transverse  coordinates  normal  to  the  z  direction  and  c  is  velocity 
of  light.  Other  symbols  bear  customary  meaning.  If  the  fields  are  mono¬ 
chromatic  with  angular  frequency  a)  and  its  time  dependence  can  be  expressed 
in  complex  notation  exp(-iu>t).  Spatial  variation  of  field  .'ctors  can  be  rep¬ 
resented  by  two  dimensional  Fourier  integral: 

U  (p,  z)  =  f  %  (k,  z)  e^ dp 

where  ^  ^  ^  ^  vector  spatial  angular  frequency.  Bold  face 

^(p,  z)  represents  field  quantity  in  space  and  italic  z)  represents  its 
Fourier  component. 

8 

In  such  a  field,  differentiation  in  time  —  is  always  equivalent  to 
multiplication  -  ico  and  Eq.  (Bl)  can  be  written  as 

^  X  (  /  5  («,  z)  -dK)  -  f  ®  (k,  z)  e^-  ^  d  k 
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Since  the  inte(;r;itiun  is  in  k  space,  the  order  of  spatial  differentiation  and 
t:  space  integration  can  be  interchanged. 


j  I  ^  5  z)  *  i  X  S{k,  z)  e^-  -  dK 

i  ,  i  u’  ,  .  i'^  •  P  . 

=  —  ®  (k,  z)  e  -  d  K 

'  c  —  -  - 

MiLwvell's  equation  for  component  fields  in  nonma^tnetic  material  (  /^  =  1)  then 
becomes 


^  5  (k,  z)  >  i  K  X  5  (k,  z)  =  —  3C  (k,  z)  (B3) 

-  -  -  c  —  - 

i  ti.' 

V  X  ^  (K.  z)  +  i  K  X  K  (k,  z)  --  -  —  »  {k,  z)  (B4) 

combinini;  Eq.  (B3)  and  Eq.  (B4). 

V  X  .TC  (»{.  z)  =  -  (  K  X  V  X  5(k,  z)  +  iKX  KX  5  (k,  z)]  (B5) 

Taking  curl  of  Eq.  (3B)  and  substituting  it  in  Eq.  (B5). 

2 

VxVx^(k.  Z)=  ~  0(K,Z)  ^  KX  KX  S  (k,z)  (B6) 

c 


Using 


D(k,  z)  =  c  :  fi  (k,  z)  +  4  (3  z) 

and  defining  equivalent  dielectric  constant  tensor 


c  =  T  +  >c»c 

=»  0  2  - 
c 


(B7) 


where  €  ^  linear  dielectric  constant  tensor 
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^{k,  z)  =  Fourier  components  of  polarization 
and 

I  =  idemfactor  or  unit  dyadic 
Wave  equation  for  an^ilar  spectrum  becomes 


V  XV  X  s  K  z)  -  <  g  :  5  K  z)  = 


4iru> 

2 

c 


^  K  z) 


(B9) 


149/150 


z)  =  Fourier  components  of  polarization 
and 

I  =  idemf actor  or  unit  dyadic 
Wave  equation  for  an^plp.r  spectrum  becomes 

4ir  ^ 

V  X  V  X  6  (k,  z)  -  €  :  z)  =  — (p  {k,  z) 

— e-—  z 

c 


(B9) 


149/150 


APPENDIX  C 

FOURIER  INTEGRAL  REPRESENTATION  OF  NONLINEAR  POLARIZATION 

The  nonlinear  polarization  of  sum-frequency  produced  in  the 

nonlinear  material  is  given  by 

P(/>,  z)  =  d:E  (p,  z)  E.  (p.  z)  (Cl) 

-  -P  -  - 

where  p  =  transverse  position  in  the  material 
z  =  longitordinal  position  in  the  crystal 
Ep(p,  z)  =  electrified  of  pump  beam 

Ej  (p,  z)  =  electric  field  of  infra-red  object  wave 

d  =  second  order  polarization  tensor 

Expressing  field  quantities  in  their  Fourier  integral  representa¬ 
tion  Eq.  (Cl)  becomes 

P(g,Z)  =  d: //  Ep  (ijyz)  ^^2) 

Using  phase -matching  condition 

K  =  K  +  K,  (C3) 

-s  -p  -i 

We  write  Eq.  (C2)  as 

P(p,z)  =  /  d  e*-8‘^{d  ;[  /  dKj  £j  («j»d)]} 

-  f  dx  e*-®  -  (d;  5  (k  ,  z)  ®  5  («  ,  z)  ]  (C4) 

— 8  =  P  “8  -i  "O 

where  5  j  (k,  z)  o  £p  (^»  z)  indicate  convolution  integral. 
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It  nuw  follows  Eq.  (C4)  that  Fourier  components  of  the  nonlineai  polariza¬ 
tion  at  si<ni-frequency  are  product  of  second-order  polarization  tnesor  and 
the  convolution  of  those  of  the  input  object  and  the  pump  angular  spectrum 
that  satisfy  the  condition  Eq.  (C3) 

(P(k  .  z)  =  d  :5  .  (k  ,  z)  »  <5  (k  ,  z)  (C5) 

-  “S  =  -  i  — s  -  p  — s 
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APPENDIX  D 
MISMATCH  FACTOR 


We  shall  carry  out  the  analysis  in  the  index  space  rather  than 
in  the  wave  vector  space.  These  two  spaces  are  related  only  by  a  constant 
factor  u)/c  in  a  way 


.  ^  A  Cl) 

n(K)  =  -  K 
0  c  - 

where  is  a  unit  vector  along  the  direction  of  K .  The  constant  factor 
ct>/c,  therefore,  does  not  appear  in  the  index  space. 

Referring  to  Fig.  12  of  the  text  the  equation  of  indicatrix  of  a 
uniaxial  crystal  expressed  in  the  Cartesian  coordinate  system,  in  whi>:h 
the  z  axis  makes  an  angle  a  with  the  crystallographic  axis  Z  and  the  x 
axis  lies  in  the  crystallographic  X  Z  plane,  is  given  by 


2 

X 


2 

cos  a 

2 

n 

e 


,  2 
sin  a 


n 


2 

+  — 2'-+-  2xz  cosot  sina 
n 


l(Dl) 


where 


n 


0 


=  index  of  refraction  along  the  optic  axis 


n  =  index  of  refraction  normal  to  the  optic  axis 


n  =  index  of  refraction  along  r  =  xSt^  ^  ^0  ^  ^^0 


Equation  of  tangent  plane  at  the  point  z  =  n,  x  =  y  =  0is 


z  H  n  cosa  sin  a 


(D2) 


If  we  define  <  to  be  an  angle  that  the  normal  to  the  tangent  plane 
(D2)  makes  with  the  z-axis,  the  equation  of  tangent  plane  can  be  written  as 
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z+  X  lull  ‘  -  11  -  G 


Cc!!!''.! riii-.;  I  'lS,  {D2/ und  (D3) 


2  .  /  1  1 

tan  ;  -  11  coSiv  sill  o  (  «  '  "o 

C  M 

11  11- 

e  0 


(D3) 

(D4) 


RCiviulion  d!  iiidicatrix  in  the  coordinate  (x,y,  z)  expressed  in  terms 


of  11-,  n  .  11.  a; id  I'.ec times 
0  e 


The  difference  between  an  arbitrary  vector  (i?  ,  V  ,  ri  ) 

X  y  z 

and  a  vector  {n  ,  V  ,  )  lying  on  the  indicatric  surface  is 

X  y  z  •  “  2 

.  ^  “r  ^  / cos^o  sina  tan^A  y  ,  , 

"z-"2  "z  ■  2t"xl  ~  -  —2  r  “21 

\  n  n.  11  /  n 

e  0  e 

Terms  containing  higher  than  second  power  in  r?  and  t?  are 

X  y 

neglected  in  arriving  at  Eq.  (D6).  The  assumption  is  valid  for  paraxial 
case. 

Eq,  (D6)  is  equal  to  the  mismatch  factor  defined  in  Eq.  (27) 
Chapter  3  in  the  text 

The  relation 


2 

cos  Q 


11 


2  2 
sin  o  tan  C. 

2  ■  2 
n»  n 


n 


2  2 
n  n„ 
e  0 


(D7) 


e  0 

can  easily  be  proven  by  using  Eq.  (D4)  and  the  relation 
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n 

0  e 


.2.2  2  2  .1/2 

(n.  sin  a  +  n  cos  a) 

0  e 

Using  Eq.  (D7)  we  finally  obtain  the  mismatch  factor 


r  2  2  2 

n  TJ  r) 

.  n  X  y 

ih  =  V  -TJ  =v  -  n  +  V  tan£  +  -  — ^ 

^zzz  X  2  22  2 

n  n  n 

^  e  0  e 
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n^n 
0  e 


,2.2  2  2  ,1/2 

(n^  sin  Qf  +  n  cos  a) 

0  e 


Using  Eq.  (D7)  we  finally  obtain  the  mismatch  factor 


il,  =  V  -  n  =  V  -  n  +  r/  tant  +  - 
^  z  z  z  ^  ^ 


r  2  2 

n  7] 

n  X 


n  n.  n 
e  0  e 
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APPENDIX  E 


HANKEL  TRANSFORM  OF 


sin(ax^) 


dX 


Hankel  transform  of  ^  is  defined 


ax^ 


H(p)  =  ^  ^  J  (xp)xdx 


ax 


Changing  the  variable  ax"^  =  y,  it  can  be  written 


(E-1) 


(E-2) 


Integration  by  parts  gives 


H(P)='^^372  -^0  J()(p7  a  ^y’^^^Si(y)dy 


(E-3) 


where  Si(y)  is  a  sine-integral  function  defined  by 


and 


„ ,  .  ry  sin  t  ...  TT 

si(y)  =  J  0  ~Y~  ■^2 


sMy)  =  -  L 


sin  t 

'y  t 
Using  the  relation 

Jq(x)  =  -  Jj(x) 

We  can  write 


dt 


(E-4) 


(E-5) 
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H(p) 


4a 


3/2  -"o 


(E6) 


=  -^372  fo  •'l  i  -^0  ’’l 

where  change  of  variable  z  =  has  been  used  for  the  2nd  term. 


Using  the  formula 


and 


/  y  si(y)  (2>/by  )  dy  =  -  [si(b)  +  1^]  /  \/b 


Iq  Jj(^^)dz  =  l 

We  finally  obtain 


(E7) 
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APPENDIX  F 

EVALUATION  OF  MISMATCH  FACTOR  FOR 
NEARLY  ISOTROPIC  CRYSTALS  U  =  0) 


Refering  to  Fig.  FI,  the  mismatch  factor  AK  is  given  by 


=  '^sz  -  Km 

/ 

(FI) 

K  =  K  cos  0  +  K.  cos  0 
sz  p  pi  i 

(F2) 

1  '2  '2 
~  K  +  K.  -  (K  0  +  KH  1 

P  1  2  p  p  ip 

From  geometry  (Fig.  3  Chapter  2) 

0  =  -T— TrT«=  o'  n 
s  s  +D/2  S  s 

(F3) 

0  =  =  0  'n 

p  p  -  D/2  p  p 

(F4) 

0.  =  — =  o'  n 

i  s.D/2  1  1 

(F5) 

Eliminating  y  from  (F3)  and  (F4) 

n  (s'  +  D/2) 

P  ~  "p  "D/2)  s 

(F6) 

Eliminating  y  from  (F3)  and  (F5) 


Usini'  0  K  /K  ,  (F6)  and  (F7).  We  can  rewrite  (F2) 
s  s  s 


Ksz  =  Kj,  .  Ki  - 


8  (S  ♦  2  ’  '^s  2  ,  p  ,  2 

- Z -  n  (p-i)** 

2  D  2 

2Kf 


D  .2 


nf(8  -  «  ) 


(F-8) 


For  the  collinear  matched  nearly  isotropic  case 


Kp  .  Ki  =  Km-  Kg 
Finally  we  obtain 

2  2,  »  D, 

-  _^sns(8 


^  r  Ki  •‘i  ] 


(F-9) 


If  p  = 


s  /  "s 


AK  =*- 


2Kj  Vnjfs-^) 


Ki 

K 


M> 


2K. 


(F-IO) 


Since 


n  (8 '  +  ^ ) 
s  2 


1 
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I 


HGURE  FI.  WAVE  VECTOR  RELATION 
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